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The purpose of this work was to examine the 8-ray spectra of Mn* and. I"* as seen in the 
Wilson cloud chamber in an effort to correlate the measurements with the existing empirically 
regarded relations of the Fermi or the Konopinski-Uhlenbeck theory. A further study of the 
y-rays which accompany these 8-rays was made. In all measurements the difficulties of cloud- 
chamber technique were critically investigated. The results appear to fit the Fermi rather 
better than the K-U theory. However, the errors of cloud-chamber measurement are seen 


to be very serious in this type of work. 


HE £-ray spectrum of Mn* was first meas- 
ured in the cloud chamber by Gaerttner, 
Turin and Crane.' Using a magnetic field of 850 
gauss, they found the spectrum to consist of a 
single K-U group having an extrapolated end 
point at 6.5 mc? (2.8 Mev). This spectrum was 
later measured by Brown and Mitchell.? The 
fields in their cloud chamber were 375 and 490 
gauss. They found that the K-U plot resolves 
itself into two groups having extrapolated end 
points at 6.8 mc? (3.0 Mev) and 3.4 mc? (1.2 Mev). 
The same spectrum was also measured in a 
B-ray spectrograph by Alichanow, Alichanian, 
and Dzelepow.* They reported the end point to 
be at 7.3 mc? (3.2 Mev), which does not agree 
with the results given above. 
It seemed to us that we might be able to 


* An abridgment of the dissertation submitted by R. H- 
B. in partial fulfillment of the requirements for the Ph.D. 
degree at New York University. 

Tt Now at the Frankford Arsenal, Philadelphia, Pennsyl- 
vania. 

1 Gaerttner, Turin and Crane, he ee Rev. 49, 793 (1936). 
as 30) V. Brown and A. C. G. Mitchell, Phys. Rev. 50, 593 

3 Alichanow, Alichanian and Dzelepow, Nature 136, 257 
(1935); Physik. Zeits. Sowjetunion 10, 78 (1936). 


- remove the uncertainties in the previous work 


by means of additional experiments. As a result 
of our work, we are able to report that the 
upper limit of this spectrum is close to the 
point found by Gaerttner, Turin and Crane ;' the 
second group announced by Brown and Mitchell 
actually does appear to be part of the spectrum ; 
in addition, we find evidence of a third group 
if the K-U theory be applied. 

The y-rays from Mn** have been investigated 
by Mitchell and Langer,‘ using coincidence 
counters to determine the energies of the 
Compton electrons ejected from an aluminum 
plate. They obtained 1.65 Mev as the y-ray 
energy, and claimed that their work shows the 
y-ray to be monochromatic. Later, Livingood 
and Seaborg® measured the energy of this y-ray 
by finding the thickness of lead necessary to 
reduce its intensity to one-half. The energy thus 
obtained was 1.2 Mev. In the spring of 1939, 


4A.G.G. Mitchell and L. M. Langer, Phys. Rev. 52, 137 
(1937). 

( he Livingood and G. T. Seaborg, Phys. Rev. 54, 391 
1938). 
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Fic. 1. The effect of the radius of curvature of the tracks 
upon the observed shape of the spectrum. If a track is to be 
measured, it must pass inside the circle of radius B, or, at 
least, become tangent to it. As shown in the lower left 
diagram, one condition for this tangency is given by 


B)? =p?+A*—2pA cos 
COS = (A?— B*—2pB)/2pA. 


From the diagram at the right, it is seen that for tracks 
whose radius of curvature is less than (A +B)/2, dmin=0; 
for tracks whose radius of curvature is greater than 
(4+B)/2, we have 


(o— B)* =p? + A*—2pA Cos omin, 
COS = (A*+B*+2pB)/2pA. 


When p is less than (A — B)/2, the track cannot enter the 
visible portion of the chamber at all. 


Norling® performed a coincidence experiment 
from which he concluded that there are two 
y-ray lines, one of which follows the most 
energetic 8-ray transition. 

We have measured this y-ray spectrum by 
examining the Compton electrons ejected from a 
thin sheet of Bakelite placed along a diameter of 
the cloud chamber. The distribution of the 
observed electrons indicates y-ray lines at about 
0.6 or 0.7 Mev and at about 1.7 Mev. There is 
also some indication of other lines not clearly 
resolved. 

As this work was being completed, Dunworth’ 
reported another measurement based on coin- 
cidence counting. From his work he concluded 
that there are two lines in the y-ray spectrum: 
one at about 0.6 Mev and the other at about 
1.7 Mev. The relative intensities given are in 
agreement with our results. It should be noted 
that these y-ray line intensities are not in accord 
with the relative population of the two 6-ray 
groups as determined either by Brown and 


Mitchell? or by ourselves. Two months later, 


*F. Norling, Naturwiss. 27, 432 (1939). 
7J. V. Dunworth, Nature 143, 1065 (1939). 


Langer, Mitchell and McDaniel® reported a 
8—vy coincidence counter measurement from 
which it may be inferred that there are at least 
two lines in this spectrum, but nothing is stated 
about their energies. A later measurement is that 
of Curran, Dee and Strothers.* They determined 
the distribution of the Compton electrons ejected 
from an aluminum window in a modified 8-ray 
spectrograph. Their results indicate y-ray lines 
at 0.91 and 2.03 Mev, the less energetic y-ray 
being twice as intense as the harder. 

The end point of the B-ray spectrum of I'*8 
was first determined by Alichanow, Alichanian 
and Dzelepow.* Our measurements agree with 
theirs as to the position of this end point. From 
the shape of the spectrum we find two groups 
having upper limits at 1.05 and 2.10 Mev. 
However, we were unable to find any y-ray of 
sufficient intensity to account for two groups in 
this spectrum. Roberts and Irvine!’ also came to 
the conclusion that if there be any y-ray accom- 
panying this radioactivity, there must be less 
than one photon for every ten electrons observed. 
Later, a y-ray of about 0.4 Mev was reported by 
Livingood and Seaborg," but neither the energy 
nor the intensity agrees with the 8-ray spectrum. 
Further work was done by Tape.” His observa- 
tions on the y-ray agree with ours; those on 
8-particles do not. 


APPARATUS 


The cloud chamber consists of an aluminum 
piston, 7 inches in diameter and 6 inches high, 
moving in a brass cylinder. The expansion ratio 
is controlled by screwing the cylinder into or 
out of the brass cup which forms the base of the 
chamber. The joint between cylinder and cup is 
sealed with Apiezon Q. To the top of the cylinder 
is fitted a removable head, sealed in place with 
Apiezon Q. 

Toward the end of an expansion, the chamber 
was illuminated by six 1000-watt clear glass 


8’ Langer, Mitchell and McDaniel, Phys. Rev. 56, 422 
_ Dee and Strothers, Proc. Roy. Soc. A174, 546 
ina. Roberts and J. W. Irvine, Jr., Phys. Rev. 53, 609 

J. J. Livingood and G. T. Seaborg, Phys. Rev. 54, 777 
F. Tape, Phys. Rev. 56, 965 (1938). 
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RADIOACTIVITY 


photoflood lamps, made specially for us by The 
Westinghouse Electric and Manufacturing Com- 
pany. The life of these lamps was in excess of 
30,000 expansions. Cylindrical lenses and suitable 
slits were used to collimate a flat beam of light 
about 1 cm high across the chamber. The pictures 
were taken along the lines of force of the mag- 
netic field. The cycle of operations was controlled 
by a commutator whose shaft also operated the 
camera through a pair of universal joints. 

The coils for producing the magnetic field have 
already been described.'* They are more efficient 
than a Helmholz-Gaugain arrangement—only 10 
amperes at 240 volts were required to produce a 
uniform field of 880 gauss throughout the volume 
of the chamber. 


EXPERIMENTAL DETAILS 


The magnetic field was measured in the usual 
way by means of a ballistic galvanometer, 
standard mutual inductance, and suitable ex- 
ploring coils of accurately known constants. The 
error of these measurements was less than 2 
percent throughout the volume of the chamber. 

The radius of curvature of each track was 
determined by reprojecting the photographic 
image to full size on a sheet of paper and match- 
ing the track to the edge of one of a series of 
brass disks. These disk gauges were prepared 
with radii from 4 to 25 cm in half-centimeter 
steps. As each track was measured, an arc was 
drawn in red pencil over its image. This method 
was found helpful in preventing the omission or 
the double tally of a track. 

In the course of this work, it appeared that 
the shape of the 8-ray spectrum was different for 
different fields. An approximate geometric solu- 
tion for this effect is shown in the legend of 
Fig. 1: ¢=the angle which the horizontal com- 
ponent of the electron’s initial momentum makes 
with the horizontal tangent to the emitting 
surface. A =inside radius of the cloud chamber 
= radius of the emitting surface = 9 cm. B=radius 
of the visible portion of the cloud chamber. 
p=radius of the horizontal projection of the 
electron track. 

If an electron is to pass into the visible portion 
of the chamber, it must cross the circle of radius 


13 R. H. Bacon, Rev. Sci. Inst. 7, 423 (1936). 
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B, or at least become tangent to it. The angles 
of emission between which this is possible are 
indicated in Fig. 1. B is a function of the position 
of the point of emission on the surface of radius 
A: this is clear from the upper portion of Fig. 1. 

The probability of an electron path of given 
radius passing inside the circle of radius B is 
proportional to (@max— min), Where the ¢’s are, 
in turn, functions of the position of the point of 


emission. The curve in Fig. 2 shows this theo- 


retical probability plotted against the radius of 
curvature of the tracks. In the same figure is 
shown the value of the ratio: 


number of tracks actually observed 
in a given momentum interval 


number given for that interval by the 
best K-U curve to fit the entire data 


Finally, for the special runs described in the 
next paragraph, the ratio: 


number of tracks observed at 
850 gauss in a given interval 


number to be expected in that 
interval from the 319 gauss curve 


is denoted by the points marked @. 
The apparent shape of the spectrum was 
found to depend also upon the position and size 


Loss of Tracks of Small Radus of Curvature from the Measured Spectrum 
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637 Gauss 0-319 Gauss 
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Radius of Curvature of Tracks cm— 


Fic. 2. Comparison of the quantity (¢max—@min) 
(smooth curve) with the ratio 


observed number of tracks 
number demanded by the best K-U curve’ 
The points denoted by © are the plot of the ratio 
observed number of tracks at 850 gauss 
number to be expected from the 319-gauss curve” 


In the case of I", the conditions applying to the calculation 
of (@max—@min) Were not constant throughout the measure- 
ment of this spectrum, so no computed curve is shown. 
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Mn£K-U Plot 
e 850 Gauss 
> 637 Gauss 
319 Gauss 
2 
1 2 3 a 5 6 4 
Energy of the Emitted Electrons me? 


Fic. 3. K-U plot of the energies of 7200 electrons emitted 
by Mn*, 


of the source. In our early work the Mn*®* was 
deposited on filter papers 3.5 X11 cm which were 
affixed to the glass ring of the chamber top. In 
our later work, filter papers 2X11 cm were used. 
These narrower sources always showed a higher 


proportion of low energy tracks. For comparison 


purposes a series of runs was made with a wide 
and a narrow filter paper on opposite sides of 
the chamber head. Each run consisted of ten 
pictures at 850 gauss followed by ten at 319 
gauss. A total of 2000 tracks obtained in this 
manner were measured. It was found that the 
yield of low energy tracks was the same for 
each source whereas the yield of high energy 
tracks was proportional to the width of the 
source. 


THEORY 


Experiments may be compared with the theory 
of Fermi’ or with that of Konopinski and 
Uhlenbeck" by a method due to Kurie, Richard- 
son and Paxton.'* The theory may be expressed 
in the form 


(1) 


where N=number of electrons having mo- 
mentum 7=Hp/1700, no=momentum of the 
fastest electron emitted, 


2raZn(1+n?)! 
fa = 


144 E. Fermi, Zeits. f. Physik 88, 161 (1934). 
‘ (1988 . Konopinski and G. E. Uhlenbeck, Phys. Rev. 48, 
16 Kurie, Richardson and Paxton, Phys. Rev. 49, 370 


(1936) 


e 850 Gauss ©6637 Gauss 


2 3 4 5 
me? 
Energy of Emitted Electrons 


Fic. 4. Fermi plot for Mn®* of the same data as used in 
Fig. 3. Although the plotted points between 3.4 and 5.1 mc* 
fall quite accurately on a straight line, it is not possible to 
resolve this spectrum further by this method, as the Fermi 
plot of the differential is a smooth curve through which it is 
impossible to pass a significant straight line. 


K=constant of proportionality, depending upon 
the total number of electrons counted in the 
experiments, k=2 in Fermi’s theory, =4 in that 
of Konopinski and Uhlenbeck. 

A plot of (N/f)"* against (1+7?)! should 
result in a straight line, depending on which 
theory fits the experimental facts. A plot of 
(N/f)* against (1+ does result in a straight 
line for many problems, at least in the middle 
energy range. In some cases, however, the 
plotted points do not arrange themselves in a 
straight line, but along a curve that suggests the 
B-rays might be,divided into two or more groups. 

It should be emphasized here that the process 
of decomposing the 8-ray spectrum into groups 
is independent of the validity of the K-U theory. 
This process was first suggested by Ellis and 
Mott.'? They found that the RaC spectrum 
could be built of five components, each having 
the same shape as the RaE spectrum, the latter 
being regarded by them as single. The end 
points and relative populations of these five 
components were correlated with the known 
energies and intensities of the y-rays accom- 
panying the RaC—RaC’ disintegration. 

Now it has been shown by several workers'® 
that, except at the end points, the 8-ray spectrum 


(1938) D. Ellis and N. F. Mott, Proc. Roy. Soc. A141, 502 

18 E. M. Lyman, Phys. Rev. 51, 1 (1937); L. M. Langer 
and M. D. Whitaker, Phys. Rev. 51, 713 (1937); J. S. 
O’Conor, Phys. Rev. 52, 303 (1937). 
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of RaE closely follows the K-U theory. There- 
fore, if the process suggested by Ellis and Mott 
be valid, and if the RaE spectrum be simple, 
then the K-U theory affords a convenient im- 
plement for effecting such a resolution into 
groups, whether the theory itself be valid or not. 
But if the assumption that the RaE spectrum 
is simple be incorrect, then it is not correct to 
interpret departures from a straight line in the 
K-U plot as evidence of group structure. In the 
following procedure nothing is assumed as to the 
theoretical validity of the K-U theory: Equation 
(1) is used merely as a convenient method of 
effecting the process proposed by Ellis and Mott. 


RESULTS 


The work on the 6-ray spectrum of Mn** may 
be divided into two parts: During the first part 
of this work we had but 50 mg of radium mixed 
with beryllium as a neutron source for activating 
the manganese, which was deposited on filter 
papers 3.511 cm. The chamber was filled with 
helium (sometimes with hydrogen) for some of 
the runs, and left filled with air for others. Runs 
were taken at four magnetic fields: 850, 637, 425 
and 319 gauss. The results of the runs with the 
three higher fields, as already announced,'® may 
be summarized as follows: 

The 850-gauss run duplicated almost exactly 
the curve of Gaerttner, Turin and Crane:! The 
observed electrons seemed to fit a single K-U 
group having an upper limit at 6.58+0.10 mc 
(2.84+0.05 Mev). The runs at 637 and 425 
gauss seemed to support the work of Brown and 
Mitchell :* The observed electrons seemed to fit 
into two groups, having upper limits at 6.55 
and 3.25 mc? (2.80 and 1.15 Mev), respectively. 
The end points given by Brown and Mitchell 
are 6.8 and 3.4 mc? (2.9 and 1.2 Mev), re- 
spectively. 

A short run at 319 gauss showed more slow 
electrons than those demanded by these two 
K-U groups. This led us to pursue the matter 
further. 

Meanwhile; we acquired the new neutron 
source containing 200 mg radium mixed with 
beryllium. We now cut our filter papers into 


19 Bacon, Grisewood and van der Merwe, Phys. Rev. 52, 
668 (1937). 


strips 2X11 cm. Runs were taken at 850, 637 
and 319 gauss. As mentioned above, the results 
of this part of the work differed from those of 
the previous summer only in that the relative 
populations of the low energy groups now 
seemed to be greater than they did at first, and 
the great excess of slow electrons observed in 
the weakest field seemed to fit a third K-U 
group having an end point at about 2.2 mc 
(0.6 Mev). See Figs. 3, 4 and 10. 


THE GAMMA-RaAys ACCOMPANYING THE 
Decay or Mn** 


The solution of NaMnQ, was irradiated in the 
regular way overnight and filtered in the morn- 
ing. The still moist filter paper was then rolled 
into as tight a wad as possible (about 5 mm in 
diameter), and placed in a rest between the field 
coils, in the midplane of the cloud chamber. 
Owing to the weakness of this source, as com- 
pared with those used by Richardson and Kurie,”° 
for instance, it was not feasible to place it as far 


425 Gauss 


os ts 2 $ 38 4 
mc —e 


Horizontal Comp of the Mi of Electrons Ejected 
by the Gamma -Rays 


Fic. 5. Observed distribution of 400 electrons ejected 
from Bakelite strip by the y-rays accompanying the decay 
of Mn®*. Upper half: Open circles, > or obtained by 
Richardson and Kurie, Phys. Rev. 50, 999 (1936). Electrons 
ejected from radiator 40 mg/cm? by annihilation radiation 
from N", 250 gauss. Closed circles, electrons ejected from 
radiator 100 mg/cm? by annihilation radiation from Cu 
isotopes excited by bombarding Ni with protons, 425 gauss. 
Lower half: Electrons emitted from radiator 100 mg/cm? by 
gamma-rays accompanying the decay of Mn®, 425 gauss. 


” J. R. Richardson and F. N. D. Kurie, Phys. Rev. 50, 
999 (1936). 
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Beta-Ray Spectrum of 
* 319 Gauss - 710 Tracks 


® 637 Gauss-654 Tracks 


N=Number of Tracks in Each Momentum interval 


Horizontal Components of the Momenta 


Fic. 6. Observed of 1300 electrons emitted 
y 


from the chamber as they did. There is thus 
introduced into the measurement of the angle 
between the observed direction of the initial 
motion of the Compton electron and the assumed 
direction of motion of the incident photon a small 
error of not more than 5° in the worst instance. 

As a check on this work, we measured the 
energy of the annihilation radiation from the Cu 
positron emitters formed by the bombardment 
of a nickel foil by protons in the Columbia 
cyclotron.”! As metallic nickel would disturb the 
magnetic field, the irradiated foil was dissolved 
in HNOs, and the small mass of Ni(NO3)2 was 
laid in the same place as had been the Mn*® 
before. The results of these measurements are 
shown in Fig. 5. For comparison, there is also 
drawn to appropriate scale the result obtained 
by Richardson and Kurie?® who used a lamina 
of about two-fifths the surface density of ours to 
measure the annihilation of radiation from N". 

The results” of our measurements, as shown in 
Fig. 5 seem to indicate a y-ray line at about 600 
or 700 kev, and perhaps one at about 1.7 Mev. 
However, the complete interpretation is not 
clear; it is probable that there are other lines 
also. For the relative intensities, we obtain 


Sum of the intensities of all lines 
below 0.7 Mev 


Sum of the intensities of all lines 
above 0.7 Mev 


21 We take this opportunity to express our appreciation 
to Professor J. R. Dunning and Dr. E. T. Booth for their 
kindness and interest. 

* Bacon, Grisewood and van der Merwe, Phys. Rev. 56, 
1168 (1939); Phys. Rev. 57, 240 (1940). 
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K-U Plot 


° 319 Gauss 


637 Gauss 


. Energy of Emitted Electron 


Fic. 7. K-U plot of the data shown in Fig. 6. 


Fermi Plot 
637 Gauss 
x 

x 

=|~ 
x 
2 3 4 © 
mc? 


Energy of Emitted Electrons 


Fic. 8. Fermi plot of the data shown in Fig. 6. 


These results agree with those of Dunworth;’ 
the rest of our analyses of these data do not. 

These estimates of the y-ray energies, how- 
ever, do not agree with those of Curran, Dee and 
Strothers.? A new measurement was therefore 
made in this laboratory by Titone,* using a 
radiator of only 25 mg/cm’, and taking ad- 
vantage of a total of 600 mg radium mixed with 
beryllium, which was at that time available. 
His results confirm ours. 


DISCUSSION OF RESULTS OF MEASUREMENTS OF 
THE RADIOACTIVITY OF 


The y-rays from Mn*® are quite intense com- 
pared with the total activity. For the arrange- 


23 L. V. Titone, essay for the degree of Master of Science, 
New York University, 1940. 
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ment of filter papers used, about one twenty- 
fifth of the number of electrons measured could 
have been secondaries ejected from the source by 
y-rays. As most of these secondary electrons are 
of low energy, it is probable that a portion of 
the “third group” really consists of secondary 
electrons. Furthermore, the filter papers were 
placed so close to the glass ring that a number 
of electrons must have been reflected back from 
the ring into the visible portion of the chamber. 
Again there is the effect of scattering in the source 
itself. All three effects tend to enhance the low 
energy end of the spectrum. 

It might be argued that the ‘‘third group” 
represents the partial internal conversion of a 
y-ray line at about 600 or 700 kev. This possi- 
bility cannot be definitely excluded, but it 
seems unlikely, as the observed unconverted 
intensity of this line is quite strong already, and 
unless there are several lines of about this 
energy in the spectrum, the possibility of internal 
conversion seems small. 

In the coincidence counting experiments, 
Norling,® Dunworth’ and Langer, Mitchell and 
McDaniel’ find that there is at least one y-ray 
line accompanying the highest energy electron 
group; and Dunworth finds that this y-ray line 
is the one having energy about 600 kev, and 
then supposes that the second 6-ray transition is 
followed by two y-rays, the 1.7 and the 0.6 Mev 
in turn. The findings of Langer, Mitchell and 
McDaniel also suggest this. 

Now, the second group of 6-rays is twice as 
intense as the most energetic group. The high 
energy y-ray, therefore, should, according to 
Dunworth’s picture be two-thirds as intense as 
the low energy y-ray. Dunworth, however, and 
we also, find that the high energy y-ray is only 
about two-fifths as intense as the low energy 
y-ray. 

Thus, in order to correlate the y-ray measure- 
ments with those of the 8-rays, it is necessary to 
postulate more than two #-ray transitions, or 
more than two y-ray transitions, or both. Hence, 
the energy level diagram given by Norling or by 
Dunworth is incomplete. 

Finally we have tried the Fermi plot in Fig. 4. 
Assuming that the distortion above 3 mc? (1 Mev) 
is small, we found by the method of least 


squares the best straight line among the points 
plotted between 3.4 and 5.1 mc* (1.2 and 2.1 
Mev). Having done this, one can postulate two 
groups of energy 5.5 and 3.4 mc® (2.3 and 1.2 
Mev). But the fit is so bad that it is impossible 
to say anything about either relative populations 
or the existence of any lower groups. 


THE RADIOACTIVITY OF [28 


The results* shown in Figs. 6, 7 and 8 may be 
summarized as follows: The observed electrons 
appear to have energies up to about 1.85 Mev. 
The K-U plot apparently shows two groups 
having upper limits at 3.06 and 5.10+0.10 mc 
(1.05 and 2.10+0.05 Mev). The supposition of 
two groups in the 6-ray spectrum of suggests 
that either the groups have different periods, or 
else that y-rays exist whose energy accounts for 
the energy difference between the two groups. 

We measured the period on a Geiger counter 
several times, following one run through twelve 
half-lives (see Fig. 9). The period was found to 
be single to within one minute. The value was 
26+1 minute. We also looked for a y-ray. The 
conclusion was that there was less than one 
photon for every 10 electrons. 


Counting Time in Minutes 


0 30 60 90 120 


Decay of 


2 

Half-Life 25.6 Min. 

100+ 

4 

4 4 
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10 T T T T T T T T T 
180 ~ 210 240 270 300 330 


Minutes after Activation Had Ceased 


Fic. 9. Half-life of 


24 Bacon, Grisewood and van der Merwe, Phys. Rev. 54, 
315 (1938). 
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Mn** \ 


e 850 Gauss 
637 Gauss 
__Best K-U Distribution 
Best Fermi Distribution 
The eleven points between 
_- the arrows fit either curve 
equally well 


1 2 3 4 me 
Comparison of Fermi and K-U Groups 


Fic. 10. The high energy end of the 8-ray spectrum of 
Mn*, drawn to much enlarged scale. It will be seen from 
this figure, that the predictions of the two theories are 
quite indistinguishable over an appreciable range. 


CONCLUSIONS 


Our results appear to lend somewhat greater 
support to the K-U theory than to the original 
Fermi theory of beta-decay. This is in accord 
with the findings of most observers. This ap- 
parent agreement with the K-U theory has been 
ascribed by Konopinski* to the distortion of the 
true spectrum by the scattering of the electrons 
in the B-ray source. In support of this explanation 
is the recent work of Tyler*® and of Lawson,’ 
who find that as the source is made thinner, 
the observed spectrum departs from the K-U 
and begins to approach the Fermi distribution ; 
against this explanation is the work of Richard- 
son and Leigh-Smith,?* who upon measuring the 
B-ray spectrum of ThC as a gas (bismuth tri- 
methyl) in the cloud chamber, find close agree- 
ment with the K-U theory, even to the large 

SE. Phys. Rev. 57, 68 (1940). 

% A.W. Tyler, Phys. Rev. 56, 125 (1939). 

27 J. L. Lawson, Phys. Rev. 56, 131 (1939). 


%H. O. W. Richardson and Alice Leigh-Smith, Proc. 
Roy. Soc. 162, 391 (1937). 
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number of very slow electrons postulated by the 
K-U theory for elements of high atomic number. 

Perhaps it might be mentioned in passing that, 
at the other end of the periodic table, Fowler, 
Delsasso and Lauritsen*® find that the K-U end 
point of the N™ spectrum agrees with the dis- 
integration data, whereas a lower end point 
does not. 

If we attempt to choose between the Fermi 
and the K-U theory by considering only the 
upper region of the spectrum, where the scatter- 
ing contemplated by Konopinski®> is probably 
very small, we find that in the case of Mn** the 
two theories can be brought into practical agree- 
ment over the range from 1.2 to 2.1 Mev. This 
is clearly brought out in Fig. 10 which shows on 
enlarged scale the data used in constructing Fig. 
3, together with the best Fermi distribution to 
fit the data between the limits mentioned. 

In view of this, it might perhaps best be said 
that the cloud’ chamber is incapable of dis- 
tinguishing between the two theories, at least in 
this case. 

Confidence in cloud-chamber measurements is 
further shaken by consideration of the disagree- 
ment between the chamber and the 8-ray spec- 
trograph reported by DuBridge and Marshall.*° 
Their work shows that the cloud chamber can 
make some electrons appear to have more energy 
than they really possess. 

The agreement between the 1.7-Mev y-ray 
line and the difference in energy of two groups 
in the Mn** spectrum is quite striking, and 
would seem to support the K-U analysis of this 
spectrum. 


— Delsasso and Lauritsen, Phys. Rev. 49, 561 
(1936). 
% 1. A. DuBridge and J. Marshall, Phys. Rev. 56, 629 
(1939). 
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At the very high temperatures and densities which must 
exist in the interior of contracting stars during the later 
stages of their evolution, one must expect a special type 
of nuclear processes accompanied by the emission of a large 
number of neutrinos. These neutrinos penetrating almost 
without difficulty the body of the star, must carry away 
very large amounts of energy and prevent the central 
temperature from rising above a certain limit. This must 
cause a rapid contraction of the stellar body ultimately 
resulting in a catastrophic collapse. It is shown that energy 
losses through the neutrinos produced in reactions between 


free electrons and oxygen nuclei can cause a complete 
collapse of the star within the time period of half an hour. 
Although the main energy losses in such collapses are due 
to neutrino emission which escapes direct observation. 
the heating of the body of a collapsing star must necessarily 
lead to the rapid expansion of the outer layers and the 
tremendous increase of luminosity. It is suggested that 
stellar collapses of this kind are responsible for the phe- 
nomena of movae and supernovae, the difference between 
the two being probably due to the difference of their 
masses. 


$1. INTRODUCTION 


NE of the most peculiar phenomena which 
we encounter in the evolutionary life of 
stars consists in vast stellar explosions known as 
“ordinary novae”’ and “supernovae.” It is now 
well established that, although these two classes 
of novae possess a great many features in com- 
mon, they are sharply separated insofar as their 
maximum luminosities are concerned. The ordi- 
nary novae, appearing at the rate of about 50 per 
vear in our stellar system, reach at their maxi- 
mum a luminosity of the order of magnitude of 
10° suns. On the other hand, the supernovae, 
flaring up in any given stellar system only once in 
several centuries, reach luminosities exceeding 
that of the sun by a factor of 10%. The inter- 
mediate luminosities have never been observed, 
and there seems to exist a real gap between these 
two classes of stellar explosions. 
The common features of novae and supernovae 
may be briefly summarized as follows: 
(1) Both ordinary novae and supernovae show 
a very similar form of luminosity curve (apart 
from the luminosity scale, of course) with a sharp 
rise to the maximum within a few days or weeks, 
and a subsequent slow decline of intensity, 
decreasing by a factor of two every four or five 
months. 
(2) In both cases the spectrum shows rather 
high surface temperature (up to 20,000°C for 


* Now at George Washington University as Fellow of 
the Guggenheim Foundation. 


ordinary novae, and probably above 30,000°C 
for supernovae), and the rapid expansion of the 
stellar atmosphere which is evidently blown up 
by the increasing radiative pressure. In the case 
of Nova Aquilae 1918, for example, the star was 
surrounded by a luminous gas shell expanding 
with a velocity of 2000 kilometers per second, 
whereas the gas masses expelled by the galactic 
supernova of the year A.D. 1054 (observed by 
Chinese astronomers) form at present an ex- 
tensive luminous cloud known as the “Crab- 
Nebula.”’ It must be noticed here that the large 
surface area of this blown-up atmosphere is 
mainly responsible for the observed high lumi- 
nosities, since the increase of the surface tempera- 
ture can only account for a factor of several 
hundreds in the surface brightness. 

(3) Whereas the “prenovae,”’ in the rare cases 
when they have been observed, represent com- 
paratively normal stars of the spectral class A 
(surface temperature about 10,000°C),' the 
“postnovae,”’ remaining after the flare-up, possess 
extremely high surface temperatures (spectral 
class O) and seem to represent highly collapsed 
configurations, such as the stars of the ‘‘Wolf- 
Rayet”’ type. 

The same evidently holds true for the case of 


1The meagerness of observational material makes it 
impossible to decide whether “‘prenovae” are located on 
the main sequence or to the left of it. There seems to be 
no doubt, however, that, as the result of ee. the 
position of the star in the Hertzsprung-Russell diagram is 
strongly shifted towards higher surface temperatures and 
smaller radii. 
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supernovae, since the star found in the center of 
the “‘Crab-Nebula,” and representing most proba- 
bly the remainder of the galactic supernova A.D. 
1054, shows the typical features of a very dense 
“white dwarf.” 

The change of state caused by these stellar 
catastrophes strongly suggests that the process 
involved here is not connected with any instan- 
taneous liberation of intranuclear energy due to 
some explosive reaction, but rather represents a 
rapid collapse of the entire stellar body as was 
first suggested by Milne.’ 

It was suggested by Baade and Zwicky,’ in 
particular application to supernovae, that such 
collapses may be due to the formation within the 
star of a large number of neutrons which would 
permit considerably closer ‘‘packing” in the 
central regions. 

It is easy to see, however, that the possibility 
of closer packing alone is quite insufficient to 
explain the rapid collapse of the star, since such a 
collapse requires the removal of large amounts of 
gravitational energy produced by contraction in 
the interior of the star. In fact, quite independent 
of whether the matter in the center consists of 
charged nuclei or neutrons, the heat produced by 
contraction must pass through the entire body of 
the star, and the rate of energy transport, de- 
pending on the opacity of the main body, will be 
the same in both cases. On the other hand, if we 
can find some way of removing instantaneously 
the heat liberated in the central regions, in spite 
of the opacity of the stellar body, the star will 
collapse with a velocity comparable to that of 
“free fall’’ independent of the kind of particles 
existing in its interior. 

The amount of gravitational energy which is 
liberated when a star of a mass M contracts from 
the original radius Rp to the collapsed radius R, is 
given approximately by 


GM? 


avecur(—-— = 


= : 1 
R. Ro R. 


Thus, for example, if a star of the mass and 
radius of our sun contracts to the size of the 


2E. A. Milne, Observatory 54, 145 (1931). 


3W. Baade and F. Zwicky, Proc. Nat. Acad. Sci. 20, 
259 (1934). 


companion of Sirius (R.= Ro/40),‘ the total liber- 
ation of gravitational energy will be of the order 
of magnitude of 10°° erg. On the other hand, the 
time of the free-fall collapse from the original 
radius Ro to any small value of the radius is given 
approximately by 


(2) 


In the case of the sun, At will be of the order 
10° sec. (i.e., about half an hour), so that the mean 
rate of energy removal necessary for such a 
collapse is about 10" erg/g sec. If we remember 
that the collapse of novae and supernovae takes 
place within a few days, we come to the con- 
clusion that the rate of energy removal in these 
actual cases may be only several hundred times 
smaller than given above. 

As we suggested in a recent publication,® this 
very fast removal of energy from the interior of 
the star can be understood on the basis of the 
present. ideas on the role of neutrinos in nuclear 
transformations involving emission or absorption 
of 8-particles. In fact, when the temperature and 
density in the interior of a contracting star reach 
certain values depending on the kind of nuclei 
involved, we should expect processes of the type 


zN4+e-—z_,N4+antineutrino (3) 


which we shall call, for brevity, ‘‘urca-processes.” 
The neutrinos formed in the above processes*® 
absorb a considerable part of the transformation 
energy (about 3), and escape with practically no 
difficulty through the body of the star. 

As we shall see later, these processes of ab- 
sorption and reemission of free electrons by 
certain atomic nuclei which are abundant in 
stellar matter may lead to such tremendous 
energy losses through the neutrino emission that 


‘The companion of Sirius possesses a mass which is 
almost equal to the mass of the sun, and may be considered 
as representing the type of the white dwarfs obtained by 
the collapse of the sun. 

5G. Gamow and M. Schoenberg, Phys. Rev. 58, 1117 


* We shall use the term “neutrino” both for ordinary 
neutrinos and antineutrinos involved in the reaction, 
since there is no noticeable difference in their behavior. 
It is also clear that one can neglect the possibility of mutual 
annthilation of these particles within a stellar body, since 
they escape from the star with practically no collisions. 
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the collapse of the entire stellar body with an 
almost free-fall velocity becomes quite possible. 

Before discussing in more detail the charac- 
teristic features of stellar collapse caused by the 
urca-processes taking place in the hot interior, we 
must develop the formulas for the dependence of 
“neutrino losses”’ on the temperature and density 
of matter. 


§2. ENERGY LossEs THROUGH NEUTRINO 
EMISSION 


Since within a contracting star the central 
temperature gradually increases from compara- 
tively low to very high values, different urca- 
processes will assume importance during the 
different stages of the contraction and collapse. 
In the very beginning, these processes will 
involve only those nuclei which can capture free 
electrons with the energy of few kilovolts, whereas 
at the height of the collapse the nuclei requiring 
the energy balance of several Mev may become 
of primary importance. Thus, it becomes neces- 
sary to have the formulas both for nonrelativistic 
and relativistic electrons. In the relativistic 
calculations we shall limit ourselves to the ex- 
treme relativistic case, since the calculations in 
the intermediate region are necessarily much 
more complicated. The values for the inter- 
mediate temperatures can be obtained with suffi- 
cient approximation by a simple interpolation. 

We shall, moreover, consider only the case of 
an ideal electron gas, which, as we shall see later, 
is of primary importance for the stars having 
large masses. It is clear that in the case of 
degeneracy the decreasing number of free energy 
levels will reduce considerably the rate of urca- 
processes, and that these processes will stop 
entirely when the electron gas becomes completely 
degenerated. 


(A) Nonrelativistic case 


Let us consider a unit mass of stellar matter 
containing n, free electrons and mz° nuclei of 
atomic number Z, which can capture these 
electrons and go over into unstable nuclei of 
atomic number Z-—-1. Since the matter in the 
stellar interior is almost completely ionized, we 
may write for the number of free electrons 


3p/mn, (4) 


where p is the density of matter and my the mass 
of the hydrogen atom.’ For the number of nuclei 
participating in the urca-process we have 


evidently 
nz°=czp/Az-mu, (S) 


where cz and Az are the concentration and 
atomic weight of the isotope in question. 

If the electron gas is not degenerated, the 
number of electrons with energy between E and 
(E+dE) is given by the Maxwellian expression 


n(E)dE EXdE. (6) 


The average life of an electron before being 
captured by the nucleus through the emission of 
antineutrinos may be calculated on the basis of 
Fermi’s theory of 6-decay, and is given by*® 


7(E) =m Ig 2-h'c*/g*nz(E—Q)?, (7) 


where g is the Fermi constant and Q the maxi- 
mum electron energy of continuous #-spec- 
trum 

The number of electrons captured by the nuclei 
in question per unit volume per unit time is 


n(E) 
N-= ——dE. 
J, 7(£) 


Using (6) and (7) we get 
2g°nzn. 
lg 2c*h4 


or 
N-= kT)? (10) 
a! lg 2cth! kT 


where the integral 


Q 
e-*(x —xo)xidx (11) 


is written in respect to a new variable 
x=(E/kT)(x0=Q/kT). (12) 


Estimating this integral by the saddle-point 


(k —Q)*EWE (9) 
Q 


? This expression presupposes the absence of hydrogen, 
since the contraction of a star can start only after all the 
hydrogen is consumed by (C—H) and (H —H) reactions. 

8’ This expression differs from the formula given by 
H. A. Bethe and R. F. Bacher [Rev. Mod. Phys. 8, 83 
(1936) ], since the latter was obtained on the basis of the 
now abandoned Konopinski-Uhlenbeck form of Fermi's 
theory. 
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method, we get 
1(Q/RT) =1(x0) = (13) 


On the other hand, the number of electrons 
emitted by the unstable nuclei (Z—1), in the 
energy interval E, (E+dE), is, according to 
Fermi 

(14) 


dNt= 


and the total number of emitted electrons per 
second per unit volume 


N+=0.152————— 
v2 r*hic* 


(15) 


In the case of equilibrium, when an equal 
number of electrons is absorbed and emitted, we 
must have N+= N~- which leads to 


= 2 16 
nz kT Vix (16) 


or, since mz-1+nz=nz° 


(17) 


2 


and 


x 


The energy taken up by antineutrinos ejected 
in the process of electron capture per unit time 
per unit volume can be easily calculated and 
found to be 


mig 2htc®? \RT 2 


while the energy of neutrinos accompanying 
electron emission is 


W® = (20) 


where \ is the decay constant of the unstable 
nucleus (Z—1). Comparing (19) and (20), and 
bearing in mind the expression for \, we get 


(21) 


so that for the total energy of neutrino radiation 
per unit time per unit volume we have 


W=W®+Wo 
(22) 


Czp 


(18) 


A 


(B) Extreme relativistic case 


In this case the energy distribution between 
the thermal electrons is given by the expression 


(23) 


which is generally correct also for the inter- 
mediate velocities. Here p, is the momentum, and 
E includes the rest energy mc? of the electron. 
The function f(R7) has a simple form only in the 
extreme relativistic case, and is then given by 


S(RT) =c-(kT)-* 
x [1 +(1 (24) 


Formula (7) for the capture of free electrons 
contains only the energy (E—Q) of the ejected 
neutrinos, which were treated relativistically be- 
cause of their negligible mass, and may also be 
used directly in the present calculations. Thus, 
combining (7) with (23) and (24) we get 


ART) -(RT)I, (25) 


where Q’=Q+mce? (i.e., the total emission energy 
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of the electron), and tivistic expression 
mc2 gm*c 2 
0 kT. 2n*h\ 
This integral can be easily calculated and is with e=E/mc*, e=Q’/mce*. Integrating we have 
given by 
, , t= (29) 
(27) 
For the number of electrons emitted in the The number of stable and unstable nuclei in 


energy interval E, (E+dE) we have the rela- equilibrium will now be given by 


= 30 
yp" [12(RT)*+ 60’°RT JL 1+ (1+  Azmy (30) 


and 
10% myp[12(kT)? +60 RT +07 + (1 +moc?/kT)? X . (31) 


zmn 


The total energy taken up by antineutrinos ejected in the process of electron capture is 


g’ngn, moc?\ 277 
m Ig 


whereas the neutrino-energy in emission 


W® =3Q0dnz-1, (33) 
where A is again the decay constant. Expressing W™ through W®,° we get 
kT 


Q 


and, after some simplifications, we obtain for the total energy of neutrino radiation per unit volume 
per unit time 


(35) 


(C) Nonrelativistic gas and relativistic $-electrons 
Of some special interest is the case when the temperature of electron gas is not yet high enough to 
treat the main bulk of particles relativistically, whereas, on the other hand, the electrons responsible 
for urca-process are those from the relativistic Maxwell distribution tail. 
It can be easily calculated that in such case the number of nuclei in Z and Z —1 states will be given 
by 
0.98 X czp 
140.98 my 
1 Czp 


140.98 my 


(36) 


and 
(37) 


® By means of (26), (32), and the equilibrium relation N*= N-. 
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For W® we have as before 


whereas W“ is given by 


M. SCHOENBERG 


W® =nz_\-d-4Q, (38) 


Wo=7.1 dal 


(39) 


Q (Q/kT+5/2)! 
Thus for the total energy of neutrino emission per gram per second we have in this case approximately 


(40) 


§3. UrcA-PROCESSES WITH VARIOUS ELEMENTS 


Since practically any stable nucleus leads to an 
unstable isobar after capturing one electron, 
there must be, altogether, many hundreds of 
possible urca-processes presenting a large display 
of all possible energies and decay constants. The 
dependence between the decay energy and the 
decay constant permits arranging all such proc- 
esses in one row, beginning with very slow 
processes, at comparatively low temperatures, 
and ending with the very fast ones requiring 
extremely high thermal velocities." 

Since our present list of the known Fermi 
elements is very limited and contains only a small 
fraction of all possible B-active nuclei, the choice 
of a particular reaction, or group of reactions, 
which would be of special importance for the 
theory of stellar collapse, presents serious diffi- 
culties. This choice can be, however, somewhat 
simplified if we remember that the importance of 
any given reaction depends not only on its 
absolute efficiency but also, to a large extent, on 
the abundance of the element in question in 
stellar material. 

As an example of urca-process which would 
take place at a comparatively low temperature, 
we may mention the case of the light helium 
isotope. The process will take place according to 
the scheme 


He?+e-—H'+antineutrino 
H*—He’+e-+ neutrino. 


The decay of H* was recently investigated by 
O’Neal and Goldhaber"™ who found 15+3 kv for 
the maximum energy of 8-particles and 7 X10~'° 
sec.—! for the decay constant (decay period being 


(41) 


” The effect of permitted and nonpermitted §-transitions 
impairs somewhat the above regularity, and must be 
we) ty: separately in each particular case. 
asm. . O'Neal and M. Goldhaber, Phys. Rev. 58, 574 


about 30 years). Since the energy of 15 kv 
corresponds to the temperature of 1.2 108°C, 
i.e., only six times higher than that in the center 
of the sun, we see that the reaction must take 
place with some speed even under solar conditions. 

However, because of the small value of \, the 
neutrino losses due to this reaction are com- 
paratively unimportant. In fact, even when most 
of the nuclei He* are transformed into H®*, the 
energy of neutrinos emitted per gram per second 
will be given by 


2 Cz erg 
W® X Xez,_ (42) 


3 g sec. 


a comparable amount being carried away by the 
antineutrinos emitted in the capture process. 
Thus, even for a star constructed entirely of 
He’ (cz=1), the energy losses are several million 
times lower than are necessary for a rapid 
collapse. There is, however, every reason to 
believe that He* is not very abundant in the 
stellar interior because of its transformation into 
ordinary helium through the capture of a thermal 
proton. The only other known £-decaying 
nucleus with low energy balance is that of the 
still less abundant radioactive element RaD 
(Emax235 kv, 421.4X10-* sec.-'), but there 
must surely be more Fermi elements of this class 
which are as yet undiscovered because of their 
comparatively weak activity. 

As an example of an intermediate case we can 
give the urca-process with the main iron isotope 
which reacts according to the scheme 


e564 e-—+Mn*6+ antineutrino 


Mn**—>Fe**+ neutrino. (43) 


This iron isotope forms probably about 12 
percent of stellar matter and the nucleus Mn°®* is 
known to possess a maximum electron energy of 


| | 
1 
| 
s 
| I 
I; 
| / 
| | 


THEORY OF STELLAR COLLAPSE 


Fic. 1. Energy-losses through the neutrino emission in 
a contracting star with M=5M(©) and w=1.7 (H com- 
pletely substituted by He). Curve A for the urca-process 
in Fe** (E=1.7 Mev; \=7.7X10- sec.-!; concentration 
10 percent). Curve B for the urca-process in He® or similar 
element (E = 15 kv; \=7 X 107 sec. !; assumed concentra- 
tion one percent). Broken line represents the average 
radiative energy losses per unit mass calculated from the 
ordinary luminosity formula for the contracting model. 
“Normal radius” represents the radius which the star had 
while still in the main sequence. 


1.7 10° ev and the decay constant of 7.710-° 
sec.! (decay period 2.5 hours).” 

This reaction, which becomes saturated at the 
temperature of 10!°°C, will remove in the form of 
neutrino radiation about 10" erg per second per 
gram of stellar matter. 

As the final example we may take ordinary 
oxygen representing about 30 percent of stellar 
material. In this case the urca-process takes place 
according to the scheme 


(44) 


The energy balance and the decay constant of 
N'® are ev and 7.7X10~ sec.—! (decay 
period 9 sec.), respectively, which leads to the 
energy losses of about 10'* erg per second per 
gram at the saturation point. 

Comparing this value with the energy losses 
necessary for the free-fall collapse, we find that 
such a collapse can take place even if only one 
percent of the stellar body participates in the 
urca-process. 

In order to get some idea of the energy losses in 
different stages of stellar contraction, we calcu- 
late numerically the rate of the urca-process with 
iron in a contracting star of mass equal to ten 


See, for example, G. Gamow, Structure of Atomic 
Nuclei (Oxford University Press, 1937), Table F. 


545 


sun masses. The temperature and density distri- 
bution within a gravitationally contracting star 
is very closely represented by the Emden’s 
polytrope with the index 3, and the values in the 
center are given by the formulas" 


T.=1.2X10-M/R, 
pc=13M/R', 


(45) 
(46) 


where Mand R are the mass and the radius of the 
star. Taking M=5-M(©)=1X10™ g, and ex- 
pressing R in terms of solar radius R(©)=7 X 10"° 
cm we get 


°C, 
385-[R/R(O) 


Assuming the data given above for the iron 
reaction, and using formulas (18), (22), (31) and 
(35), we get for the energy losses through neutrino 
emission (per second per gram) the values 
represented in Fig. 1. Since the luminosity of a 
contracting star, as the function of its mass and 
radius, is given by the expression" 


[R/R(O) erg/sec. 


we see that, even if the urca-process is limited to 
one percent of the stellar mass, neutrino losses 
become comparable to radiative losses at 
R=0.1R(©). For R=0.01R(©) the losses 
through neutrino emission will reach the tre- 
mendous value of 10" erg/g sec. 


(47) 
(48) 


(49) 


§4. DyNAMICS OF THE COLLAPSE 


Since the hydrodynamical equations describing 
the collapse of a stellar body, cooled from the 
center by neutrino emission, are necessarily very 
complicated, we shall limit ourselves in the 
present article to the entirely qualitative 


discussion. 
It must be clear first of all that in different 
stages of stellar contraction the most important 


13 See, for example, A. Eddington, Internal Constitution 
of Stars (Cambridge University Press, 1926). We take here 
the mean molecular weight of stellar matter to be equal to 
1.7, which corresponds to 35 percent helium, and 65 percent 
of heavier elements. 

“4 A. Eddington, reference 16. In this expression, as well 
as formulas (45) to (48), we assume » = 1.7, and the opacity 
equal to that of the Russell mixture, forming 65 percent 
of stellar matter (the rest being helium). 
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role is played by urca-processes with different 
elements. In the very beginning, when the central 
temperature is comparatively low, neutrino 
cooling can be produced only by the elements 
with the low value of energy balance. As we have 
seen above, energy losses due to such elements 
are rather small, even at the saturation point, so 
that one may expect here only a slight acceleration 
of normal gravitational contraction. However, as 
the central temperature rises, new and more 
powerful urca-processes are introduced, the rate 
of contraction increases, and there results a 
catastrophic collapse. 

It is important to keep in mind that, since 
urca-processes take place in the regions of highest 
temperature, the energy flow towards these 
points, either through radiative or through 
convective mechanisms, is entirely excluded. 

Thus, in order to maintain a central pressure 
necessary for the support of the outer layers, the 
star must develop in its interior a rapidly growing 
central condensation. The heat which will be 
generated in compressing the extra material into 
this small central region will be rapidly removed 
by neutrino emission, and the whole process will 
be somewhat analogous to the condensation of 
water vapor contained in a vertical vessel with 
the bottom cooled by liquid air. 

However, whereas in the above example the 
vertical vapor column, unsupported from be- 
neath, will simply fall down as a whole, gas 
masses forming the body of a collapsing star will 
be strongly heated by compression. 

Since this heat cannot escape towards the 
center (except from the layers in the immediate 
vicinity of the urca-region), part of it must 
remain in the collapsing body and part be radi- 
ated from the surface, thus increasing the stellar 
luminosity. It is not difficult to see that, under 
such conditions, inner parts of a stellar body will 
continue to move towards the center, whereas the 
outer regions will begin to expand. 

Such an expansion of the outer parts of the 
star, resulting in the decrease of their opacity and 
the increase of radiating surface, must necessarily 
lead to a rapid rise of luminosity which we ob- 
serve as the ‘‘flare-up”’ in the nova phenomena. It 
must be noticed that this “flare-up” probably 
takes place as soon as the collapse of the interior 
has progressed to a high degree. At a certain 
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stage the rapidly increasing radiative pressure 
overbalances the weight of stellar atmosphere, 
and large masses of gas begin to flow outside 
forming the luminous shell characterizing all 
nova and supernova explosions. We must stress 
here the point that, according to the proposed 
theory, the increase of luminosity during stellar 
collapse is an entirely secondary phenomenon, so 
that in the first approximation these changes of 
luminosity may be easily neglected. 

Let us turn now to the question of the ultimate 
fate of a collapsing star, and the possible reasons 
for the difference between the ordinary novae and 
supernovae. 

It was shown by Chandrasekhar" that, whereas 
the stars with mass smaller than 5.75./(©)/u? 
(=2-M(©) for the molecular weight »=1.7) 
possess stable configurations with a certain mini- 
mum radius (white dwarfs), the stars having the 
largest masses are subject to unlimited contrac- 
tion.'® In view of this fact, it is natural to expect 
that, for the star masses smaller than the above- 
given critical value, the collapse will occur on a 
much smaller scale than for heavier bodies. The 
increasing condensation in the interior of such 
light collapsing stars will sooner or later lead to 
electron degeneracy which will slow down and 
finally stop all urca-processes. For these small 
scale collapses we should expect that the final 
state will be represented by a white dwarf of a 
mass comparable to that of the original star, and 
that the amount of material thrown out in the 
form of a gas shell will be comparatively small. 
This expectation seems to be in good agreement 
with observations of ordinary nova phenomena, 
where the expelled gas shell takes up only about 
one-hundredth of a percent of stellar mass, and is 
completely dissolved in space several years after 
the explosion. 

On the other hand, the stars possessing a mass 
larger than the critical one will undergo a much 
more extensive collapse, and their ever-increasing 
radiation will drive away more and more material 
from their surface. The process will probably not 


% See, for example, S. Chandrasekhar, Introduction into 
the Study of Stellar Structure (Chicago University Press, 
1939), Chapter IV. 

16 This results from the fact that in very heavy stars 
the velocities of electrons in the generated Fermi gas 
approach the velocity of light, after which such a gas is 
unable to support the weight of the stellar body. y 
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stop until the expelled material brings the mass of 
the remaining star below the critical value. This 
process may be compared with the supernovae 
explosions, in which case the expelled gases form 
extensive nebulosities (as ‘‘Crab-Nebula”’ or 
“Filamentary Nebula’), apparently assuming 
this permanent state after the explosion. This 
point of view seems also to acquire some confir- 
mation from the fact that the relative number of 
novae and supernovae occurring in stellar sys- 
tems is of the same order of magnitude as the 
relative number of stars having masses smaller 
and larger than 2- M(©). 


CONCLUSION 


In the present article we have developed the 
general views regarding the role of neutrino 
emission in the vast stellar catastrophes known to 
astronomy. It must be emphasized that, while 
the neutrinos are still considered as highly 
hypothetical particles because of the failure of all 
efforts made to detect them, the phenomena of 
which we are making use in our considerations 
are supported by the direct experimental evi- 
dence of nuclear physics. In fact, the experiments 
of Ellis and Wooster’ and of Meitner and 


Orthman" leave no doubt that the energy balance 
does not hold in the processes of radioactive B-trans- 
formations, and all later evidence on this subject 
strongly indicates that this ‘‘disappearance of 
energy” occurs in such a way that it appears to be 
carried away by particles of almost unlimited 
penetrability. 

Whereas the fundamental ideas of the proposed 
theory are very simple and the physical part of 
calculations pertaining to the rate of neutrino 
emission can be easily carried out on the basis of 
existing formalism, the problem of the dynamics 
of the collapse represents very serious mathe- 
matical difficulties. It is to be hoped that these 
difficulties will be overcome by the choice of some 
suitable simplified model. 

It is our pleasant duty to express our thanks to 
Drs. S. Chandrasekhar, E. Fermi, E. Teller and 
M. A. Tuve for helpful discussion and interest in 
the problem. One of us (M. Schoenberg) is 
grateful to the John Simon Guggenheim Memo- 
rial Foundation for the grant which enabled him 
to carry out this research. 


17 C. D. Ellis and W. A. Wooster, Proc. Roy. Soc. A117, 
109 (1927). 
(1930) Meitner and W. Orthman, Zeits. f. Physik 60, 143 
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Measurements of the meson lifetime were made by comparing the lead absorption curves 
of the cosmic radiation at two different elevations, and by compensating the air column at the 
higher level by a layer of graphite placed on top of the counter telescope. The method does 
not require any knowledge of the energy distribution or of the height of the producing layer 
of the mesons. All systematic errors are greatly reduced by the differential method used and 
can be shown to be insignificant. The value of the lifetime comes out to be appreciably shorter 
than in previous determinations, viz. (1.25+0.3) X 10~-*(uc?/108 ev) sec. 


1. THE METHOD 


FTER the instability of the mesons was 
suggested by Yukawa! from theoretical 


1H. Yukawa, Proc. Phys. Math. Soc. Japan 20, 319 
(1938). 


reasons, experimental evidence of various kinds 
for this effect was discussed by several authors, 
notably by Euler and Heisenberg.? A more 


*H. Euler and W. Heisenberg, Ergeb. d. exakt. Natur- 
wiss. 17, 1 (1938). 
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recent survey of all the earlier observations has 
been given by Rossi.* While there seemed to be 
good support for the decay hypothesis with 
values of the lifetime of two to four micro- 
seconds, it could by no means be said that it was 
definitely confirmed. Furthermore, all previous 
methods necessitated assumptions regarding the 
energy distribution of the mesons or the height 
of their production or both, in order to obtain a 
numerical value for the lifetime. This holds true 


YA; 
Y 


Fic. 1. Experimental arrangement. 


also for all the more recent determinations of the 
lifetime, which have been published since.*~? 
We have tried, therefore, to set up an experiment 
which would be independent of such additional 
data and which would be free from objections 
which might cast doubt on the reality of the 
effect. 

The only feasible method to select a definite 
energy range in counter experiments is by 
differential absorption. The mesons seem to 
follow quite closely a range absorption law‘ in 
good agreement with Bloch’s formula for energy 
loss by ionization. Every point on an absorption 
curve corresponds, therefore, to a definite energy. 

The procedure adopted was the following one. 
The absorption curves in lead were taken at two 
different elevations. In place of the absent air 


3B. Rossi, Rev. Mod. Phys. 11, 296 (1939). 

4M. A. Pomerantz, Phys. Rev. 57, 3 (1940). 

5 B. Rossi, Norman Hilberry and J. Barton Hoag, Phys. 
Rev. 57, 461 (1940). 

®M. Ageno, G. Bernardini, N. B. Cacciapuoti, B. 
Ferretti and G. C. Wick, Ricerca Scient. 10, 1073 (1939); 
Phys. Rev. 57, 945 (1940). 
( + V. Neher and H. G. Stever, Phys. Rev. 58, 766 

1940). 


column at the higher level a layer of graphite of 
equal mass was put on top of the counter tele- 
scope. As the absorbing power of carbon and air 
is practically the same, according to Bloch’s 
theory, the only difference in the two series of 
observations consists in the compression of the 
air mass into the smaller volume of the graphite 
where the mesons have less chance for sponta- 
neous decay. If N(z) is the number of counts per 
unit time at a lead thickness z, then dN/dz will 
be proportional to the number of particles corre- 
sponding to the energy of particles of range z in 
lead. The attenuation factor A, measuring the 
fraction of particles surviving at the lower level 
I is therefore given by 


A(E) = (ON,/dz 0N11/02) z=2(E)» 


where the indices J and JI refer to lower and 
upper level, respectively. 

Recently Fermi® has pointed out that the 
density of packing of the absorbing material has 
a direct influence on its stopping power. The 
effect is due to partial screening of the contribu- 
tions of atoms at some distance from the 
passing particle. Formulas for this effect have 
been developed by Fermi and by Halpern and 
Hall.® A correction for this effect can be easily 
applied in our case, and it turns out to be 
quite small. 

A preliminary run with the described arrange- 
ment” gave, simultaneously with the somewhat 
similar measurements by Rossi and _ collabo- 
rators,> the first unambiguous proof for the 
decay effect. In the present series we have tried 


TaBLe I. Counting rate for different lead absorber thicknesses, 


DuRHAM Mt. 


cT. HR. CT./HR. cT. HR. CT./HR. 


0} 8291 198 41.9403 | 4179 66.5 62.8+0.6 
3| 7601 235.25 32.340.3 | 2989 70.25 42.5+0.5 
6; 7980 271.5 29.4+0.2 | 4226 113 37.4+0.4 
9 
2 
5 


2412) =81.33 29.6+0.4 
5523 196.25 28.14%0.2 | 4627 135.75 34.1+0.3 
3793 140.75 26.940.3 | 3606 115.25 31.3+0.3 
3913 150.75 25.95+0.3 | 4129 141.5 29.2+0.3 
2563 99.75 25.7+0.3 


8 E. Fermi, Phys. Rev. 56, 1242 (1939); 57, 485 (1940). 

®O. Halpern and H. Hall, Phys. Rev. 57, 459 (1940). 

10 W. M. Nielsen, C. M. Ryerson, L. W. Nordheim and 
-K. Z. Morgan, Phys. Rev. 57, 158 (1940). 
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to hold all corrections for side effects to a mini- 
mum and we hope to have obtained a more 
reliable determination of the lifetime than here- 
tofore available. 


2. EXPERIMENTAL PROCEDURE AND DATA 


The absorption curves of the cosmic radiation 
of vertical incidence were measured for lead from 
0 to 15 inches, in Durham, North Carolina 
(elevation 125 meters, mean barometric pressure 
75 cm Hg) and on Mount Mitchell, North 
Carolina (elevation 2040 meters, mean baro- 
metric pressure 60 cm Hg). In the latter observa- 
tions, the vertical air column of 15 cm Hg=203 
g/cm* between the elevations of Durham and 
Mount Mitchell was compensated for by a 
layer of graphite of a thickness of 200 g/cm?. As 
shown in Fig. 1, the graphite filled completely 
the sensitive cone of the counter telescope. 
Observations were made at the higher elevation 
in the Park Warden’s cabin on Mount Mitchell. 

An additional layer of wood was placed on the 
roof of the small shack on the Duke University 
campus where the observations at the lower 
elevation were performed. In this way the 
amount of wood above the apparatus in both 
cases was of the same thickness of 14 inches. 

Fourfold coincidences were recorded between 
counters arranged as shown in Fig. 1 with 
various thicknesses of absorbing lead between 
them. The counters, the shielding lead at the 
sides, and the absorbing lead layers between, 
were all supported on an angle iron frame not 
indicated in the figure. The same frame was, of 
course, used in both series of observations, and 
in both cases the axis of the counters made 
approximately the same angle with respect to 
the east-west direction. The three bottom 
counters were shielded by four inches of lead on 
each side. The upper counter was shielded by 
two inches of lead on each side as indicated in 
the figure. The pieces of shielding lead and the 
absorbing layers between counters were approxi- 
mately three times the active length of the 
counters. 

The counters were approximately 5 cm in 
diameter and had an effective length of 20 
centimeters. They were prepared and filled with 
hydrogen in a manner somewhat similar to that 


suggested by Shonka." Tests showed that the 
counters had an efficiency of 95 percent or more. 
The counter and recording circuits will be pub- 
lished elsewhere. In a series of measurements of 
this kind it is quite essential that the over-all 
efficiency of the counting apparatus be constant 
over both series of measurements. In addition to 
the frequent checks on counting rates of the 
individual counters we have assured ourselves of 
the reliability of the apparatus by repeating the 


CTS PER WHOUR 


ZT IN INCHES PB 


Fic. 2. Counting rate as a function of lead absorber 
thickness for Durham (curve I) and Mount Mitchell 
(curve IT). 


low altitude observation for two absorbing thick- 
nesses of 6” and 15” after the observations on 
Mount Mitchell were completed. These repeat 
measurements are the lower entries of the 
Durham data for the corresponding thicknesses 
as shown in Table I. The two series of observa- 
tions made in Durham which bracket the Mount 
Mitchell observations agree well within the 
statistical errors of the observation. The measure- 
ments here reported are also in very satisfactory 
agreement with the results reported earlier.'® 
The complete 1940 data-for both series of 
observations are shown in Table I and the 
absorption curves are plotted in Fig. 2. 


3. EVALUATION OF THE LIFETIME 
The energy lost by ionization in an air column 
of 1900 m is of the same order of magnitude as 


1 See J. B. Hoag, Electron and Nuclear Physics (D. Van 
Nostrand, New York, 1938), p. 432. 
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in about 10 inches of lead. In the evaluation of 
the lifetime in terms of the observed attenuation 
factor, the variation of energy along the path 
has to be taken into account. The probability 
of decay w in an interval dx is 


wdx = (dx/cr)(uc*/pe), (2) 


where r is the lifetime at rest, u the rest mass and 
p the momentum of the meson. We assume for 
the loss of momentum a linear law, i.e., the 
momentum given at a point x backwards from 
the point of observation will be 


pc = poc +x, (3) 


where the specific loss in air 8 can be considered 
as constant for kinetic energies of the order of 
the rest energy and over. The attenuation factor 
becomes then 


A(poc) =exp (-f ) 


1 poctBzo 1 yc? 
B poc CT pe 


poc ) uc*/erB 
= 4) 
poc+Bxo 


In this formula fo is the final momentum of 
the particle and Bx» the amount of energy lost 
by ionization in the air column. In (4) the 
variation with altitude of the density of the air 
has been neglected. We have not introduced this 
refinement, as it only would give a correction of 
the order of one percent in A. 

Solving (4) for r we obtain 


=— lg lg —. 5) 
cB poc A 


ue? poc + Bxo / 1 

The specific energy loss in air for particles with 
a kinetic energy slightly larger than their rest 
energy is B=1.9X10° ev/g, i.e., 2.3108 ev/km 
air. With x9=1.92 km and referring 7 to a 
standard mass of 10* ev and measuring it in 


microseconds, we obtain 


uc? 
)x10-* sec., (6) 
108 ev 
+4.4 1 
y= 1.45 (7) 
pe A 


where pc is expressed in units 10° ev. 
The energy loss in lead of a meson is = 1.2 X 10° 
ev/g, i.e., 0.35 X 108 ev/inch. For kinetic energies 


TABLE II. Attenuation factor, average energy and lifetime 
for various intervals on lead absorption curve. 


INTERVAL 


INCHES PB 1/A Po py 
Oto 3 2.12 
3to 6 1.76 2.45 2.62 
6to 9 2.64 3.50 1.22 
9 to 12 2.35 4.55 1.14 

12 to 15 2.13 5.6 1.12 


greater than 108 ev the pc scale as a function of 
lead thickness is 


pc = (0.9+0.35z) X 108 ev. (8) 


The end correction has been obtained from the 
range-energy relation as given by Bloch’s for- 
mula for particles with mass ~ 200 m. 

Table II gives the values of 1/A, i.e., the ratio 
of slopes of curves I and II for various intervals 
of lead thicknesses together with the average 
value of pc and the resulting y (lifetime in 
microseconds). 

The values of y agree very well for the intervals 
between 6 and 15 inches. The interval between 3 
and 6, however, seems to give a much larger 
value. This behavior finds its natural explanation 
in the role of the soft component for not too 
large thicknesses of the lead absorber. An in- 
spection of the absorption curves of Fig. 2 
shows a sharp upturn in the interval in question. 
This suggests that the soft component makes a 
noticeable contribution, of the order of 3 to 4 
percent, to the counting rate at 3 inches, as 
indicated by the dotted extrapolation lines in 
Fig. 2. An effect of this magnitude is indeed to 
be expected from the theory of the soft com- 
ponent. It can be seen from an estimate by 
Heitler” that a fraction of order of 75 of the 


2 W. Heitler, Nature 140, 235 (1937). 
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soft radiation should be able to produce co- 
incidences with 80 g/cm? lead interposed, and 
this would give approximately 1/30 of the total 
counting rate. At 6 inches, however, the con- 
tribution of the soft component should be quite 
negligible. The correction for the soft component 
might even be larger for curve I (Durham) 
since electrons coming from the meson decay in 
the air column would give an additional con- 
tribution. If we subtract from the counting rate 
at 3 inches on both curves the same amount of 
1.4 counts per hour (1/30 of 42) we would bring 
the value of y down to 1.6, and the sketched 
extrapolation of the curves would bring it to 1.35. 
The effect of the soft component is thus certainly 
of the order of magnitude required to remove 
the above discrepancy. The exact value of the 
correction can, of course, not be determined, but 
we can safely discard the point at 3 inches for 
the evaluation of the lifetime. 

For the intervals between 6 and 15 inches 
the value of y is practically constant, while the 
attenuation factor itself varies appreciably (com- 
pare Table II). This behavior means a confirma- 
tion of the Lorentz time factor, i.e., the energy 
dependence of the mean free path against decay. 

In order to obtain an estimate of the probable 
error introduced by the errors in the counting 
rates we can proceed as follows. In view of the 
smoothness of the curves we can consider the 
whole interval between 6 and 15 inches. The 
average value of pc in this interval is 4.4108 
ev.8 With the probable errors listed in Table I 
we obtain for this interval 1/A = 2.38+0.32. The 
corresponding values of y are 


2.70 
1.01. 


2.38 
1.16 


1/A 2.06 
1.38 


We feel safe, therefore, to give as a final result, 
apart from systematic errors to be discussed in 
the next section 


vy =1.2+0.3. 


‘8 This corresponds to 10 inches Pb. The most consistent 
way to determine the average momentum for a finite 
interval is to take the geometrical mean 


551 


5. SYSTEMATIC ERRORS. COMPARISON WITH 
OTHER OBSERVERS 


Our results are not very sensitive against 
changes in the assumed values of stopping power 
of air and lead, and corrections for such changes 
could easily be applied. It might, however, seem 
somewhat more doubtful to identify the energy 
scale on the lead absorber with the values from 
the ionization losses alone. The correct procedure 
would be, of course, to take the geometrical 
average of the energy of the distribution taken 
out by all occurring processes in a unit layer of 
lead. Any additional effect besides ionization, 
i.e., scattering, nuclear absorption and such, 
would increase the meson energy and therefore 
tend to lower the value of the lifetime. In view 
of the behavior of the geometrical mean this 
would, however, introduce a serious deviation 
only if the number of particles of high energy 
absorbed would be at least comparable with 
those lost by ionization. 

In our treatment we assume that no mesons 
are generated between Mount Mitchell and 
Durham. We also assume that there is no 
appreciable admixture of protons. 

There are three differences in the operation of 
the counter telescope at the two locations at 
different altitude, which might give rise to 
systematic errors and which have to be dis- 
cussed. 

(I) The mesons, before entering the telescope, 
emerge from different materials.—It will be much 
more likely that a meson is accompanied by soft 
secondaries when coming from a dense material. 
It might happen therefore, that the meson misses 
the first counter, which, however, responds to a 
secondary. As the probability for such an event 
is independent of the lead absorber between the 
counters this effect would give a constant per- 
centage correction to the curve II (at the higher 
elevation). It would reduce, therefore, the value 
of 1/A by the same fraction. A one percent 
change in 1/A would produce only 1.5 percent 
change in the lifetime. This effect is therefore 
wholly negligible, even if we should allow one 
percent for it as Rossi does. 

(Il) The soft component is of different intensity 
outside the cone of the telescope; Auger showers.— 
Owing to the heavy shielding of the counter 
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telescope only showers of very high energy could 
possibly give a contribution to the counting rate. 
The rate of increase of large showers and bursts 
from sea level to an altitude of 2000 m is about" 
4:1. Therefore, in case the lead filters between 
the counters have a noticeable effect on the 
chance of recording such showers a correction 
for these counts must be made. 

In our preliminary measurements,'® observa- 
tions were made without any side shielding and 
with side shielding of the lowest counter. The 
reduction of the counting rate (amounting to 
one count per hour) by the shielding was prac- 
tically independent of the lead thickness be- 
tween the counters, and we concluded, therefore, 
that the large air showers did not affect appre- 
ciably our lifetime determination. To make sure 
of this for our present experiment, we made 
measurements in Durham with one of the 
counters placed out of line as shown in Fig. 1. 
The counting rate then was very low, 0.13+0.03 
per hour with 6 inches Pb between the counters, 
and 0.09+0.02 per hour with 15 inches Pb. It 
seems doubtful whether the slight dependence 
on lead thickness is at all significant. However, if 
we take these figures as a true measure of the 
contributions of large sidewise showers in Dur- 
ham and admit an increase by a factor 4 for the 
rate on Mount Mitchell, the correction to the 
lifetime would only be one percent, and if we 
take the limits 0.16 count per hour at 6 inches 
Pb and 0.07 at 15 inches Pb, the correction would 
still be only 3 percent. We feel sure, therefore, 
that the effect of large showers can be entirely 
neglected compared to the statistical errors. 

(III) The stopping power of the air and carbon 
columns of equal weight is not quite the same.— 
The differences due to the different nuclei in 
ionization losses and probably also nuclear ab- 
sorption if any, should be negligible. As already 
mentioned, however, the packing density itself 
has some influence on the stopping power, as 
discussed by Fermi.® It is particularly important 
to consider this effect carefully as it cannot be 
distinguished readily from a true decay. 

Formulas for the change in stopping power 
due to the density of the traversed material have 
been given by Fermi® and by Halpern and Hall.® 


4 Compare for instance, Norman Hilberry, Chicago 
meeting, Phys. Rev. 59, 112A (1941). 


In sufficient approximation, the specific energy 
loss is diminished compared to the values given 
by the Bloch formula by 


>1. (9) 


dE e—1 e—1 
A—= lg 
dx mv? L 1—v?/c? 1—v?/c? 
Here N is the number of electrons per unit 
volume and ¢ an effective dielectric constant. 
According to Halpern and Hall ¢ is given by 


Ne? 
e—1= 


(10) 
ev =3.5Z X10-™ sec.-! 


where vy is the geometrical average over the 
atomic frequencies. This formula gives a good 
approximation if 


Ne?/xrmv? > 1 —v?/c? 


for the largest frequencies occurring. Otherwise 
the effect is reduced considerably. 

A correction of this type disturbs, of course, 
the balancing of the air column for curve I by 
the carbon for curve II. We can restore this 
balance by shifting the curves against each other, 
i.e., by comparing the slopes not at the same 
lead thickness, but at points where the total 
energy loss of the mesons just stopped is the 
same. As our curves show only a very small 
curvature in the range used (6 to 15 inches Pb) 
the error committed in neglecting the Fermi 
effect can only be rather small. 

Expression (10) gives for carbon «= 1.13. The 
total deficiency in stopping power of 200 g of 
carbon for particles with 1/(1—v?/c?)!=E/ye? 
=6.6 is then calculated from (9) to be 0.108 
10° ev, which corresponds to 0.31 inch of Pb. 
This is, therefore, the amount the two curves 
have to be shifted against each other. The above 
value for v corresponds to the average energy in 
the carbon absorber of those particles which are 
stopped by 10 inches of lead. The differences in 
counting rate for the three intervals used are 3.3, 
2.8, 2.18. The second differences are thus 0.5 and 
0.67, or an average of about 0.6 per three 
inches on a three-inch interval, i.e., 0.6 per inch 
shift for a nine-inch interval. The counting rate 
difference of 8.2 between 6 and 15 inches has, 
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therefore, to be corrected by 0.31 X0.6=0.19 to 
8.04, giving an 1/A of 2.32 in place of 2.38 as 
formerly. This corresponds to an increase in the 
lifetime y by 3 percent only. The comparatively 
small value of the dielectric constant ¢« seems to 
agree fairly well with the measurements of Crane 
and collaborators® which showed only a rather 
small deviation from the Bloch formula in 
carbon for electrons with E/mc?=20. Even if 
we would take «=2, which value is very im- 
probable, the correction to be applied to our 
value of y would be only 11 percent. 

All the systematic effects discussed here give, 
thus, corrections much smaller than the probable 
limits of error due to statistics. We believe, 
therefore, 


7 =(1.25+0.3) (uc?/108 ev) X 10~ sec. 


to be a fair representation of our results. 

Our value for the lifetime is thus considerably 
shorter, about one-half, than heretofore assumed. 
A direct comparison with previous observations 
can only be made with those of Rossi et al. We 
obtain a mean free path 


6x 37.4 


L= = 1.92 km ——=(8.0+0.8) km 
2.94 


lg 1/A 


for mesons which can traverse 15.2 cm Pb, as 
compared to Rossi’s value of (9.4+1.6) km for 
particles which can traverse 12.7 cm Pb. Our 
value is somewhat smaller, the more so in view of 
the thicker lead absorber, but well within the 
limits of error. This means that the average 
energy of the mesons must be definitely larger 
than the value calculated by Rossi from Blackett's 
energy distribution. 


*H. R. Crane, N. L. Oleson and K. T. Chao, Phys. 
Rev. 57, 664 (1940). 
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Similar remarks apply to the failure of other 
observers to find a lifetime as short as one 
microsecond by inclination measurements and 
other methods. All these measurements imply a 
knowledge of either the height where the mesons 
are produced, or their average energy or both. 
It is generally assumed that most of the mesons 
are produced near the maximum of the soft 
component, but apart from a_ hypothetical 
genetic connection, there are no facts in support 
of this assumption. The observations of Schein 
and collaborators'® seem to indicate, on the 
contrary, that the production of mesons extends 
to considerably lower altitudes. The new value 
of the lifetime will, however, be of importance 
for the discussion of the temperature effect and 
of the fraction of the soft component due to 
meson decay and similar questions. 

Two of the authors (W. M. N. and L. W. N.) 
wish to acknowledge grants from the Research 
Council of Duke University which have made 
possible these and earlier observations on meson 
lifetime. We wish also to acknowledge the kind- 
ness of Mr. T. W. Moore, Superintendent of 
State Parks, in making available the facilities on 
top of Mount Mitchell. We are also indebted to 
Mr. Robert Wilson and Mr. Charles Thompson, 
weather observers on Mount Mitchell, and par- 
ticularly to Mr. Ed. Wilson, Park Warden, for 
their very helpful cooperation. We are also in- 
debted to the National Carbon Company for 
making available the large graphite blocks. One 
of us (K. Z. M.) wishes to take this opportunity 
to acknowledge a grant from the North Carolina 
Academy of Science for aid in a previous in- 
vestigation."” 


16M. Schein, E. O. Wollan and Gerhart Groetzinger, 
Phys. Rev. 58, 1027 (1940). 
( x8): M. Nielsen and K. Z. Morgan, Phys. Rev. 54, 245 
1938). 
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Theoretical values for the intensity of the soft component due to meson decay are compared 


with the experimental data. 


HE lifetime of the meson has been de- 
termined by the Duke University group! 
to be 


ev) sec; y=1.25+0.3. 


We believe this value to be the most reliable one 
obtained so far. The considerably smaller value 
compared to those reported previously? necessi- 
tates a reconsideration of the role of the decay 
part of the soft component. It will be seen that, 
with the data available at present, a rather 
critical situation is reached which might possibly 
lead to important consequences as to the nature 
of the meson decay. 

The intensity of the soft component produced 
by the meson decay depends essentially on the 
lifetime of the meson and the fraction a of the 
total mesonic energy given to electrons and 
quanta. The counting rate in one of the usual 
arrangements will depend sensitively on the 
energy distribution of the electron component in 
equilibrium with the mesons. The corresponding 
problem of the theory of showers, i.e., the 
determination of the electron track length in 
showers in their dependence on energy has been 
solved in reasonable approximation by Tamm 
and Belenky.* Their result can be expressed as 
follows: 

A single electron of energy Ey produces a total 
track length with energy larger than E 


Ey 
E). (1) 


1 See the preceding paper by W. M. Nielsen et al., Phys. 
Rev. 59, 547 (1941). 

2 The only other measurements comparable to ours, by 
B. Rossi, Norman Hilberry and J. Barton Hoag, Phys. 
Rev. 57, 461 (1940), also led to a rather short lifetime. 
Our value for the ratio 7/y is still only one-half of theirs. 
assh Tamm and S. Belenky, J. Phys. U.S.S.R. 1, 177 


Here £ is the specific energy loss and 


* 2.4E 
e(Eo, E)= ——dx; e= » (2) 
x 


€ 7 


where £; is the energy at which ionization and 
radiation losses for electrons are equal. For 
E.>E the upper limit of the integral (2) can be 
taken as infinite and ¢ can then be expressed by an 
exponential integral function. For air ev) 
the function g(«, E) has the following values: 


Eev | | 4X10* | 10X108 25x 108 {100% 108 
e | 1 | 08 | 067 | OS | 0.25 


The limiting value g(0)=1 expresses the fact 
that all energy will finally be dissipated by 
ionization, whatever multiplication processes 
may have taken place before. 

The probability of decay of a meson of energy 
E, per unit path is* 


If we assume that every meson gives the fraction 
a of its energy to electrons or photons of compa- 
rable energy and that f(£,)dE, is the distribution 
function of the mesons, then the total energy 
given to the electrons per unit length is 


E,f(E,) uc? | 
AEv=a f dE,=a—N, (3) 
ct cr 


where N is the total number of mesons. Inserting 
this for Ey into (1) and dividing by N, we obtain 
the ratio of electrons with energy larger than E to 


‘The number of mesons of kinetic energy ~wuc* and 
less is rather small. In the following we neglect corrections 
due to the difference between - and E and also those due 
to the continuous energy distribution of the decay electrons 
produced by mesons of definite energy. 
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the total number of mesons 


Net uc* 
r= =a— 9(E) (4) 
mes T 


and with our value for yc?/cr 
r=1.3—9(E). (5) 
Y 


The value for E to be inserted into (5) depends 
on the characteristics of the counter telescope. 
For the usual arrangements a lower limit between 
4 and 10X10 ev might be assumed which would 
give g¢~0.7 to 0.8. 

The same ratio (5) with g=1 (i.e., all electronic 
energy used up in ionization) will also hold for the 
total ionization® produced by the decay electrons 
and the hard component. 

The empirical situation regarding the ratio of 
soft to hard component seems still to be rather 
ambiguous. Counter telescopes operated at sea 
level in open air give a ratio ~4 : 3 for counts 
without any absorber and with 10 to 15 cm Pb.° 
This means a ratio 720.33 for the total soft 
component and the most absorbable part of the 
mesons to the mesons which can penetrate the 
lead absorber. The chief uncertainty of this 
method consists in the possible neglect of coher- 
ence between rays of the soft component which 
might lead to an underestimate of the soft 
intensity. Measurements with thin-walled ioniza- 
tion chambers without and with lead shielding 
give a considerably higher ratio’ r~0.5 to 0.6. 
None of these measurements have, however, 
been made under ideal conditions for the present 
purpose. The chief uncertainties in such measure- 
ments lie in the effects of the walls and of the 
radioactivity of the surroundings which may 
depend on the shielding and probably tend to 
increase r. The values for r are, of course, for the 


5B. Rossi, Phys. Rev. 57, 469 (1940) has used this 
argument already with slight refinements which, however, 
give rise only to insignificant corrections. The values used 
by him correspond to y =2.5. 

® Compare for instance the Durham curve of the pre- 
ceding paper which has been taken under the best possible 
conditions. See also G. Bernardini et al., Phys. Rev. 58, 
1017 (1940). 

7H. Schindler, Zeits. f. Physik 72, 625 (1931); H. 
Hoerlin, Zeits. f. Physik 102, 652 (1936); J. C. Street and 
R. T. Young, Phys. Rev. 46, 832 (1934); 52, 552 (1937). 


total of the soft component including knock-on 
and primary electrons. Attempts to resolve this 
total ratio into the various contributions seem 
still to be rather uncertain. An estimate by 
Euler® seems to show that the effect of primary 
electrons is at sea level only a few percent of the 
total. The percentage of knock-on electrons is 
about 10 percent of the hard component ac- 
cording to the best estimate given by Tamm and 
Belenky.* 

The principal hypotheses regarding the meson 
decay and their consequences are as follows: 

I. The meson has the spin $ and decays into an 
electron and a photon. In this case a=1. With 
= 1.25, r would be ~0.8 for counters and ~ 1 for 
the total ionization. These values seem to be 
quite incompatible with the empirical data, and 
this hypothesis can be definitely ruled out. 

II. The meson has the spin 1 and decays into 
an electron and a neutrino. In this case a=} and 
r~0.4 for counters and ~0.5 for total ionization. 
These figures give about twice the amount of 
r~0.2 to 0.25 which follows from the counter 
experiments, and which has been assumed for the 
decay component by previous investigators.® ® 
The only way to remove this discrepancy would 
be to assume a rather high degree of coherence 
between the rays of the soft component and 
possibly of the soft and the hard rays.°® 

On the face of the present evidence it does not 
seem assured that even a=} can be admitted. 
Smaller values of a would mean, of course, that 
more complex modes of disintegration would 
have to be assumed with simultaneous creation 
of a number of neutrinos. One possibility would 
be the disintegration of a meson with spin } into 
an electron and two neutrinos (a= 4). Or else the 
meson disintegration would have to be considered 
as a kind of explosion and not as a simple 8-decay 
process. In this case, of course, no information 
regarding the spin of the meson could be ob- 
tained. It seems highly important in any case 
that this situation should be clarified by further 
experiments. 

®H. Euler, Zeits. f. Physik 116, 73 (1940). 

*An indication for such a coherence effect might be 
found in the fact that a comparison of counter and ioniza- 
tion-chamber measurements tends to give higher values 
for the specific ionization of cosmic rays than direct 


countings of droplets; compare the summary by Thomas 
H. Johnson, Rev. Mod. Phys. 10, 209 (1938). 
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To determine the exchange properties of the neutron- 
proton interaction, the existing knowledge of the neutron- 
proton force in states of even parity, gained from the 
properties of the deuteron ground state and low energy 
scattering experiments, must be supplemented by informa- 
tion concerning the interaction in odd parity states, 
information which can be obtained only by observations 
on high energy neutron-proton scattering and deuteron 
photo-disintegration by energetic y-rays. Calculations 
have been performed for three types of interactions with 
the purpose of testing the sensitivity of such experiments 
to variations of the exchange operator dependence of the 
interaction. These interactions are analogous, in isotopic 
spin dependence, to the potentials predicted by three 
forms of current mesotron theory: (I) “Symmetrical” 
(II) “Charged” (IIT) ‘“‘Neutral.”” With the interaction in 
even parity states described by rectangular well potentials 
with constants adjusted to fit the binding energy and 
quadripole moment of the deuteron, and the cross section 
for slow neutron-proton scattering, each of these three 
potentials makes a definite prediction concerning the 
interaction in odd parity states. The scattering calculations 
were performed for a neutron energy of 15.3 Mev. The 
results for total cross sections and angular distributions 


in the center of mass system are as follows: 
(I) «=0.621 cm?, 

~(1 —0.080 cos +0.077 cos* 3) 
(II) o=0.666 x 10-** cm?, 

o(8)~(1 +0.126 cos +0.042 cos? 3) 
(III) o«=0.983 x 

o(3)~(1+0.932 cos +0.457 cos? 8). 


The energy of the Li+H y-rays (hw=17.5 Mev) was 
adopted for the computations on photo-disintegration. 
The three theories under discussion predict the following 
electric dipole total cross sections and angular distributions 


(1) =0.768 X 10°77 cm’, (8) ~ (sin? 8 +0.015) 
(II) o=0.723X10-" cm?, o«(8%)~(sin? 8 +0.07/) 
¢=0.376XK10-" cm’, ~(sin? 8 +0.36). 


Calculations have also been performed for the small cross 
sections arising from magnetic dipole and electric quadri- 
pole absorption. The spherically symmetrical term in the 
electric dipole angular distribution is a consequence of 
the noncentral forces invoked to explain the deuteron 
quadripole moment. High energy photo-disintegration 
angular measurements thus constitute the most sensitive 
test of both the isotopic spin dependence of the neutron- 
proton interaction and the existence of noncentral forces. 


A DETAILED empirical determination of 
the neutron-proton interaction is funda- 
mental to the development of a mesotron field 
theory adequate to explain nuclear phenomena. 
A previous paper! has been devoted to the 
discussion of the neutron-proton interaction in 
states of even parity, based on properties of the 
ground state of the deuteron and low energy 
collision phenomena. Such information should 
serve to specify the spin coupling between the 
mesotron field and nuclear particles. To investi- 
gate the charge properties of the mesotron field, 
_ it is necessary to study the exchange nature of 
the neutron-proton interaction, or equivalently, 
the interaction in states of odd parity. The 
requisite information can be obtained only by 
consideration of high energy continuum states 


* On sabbatical leave from Brooklyn College, Brooklyn, 
New York. 

1 W. Rarita and J. Schwinger, Phys. Rev. 59, 436 (1941). 
This paper is referred to as NPI in the text. 


of the neutron-proton system, realized in experi- 
ments on neutron-proton scattering and the 
photo-disintegration of the deuteron. It is the 
program of this paper? to discuss the experi- 
mental implications of several force laws, and 
guided by their characteristic predictions, to 
determine the most critical experimental test of 
an assumed interaction. 


EXCHANGE INTERACTIONS 


The calculations of this paper will be per- 
formed for three types of exchange forces with 
the purpose of testing the sensitivity of experi- 
ments on high energy neutron-proton scattering 
and deuteron photo-disintegration by energetic 
y-rays to the exchange properties of the neutron- 
proton interaction. The adopted interaction 


?Some of the results of this work were published in a 
short paper written in collaboration with H. A. Nye: W. 
Gai . Schwinger and H. A. Nye, Phys. Rev. 59, 209 
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energies are: 

(1) t2{Jo(r) Siz}, 

(II) Jo'(r) + Ji (r)o1- 
(1) 

V=— (r) + Si" S12}, 


Si. a2 


expressed in terms of isotopic spin operators 
“1, t2, the spin operators o;, o2 and the relative 
position vector r. The particular isotopic spin 
dependence of the various potentials is motivated 
by current theories, the potentials (I), (II) and 
(III) being respectively analogous to the force 
laws predicted by the symmetrical, charged and 
neutral mesotron theories. It should be stressed, 
however, that any evidence offered on behalf of 
one of these potentials should not necessarily be 
construed as support for the corresponding 
mechanism suggested by the present, incomplete 
mesotron theories. Our previous paper (NPI) 
constituted a discussion of the interaction in 
states of even parity based on the simplifying 
choice of rectangular well potentials for the 
various coupling terms, of equal range but 
different depths, i.e., 


V (even) = — (2) 


with g and y constants and J(r) a rectangular 
well function of depth Vo and range ro. This 
specification of the interaction in even states, 
combined with the special isotopic spin de- 
pendence of potentials (I), (II) and (IID), 
provides a complete determination of these 
several interactions. The elementary isotopic 
spin dependence of these interactions permits a 
simple representation of the interaction in odd 
states in terms of that in even states, viz. : 


(I) V (oad) = —} coven): 
V (oad) =—3! V woven): 
(II) (odd) —*V ceven); 
(3) 
1V (oad) =-! V even); 
(IIT) 3 V waa) =3 


V coda) =! V (even): 


The various odd interactions thus depend 
strongly on the total spin in view of their 
relation to the coupling energies in the singlet 
and triplet even states: 


1V (even) = — {1—2g}J(r), 
(even) = {1+75S12} J(r). 


The constants describing the potential wells 
have been evaluated in NPI and are rewritten 
for reference : 


g=0.0715, y=0.775, V»o=13.89 Mev, (5) 
ro = 2.80 X 10-" cm. 


(4) 


SCATTERING OF FAsT NEUTRONS BY PROTONS 


The discussion of NPI was confined to the 
treatment of the neutron-proton scattering 
associated with the *S,+D, continuum state, a 
valid approximation at low energies. We wish 
to consider neutrons of sufficient energy that the 
and *P,+*F, states are also perturbed. 
The usual scattering theory must be extended 
to describe this process, for the interactions, and 
therefore the radial functions, differ in the 
various states of total angular momentum J. 
To develop such a theory it is convenient to 
have simple operational representations of the 
*P,, and wave functions. The spin and 
angular dependence of the *Py and *P; wave 
functions may be expressed symbolically by 


That these forms do indeed correspond to triplet 
state functions of the correct total angular 
momentum and magnetic quantum number m is 
immediately evident from their spin symmetry 
and rotation properties. The P state character 
is verified by noticing that they are spherical 
harmonics of the first order. To perform the 
scattering analysis, it is necessary that the 
P component of an incident plane wave 
exp [ik-r viz: 
i(kr) k-r 


34 
kr kr 


(k=(ME/h*)', g:(p) ~ —cos p) 
be decomposed into its *Po, *P:, *P2 constituents. 


and 


(6) 


_| 
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The *P, and *P; components are: 


The latter assumes the simpler form 
gil r 
~im— 
2 kr2 r 


when the total spin is quantized in the direction of wave propagation. The remainder of the P wave 
function (6) is then the *P, component of the plane wave. The P wave function describing the 
system under the influence of a coupling is obtained by replacing the radial function of free motion, 
gi, by the separate radial functions vo, v; and v2 in the *Po, *P; and *P:2 constituents of the P wave 
function. The complete P function then obtained is: 


e'yo(r) 1 r 3. (r) 1 
kr 2 r 2 y @ r 
1 r 3 1 r 
kr kr 2 r 2 2 r 


The phase factors e,1.2 are inserted to guarantee that the wave function satisfies the boundary 
condition of a scattering problem, in virtue of the asymptotic forms of the P radial functions: 
Vo,1,.2~ —cos (kr+70,1,2). It should be remarked that this treatment of the *P2 state neglects its 
interaction with *F». 

The spherically diverging wave describing the scattered particles becomes asymptotically : 


3e'": sin sin (no — n2)—(@1 — 2) 
kr kr 2 r 


3 1 r 
sin (m— m2) (ert) (9) 
r 


The cross section for scattering through an angle # into unit solid angle is obtained by computing 
the squared absolute value of the coefficient multiplying e*"/r and averaging the result over the 
three magnetic sublevels associated with the quantum number m. This evaluation is simplified 
by noting that 

x0" = (K/k)x1", (10) 
if k is the quantization direction. The amplitude of the scattered wave (9) is then of the form of a 
spin operator A acting on x:”. The absolute squares of the amplitudes summed over the triplet 
states is simply the diagonal sum of A+A minus the absolute square of the diagonal matrix element 
of A in the singlet state. The angular scattering cross section, thus evaluated, is: 


= (1/k*)[ (sin? no+3 sin? sin? n2) cos? 3 
— (3 sin? (no— m2) +2 sin® (m1—2))(3 cos? 3—1)]. (11) 


The total cross section for the P scattering is composed additively of the cross sections for the several 
J states, multiplied by appropriate statistical weights. 


3 5 
sin? yo-+— sin® sin? (12) 
R?L9 9 9 
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To include the interference between *P and *S, scattering, we must add to (9) the spherically 
diverging *S; wave’ 


eikr 1 
(13) 


r 2tk 
where 


is the complex *S; phase shift associated with the magnetic substate m. By a spin summation calcu- 


- lation similar to that above, we obtain the interference contribution to the angular cross section: 


1 1 3 5 1 ( 
cos sin? sin? nts sin® (sin no sin (29 — no) +sin ne sin (2x — 
(sin m sin —m)-+sin m2 sin (15) 


The interference between *D, and ‘P is also easily evaluated, but the rather complicated expressions 
will not be written since these interference terms are not numerically significant. A formula equivalent 
to (11) has been derived by Kittel and Breit.‘ Further our interference formula (15) reduces to theirs 
when the *S,; phase shifts are considered real and independent of m (ordinary potential). 

To determine the effective potentials in the various *P states, the spin-spin operator S;. must 
be evaluated. The states *Py and *P; are unmixed and must therefore be eigenfunctions of Sj». 
It is immediately verified from the explicit forms of these functions that the eigenvalues are *P»: 
Si2= —4, *P: Sip=2. Indeed these are the only eigenvalues of S12 since (Si:2+1)?=9 in triplet states. 
Since we shall neglect the coupling of *P, with *F; the effective potential in the *P, state is obtained 
by inserting the diagonal value of S2. This value is most easily obtained by noting that the diagonal 
matrix element of S12 in a *P,y state multiplied by its statistical weight, 2/+1, and summed over 
the three possible J values must vanish since this sum equals the diagonal sum of Sj. Therefore 
the required number is —}{[1X —4+3X2]=—#2. Thus, within the range of interaction, the po- 
tentials *V(.aa) of the ‘‘neutral”’ theory (III) are *Po: 29.16 Mev, *P;: —35.41 Mev, *P2.: —9.58 Mev. 
The potentials of the ‘“‘charged”’ theory (II) have the reversed sign, while those of the “symmetrical” 
theory (I) are obtained from (II) through division by 3. 

Explicit calculations have been performed for a neutron energy of 15.3 Mev. The *P phases 
no, 1,2 were evaluated for the three potentials with the results indicated below 


No m 
(1) 0.0745 — 0.0545 —0.0175 
(11) 0.5310 — 0.1146 — 0.0462 (16) 
—0.1025 0.9953 0.0731. 


It may be noted that the signs of these phases are in accordance with the attractive or repulsive 
nature of the corresponding potentials. The P angular cross sections computed from Eq. (11), 
with the phases appropriate to the three potentials are: 


(1) = (1/k*) [0.0038 +0.0045 cos? 3], 
(II) =(1/k*) [0.1028 — 0.0020 cos? 8], (17) 
(III) o(8) = (1/k*)[0.4871+0.6871 cos? 3] 
with 
4x /k? =0.6823 X 10-** cm’. (18) 


3 The notation conforms to that of NP 


I. 
‘C. Kittel and G. Breit, Phys. Rev. 56, 744 (1939). 
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The four complex phases describing the *S,+*D, scattering have been evaluated by the procedure 
outlined in NPI, with the following numerical results: 


ko = — 1.6618, Ko = — 1.6735, x2 = — 0.0035, = — 0.0152, 


(19) 
= —0.0811, = =0.0532. 
The *S,+*D, angular cross section computed from Eqs. 31, 35 and 37 of NPI is then: 
= (1/k?)[0.9764+0.0464 cos? ]. (20) 


The complete triplet differential angular cross sections, including *P, *S,+*D,, and *S,, *P interference 
scattering, become : 


(I) o(8)d2 =0.6804[0.9830+0.0017 cos +0.0510 cos? #](d2/4m) X 10-24 cm?, 
(II) o()d2 =0.7464[0.9865 +0.1931 cos 8 +0.0406 cos? #](d2/4m) X 10-24 (21) 
o(8)dQ = 1.1654[0.8569+0.8484 cos +0.4294 cos? #](d2/4r) X 10-24 cm?. 


The angular factors in brackets are so normalized that their average over all directions is unity; 
the initial multiplicative numbers thus represent the total cross section in 10-*4 cm? units. 

These triplet cross sections must be appropriately averaged with singlet cross sections to give 
numbers suitable for comparison with experiment. Within the range of interaction, the 'P interaction 
energy of the “‘neutral’’ theory (III) equals that in the 'S state: —11.90 Mev. The interaction 
energy of ‘‘charged” theory (II) has the opposite sign, while that of the ‘“‘symmetrical’’ theory (I) 
is three times larger than (II). The 'S+'P scattering cross sections were computed by the usual 


scattering theory, with the following results: 


(I) o(8)d2 = 0.4445[0.9392 — 0.4380 cos +0.1824 cos? ](dQ/4r) X 10-74 
(II) a(3)d2=0.4236[0.9855 — 0.2402 cos +0.0437 cos? 3 ](dQ/4mr) X10-*4 cm?, (22) 
(IIT) +0.4981 cos d+0.1338 cos? 3 ](dQ/4r) X 10-*4 


The total cross sections obtained by combining 
the singlet and triplet cross sections in the ratio 
of statistical weights are for the three theories: 


(1) o =0.621 X10-%4 cm?, 
(II) =0.666 X 10-24 cm?2, (23) 
(111) = 0.983 X 10-24 cm?, 


while the angular distributions are of the form: 
o(8)~(1—0.080 cos 8 +0.077 cos? 8), 
(II) o(8)~(1+0.126 cos 3+0.042 cos? (24) 
(III) o(8)~(1+0.932 cos 3+0.457 cos? 8). 

The only measurements which have been 
performed at this energy are those of Salant, 
Roberts and Wang, and Salant and Ramsey,‘ 
who employed the energetic neutrons emitted in 
the Li+D reaction. The scattering cross section 


5E. O. Salant, R. B. Roberts and P. Wang, Phys. Rev. 
55, 984 (1939). E. O. Salant and N. Ramsey, Phys. Rev. 
57, 1075 (1940). 


obtained by these authors is o=(0.66+0.07) 
X10-* cm’, referring to a neutron energy of 15 
Mev. This result is consistent with both theory I 
and theory II, but would appear to definitely 
exclude the potential IIT. 

Recently, Powell, Heitler and Champion® 
have measured the angular distribution of the 
recoil protons produced by collisions with the 
8.7-Mev neutrons emitted in the B+D reaction. 
Their data, although consistent with spherically 
symmetrical scattering in the center of gravity 
system, do not preclude small deviations from 
isotropy. In particular, the angular distribution 
predicted by the “neutral’’ theory, the most 
asymmetric of the three: 


~(1+0.229 cos 3+0.159 cos? 8) (25) 


cannot be considered in definite disagreement 
with the observations. 


°C. F. Powell, H. Heitler and F. C. Champion, Nature 
146, 716 (1940). 
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THE PHOTO-DISINTEGRATION OF THE DEUTERON 


This section will be devoted to a discussion of 
deuteron disintegration produced by electric 
dipole, magnetic dipole and electric quadripole 
y-ray absorption. The cross section for the 
absorption of a y-ray of energy hw incident in 


e hw 
hc Mc? 3 m, m’ 


The required average over the initial magnetic 
substates and summation over the final substates 
is indicated in the formula. The three absorption 
processes under discussion are thus described by 
the matrix elements of er/2, the electric dipole 
vector, (eh/2Mc)M, the magnetic dipole vector, 
and err/4, the quadripole moment dyadic of 
the system. For completeness we shall list all 
the various final states associated with the three 
processes : 


Electric dipole : 3P,+ 


Electric quadripole: *S,+'D,, *D3+%Gs, 
Magnetic dipole : “De, *Si+*D1, *De. 


It should be noted that both triplet and singlet 
states may serve as final states for the magnetic 
dipole absorption process. The total photo-cross 
section is composed additively of the separate 
cross sections for these processes. In treating 
the angular distribution of the disintegration 
products, however, it is necessary to include 
interference between states of the same total 
spin; triplet and singlet states do not interfere, 
of course. At the high energies contemplated in 
this paper, the photoelectric disintegration is 
the dominant process and we shall include only 
the interference of the *P states. We shall also 
calculate the cross sections for the electric 
quadripole and magnetic dipole disintegration, 
chiefly for their theoretical interest. 

The initial state function ¥; is the ground 
state wave function of the deuteron (NPI 


Eq. (2)): 


The final state of the electric dipole transition is 
described by the Eq. (J8), but with the modifi- 
cation that the phase factors e*%,1,2 are replaced 


1 h 
2 2Mc 


the direction of the unit vector x with unit 
polarization vector e, inducing a transition of 
the deuteron system from an initial state WV; to 
any final state W, describing emission of the 
disintegration products within the solid angle 
dQ about the propagation vector k, is expressed by 


(26) 


by e~'%,:.2 to satisfy the asymptotic boundary 
condition of a plane wave exp [ik-r] and a 
converging spherical wave,’ appropriate to this 
problem. In evaluating the matrix element of 
the electric dipole moment, the operator Sj. 
must be replaced by the numbers —4, 2, —? 
when acting on the *Po, *P; and *P: wave 
functions, respectively. The further technique 
of calculation is identical with that employed in 


the treatment of scattering. The resultant - 


differential cross section, restricted to unpolar- 
ized y-rays, is then 
e? 


o(3)dQ=—- — — ) 
3 hc hk 


1 

——(3 cos? 
72 

cos (no— n2) 


dQ 
cos (ni—2))}—, (28) 
4dr 


expressed in terms of 3, the angle of emission 
with respect to the direction of the y-ray, and 
the three radial integrals 


J roo(u—2hw)dr, 


rv\(u+2-‘w)dr, (29) 


J ho Jar 


7 Cf. N. F. Mott and H. S. W. Massey, The Theory of 
Atomic Collisions (Oxford University Press, New York, 
1933), p. 258. The eo of this reference applies to 
ordinary entials. The required generalization is given 
in Appendix IT. 
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The total photoelectric disintegration cross 
section is simply : 
re? Mw il 


— — (30) 
hk 9 


It should be remarked that there is no inter- 
ference between *P,) and *P, states, for in 
virtue of their property of being eigenfunctions 
of Siz associated with different eigenvalues, the 
spin scalar product vanishes. The formulae (28) 
and (30) reduce correctly to the Eqs. (75) and 
(76) of NPI when the particles in P states are 
regarded as unperturbed vo, 1,2= 1. 

The photomagnetic transition to the singlet 
state including the interference of 'So and 'D, 
has already been discussed in NPI (Eqs. (79) 
and (80)), but subject to the approximation of 
disregarding the effect of forces on the 'Dz wave 
function. This approximation is easily removed ; 
the modifications required are the replacement of 
g2 by the exact wave function wo( ~ — sin (kr +42)) 
in the radial integrals, and the substitution of 
cos (59—52) for the cos 49 occurring in the 
interference term of NPI Eq. (79). 

The treatment of both electric quadripole and 
magnetic dipole transitions to the triplet state 
need not be more than approximate, for these 
effects are quite small. In the electric quadripole 
transition we shall disregard the effect of spin- 
spin forces on the D states, describing them by 
ordinary potentials. We shall further neglect 
the transitions between the continuum D states 
and *D, portion of the ground state. The total 
cross section for this process is then: 


) “| (f 
240 hc Mc? \|Ey|7 k 0 


+( J |, (31) 


wherein the u; and w; are, respectively, the *S 
and *D continuum radial wave functions. The 
magnetic dipole transition to the triplet state is 
of some interest, for its very existence is a 
consequence of the spin-spin forces, insofar as 
these forces affect the ground state of the 
deuteron. The cross section, computed with an 
approximation regarding the continuum D states 


similar to that above, is simply: 


The numerical evaluations of these cross 
sections and angular distributions have been 
performed at the energy of the y-rays emitted 
in the Li+H reaction: Aw=17.5 Mev. The total 
cross sections for photoelectric dipole absorption 
are computed to be: 


(I) cm? 
El. dip. (II) o=0.723X10-? cm? (33) 
(III) o=0.376X 


To explain the small cross section obtained from 
the ‘‘neutral” theory (III) it is necessary to 
observe that a strong attractive force in the 
final state diminishes the electric dipole matrix 
element, when the wave-length of the dissociated 
particles is comparable with the size of the 
deuteron. This effect, which is particularly 
enhanced in the *P,; state of the ‘‘neutral”’ 
theory, accounts for the diminution of the cross 
section. For comparison, we include the total 
cross section computed by treating the disso- 
ciated particles as unperturbed: 0.774X10-* 
cm*, and that obtained by describing the 
interaction in the P state by an ordinary 
attractive potential : 0.529 cm’. 

A specific characteristic of the spin-spin forces 
is the modification of the usual electric dipole 
sin? # angular distribution by the introduction 
of a spherically symmetric term. Proceeding 
from Eq. (28), the angular distributions have 
been evaluated for the three theories: 


(I) o(8)~sin? }+0.015 
El. dip. (II) o(#)~sin? 8+0.077 (34) 
(III) o(8)~sin? }+0.36. 


An appreciable isotropic term will occur when 
there is a large disparity in the values of the 
radial integrals Io, J1, J2; a strong attractive 
force in one of the states most readily produces 
such an effect. Although this situation is realized 
in the *P» state of the ‘‘charged” theory (II), 
the small statistical weight of this state makes 
the effect less important than that produced by 
the strong attraction in the *P, state of the 
“neutral” theory. 
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The interference between the 1S and 4D states 
resulting from photomagnetic absorption to the 
singlet state produces an appreciable departure 
from spherical symmetry. The angular distribu- 
tion predicted by the modified Eq. (79) of NPI: 


Mag. dip. (1-30) : o(8)~ (1—0.57 cos? 8) 


represents a significant reduction® in the forward 
intensity. However the total cross section for 
. the process: 


Mag. dip. (1-0): ¢=0.95 X cm?, 


is at most three percent of the electric dipole 
cross section. 

Although the cross sections for magnetic 
dipole and electric quadripole transitions to the 
triplet state are negligibly small, viz. : 


Mag. dip. (1-+0) : ¢=0.048 10-* cm? 
El. quad. o=0.074X 10-** cm?, 


it is of interest that the magnetic process, a 
consequence of the small *D,; component of the 
ground state, is comparable in magnitude with 
the electric quadripole process, which originates 
primarily in the principal component, *S;, of 
the ground state. 


CONCLUSION 


The two types of measurements which may 
be made on neutron-proton scattering and on 
the photo-disintegration of the deuteron concern 
the total cross sections, and the angular distribu- 
tions of the scattered or disintegrated particles. 
It is clear from the preceding discussion (cf. 
(23) and (24)) that a total cross section measure- 
ment of not too great accuracy can distinguish 
between the “neutral” theory and the two 
theories I and II. However, to decide between 
the latter theories by total cross section measure- 
ments would require difficult, precise, observa- 
tions. We wish to stress that valuable informa- 
tion can more easily be obtained by angular 
distribution studies, in particular, by the investi- 
gation of the angular distribution of the photo- 
disintegration particles. The angular distribution 


§ It may be noted that the direction of the deviation from 
spherical symmetry involves the sign of the deuteron 
quadripole moment. If the moment were negative, the 
intensity would be increased in the forward direction. 


of the protons recoiling under neutron impact 
will also differentiate easily between theories III 
and I, II (cf. (24)), but again the predictions of 
the latter theories are too similar to be resolved 
by anything short of a painstaking measurement. 
Although the situation would appear to be 
identical for the angular distributions describing 
photo-disintegration (Eq. (34)), it is relatively 
simpler to perform and analyze such experiments, 
for the energy of the disintegration products is 
essentially independent of angle, in contrast 
with the great variation occurring in neutron- 
proton scattering. Furthermore, photo-disinte- 
gration angular measurements test not only the 
isotopic spin dependence of the interaction, but 
in addition the presence of the spin-spin coupling 
energy, for the existence of photoelectric emission 
in the forward direction is a consequence of only 
such noncentral forces. 

We thank Professor J. R. Oppenheimer for 
his interest in this work. 


ApPENDIXx I 


The discussions of this paper and our previous paper 
NPI have been based, mathematically, upon special spin 
techniques which are peculiarly adapted to the treatment 
of states with zero or unit angular momentum, but which 
cannot be extended to states of higher angular momenta. 
We shall sketch a general method;which employs the oper- 
ational methods developed for the study of the states of a 
mesotron with unit spin.* The basic idea is that the three 
functions required to specify the states of a system with 
unit spin may, in a suitable representation of the spin 
matrices, be regarded as the components of a space vector. 
The three types of normalized spin angular functions 
(L=J, L=J+1) now appear as vector operators acting 
on the spherical harmonic Ps": 


®,2,=[J(J+1) 

The operator S;2, when written in terms of the total spin: 

r 2 
-4, 
is easily translated into the vector notation: 
rir 


The matrix elements of S,. are obtained immediately by 


*H. C. Corben and J. Schwinger, Phys. Rev. 58, 953 
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elementary vector algebra. Thus 


(r-L=0) 


—2(J+2 J(J+1)}} 
We may note that the diagonal matrix elements Sy, 
satisfy two simple sum rules 
(1) (with the exception of J=0) 


(2) (2J+1)Syi=0. 


Further, the nondiagonal matrix element, Ny, which con- 
nects the states L + J+1, is represented in terms of either 
diagonal element by 
(3) NPH=9—(Sy, 4+1)*. 
Sum rules (1) and (3) are an immediate consequence of 
the simple algebraic property of Si2: 

Sy? =8—2S)2. 
Sum rule (2) has been mentioned briefly in the text. Its 
proof is elementary: 
2 (2J+1)Ss1= (LJm|Si2|LJm) 


= =0 


for both orbital and spin quantum number sums vanish. 


APPENDIX II 


The photo-dissociation calculations of the text are 
based upon a perturbation formula in which the final wave 
function, describing the dissociated particles, satisfies the 
boundary condition that it represent a plane wave plus a 
converging spherical wave. Although this theorem is well 
known for ordinary potentials, it requires justification 
when the potential involves spin operators in such manner 
that the interaction differs for the various states of total 
angular momentum J. 

The required proof is easily given with the aid of pro- 
jection operators Oz,’ which decompose a wave function 
associated with orbital angular momentum L into its J 
components. More precisely, the projection operators are 
defined by their action on normalized spin angular func- 
tions of total angular momentum J, F,/™: 


Such operators are easily constructed from the operator J*. 
The wave function describing scattering is then constructed 
from the plane wave: 


exp pc, r)x:", 


by isolating the J component of the L wave: O17 P1(k,r) x.” 
and replacing the radial function of free motion, gz(kr), 
by the appropriate radial function 


exp [#527 Juz/(r), (ux4(r)~sin 
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Hence 


If a time dependent perturbation 
+ 

act on the deuteron system in the initial state W;, the 
function describing the outgoing wave, obtained by the 
usual procedure, is: 
M 

x (Fim 
indicating explicitly the spin coordinate o of the outgoing 
particles. (For electric dipole absorption, V(r) =—ie 
X (2xhw)te-}r.) Upon inserting projection operators, 


Ile! 

exp me) 
, 
= [= exp F,/™(r’), 

the outgoing wave becomes 

XE 0)*, Vir Wale’ 


The functions F,/"(r,¢) are obtained from the set 
P,™1(t)x:"s(¢) by a unitary transformation, and therefore 


= 2 
™L™s 


Plt, 


by the spherical harmonic addition theorem and the 
fact that 


= 
The amplitude of the outgoing wave, viz.: 


M eikr-iBtih 


is thus in the form of a matrix element with the final wave 
function: 


I(r! 


which has the form of the scattering wave function de- 
scribing propagation in the direction of r, the direction of 
observation, with magnetic quantum number ¢, the spin 
coordinate of the outgoing wave; but with the phase factor 
exp [—16z,/]. It will, therefore, describe a plane wave and 
a converging spherical wave. It is of some interest that ¥, 
is obtained from Wscatt by the combined operations of 
time and space reflection. 
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PHYSICAL REVIEW 


The Near Infra-Red Spectrum of Water Vapor 


Part I. The Perpendicular Bands v. and 2v. 


HARALD NIELSEN 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
(Received February 10, 1941) 


The rotation-vibration bands in the water vapor spectrum lying in the bolometric infra-red 
have been completely remeasured. In this paper the measurements are reported on the bands 
v2 and 2y2. An analysis of the complicated rotational structure has been carried out. From the 
interpretation of the data the value of the normal frequency »2 and the anharmonic constant 
U22 have been calculated as well as the values of the moments of inertia effective in these 
two vibration states. The frequency v2 is found to take the value 1653.8 cm™ and the anhar- 
monic constant U2,;=— 19.25 cm~'. For the states vz and 2», the effective moments of inertia 
take the values J,“ =0.905 g J,“ =1.902 g cm?, J, =3.054 g cm? 
and J, =0.775 X g cm?, J, = 1.895 X 10-™ g cm?, J, =3.095 x 10- g respectively. 
The values arrived at for the above constants serve to confirm the theory of the water vapor 


molecule. 


I. INTRODUCTION 


HE infra-red spectrum of water vapor has 

probably been studied more extensively 
than that of any other molecule. High dispersion 
measurements had been made by Sleator! 
already in 1918 and by Sleator and Phelps? in 
1925 and while their resolution was insufficient 
to permit a rotational analysis, they did serve 
to reveal the character of the bands and facili- 
tated their identification. So many bands have 
been measured and identified in the water vapor 
spectrum that after a detailed study of the 
vibration problem, Bonner* was able to evaluate 
the harmonic and anharmonic constants of the 
potential energy of the molecule with consider- 
able accuracy. An attempt to analyze the 
rotational structure was first made by Mecke* 
who successfully identified many lines in the 
bolometric bands remeasured by Plyler and 
Sleator® and in the bands in the photographic 
infra-red. Mecke did not compensate for the 
centrifugal stretching of the molecule in his 
work so that his assignment of lines can be 
considered reliable only in cases when the 
transitions are between levels where the quantum 


1W. W. Sleator, Astrophys. J. 48, 125 (1918). 
(1925) W. Sleator and E. R. Phelps, Astrophys. J. 62, 28 
oF. S. Bonner, Phys. Rev. 46, 485 (1934). 
*R. Mecke, Zeits. f. Physik 81, 313 (1933). 
(19st) K. Plyler and W. W. Sleator, Phys. Rev. 37, 1433 
1). 


number of total angular momentum, J, has 
small values. 

The most convincing progress toward a satis- 
factory interpretation of the water vapor spec- 
trum is contained in two recent papers, one by 
Randall, Dennison, Ginsburg and Weber® and 
the other by Darling and Dennison.’ The first 
of these deals with high dispersion measurements 
on the pure rotation spectrum of water vapor. 
The authors have identified nearly all the 
observed lines with quantum transitions made 
by the molecule between the rotation energy 
levels associated with the normal vibration state. 
It has thus been possible for them to deduce 
the values of rotation energy levels in this 
vibration state. The second paper is essentially 
a re-evaluation of the vibration-rotation data 
for the molecule which has enabled the writers 
to assign values to the moments of inertia 
effective in many vibration states and to evaluate 
the constants of the potential energy function. 
In general, the results are consistent with the 
theory for nonlinear triatomic molecules as first 
set forth by Shaffer and Nielsen* and subse- 
quently verified by these authors, that the 
effective moments of inertia are, to the second 


* H. M. Randall, D. M. Dennison, N. Ginsburg and L. R. 
Weber, Phys. Rev. 52, 160 (1937). 
7B. T. Darling and D. M. Dennison, Phys. Rev. 55, 128 


(1940). 
( 909) H. Shaffer and H. H. Nielsen, Phys. Rev. 56, 188 
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77 


Fic. 1. Rotation vibration lines in the 6.264 water vapor band. 


GALVANOMETER DEFLECTIONS ———>— 


order of approximation, linear relations of the’ *for this is that, when measurements have been 


vibration quantum numbers. 

The bands for which the rotational analysis 
has been least convincing are the ones lying in 
the region from 1.5y to 7.0u. A principal reason 


made, at least in the case of the fundamental 
bands, far too much water vapor has been 
present. This condition has prevailed because it 
has been necessary to make the measurements 
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TABLE I. Rotation-vibration lines in the 6.26 water vapor band. 


Line Tran- | Line Tran- | Line Tran- | Line TRan- 
No. » Oss. Cate. | No. Oss. Cate. sition | No. Ops.» Cate. sition | No. » Oss. Cate. SITION 


200(a) 1490.0 1490.0 15405 15404 16383 16378) (31-2 
(6) 1491.1 1491.1 179 1638.1} —3: | 119(a) 1696.0 1696.0 \5s—4-s 
199 14962 14961 43—5.5 15424 44—4.2 1638.5} (b) 16964 1696.6 6; 
198(a) 14992 14992 |178 1543.6 15435  5a—h |1 1641.2  1642.2(W) | 118 16978 1607.7 
14994° 15002 64—64/177 1545.3 | 1643.7 17005 (17001 
197 1502.1 ? | 176(a) 15498 1549.9 | 139 16463 16462 % —3y 1700.5 
1 1505.8 1550.3 1388 1648.1. 16482 
194 15086 1508.6 4y—5.g| ©) 15583 15583 1653.0} (6s—5 | 1706.7 
193 1509.8 () 1558.9 15587) 4; (6) 16535* 16533 2a—11/113 17108 17107 
15126 5s —5s 1558.9} (c) 16539° 16534 1710.8} \5: 
15129 6; —6, |173 1550.6 1550.6 | (d) 16544~ 1654.5 54 —G2 | 11%(e) 1713.6 
192 15124 18122 |172 15603 15605 | 185(a) 1654.9 17146 17146 6 
15128 63—61/171 15651 15653 le —2s| (6) 16554 16554 5 —h |111 1715.80 
15126  5a—4s | 169 1569.0 1569.2 2% —31|133 1661.7°(W)16002 61—5s | 110(a) 17178 17178) 
191 1514.9 6s — 6, | 168 15704 | 132 16682 16831 % —24 17178} \61—5s 
190 1516.5 1516.5 —5: | 165(a) 15748 129 1669.0 16687) f62—64| (6) 17189 17189 3 —26 
189 15174 15170 4 | (6) 15762 15763 1660.0} \42—30 | 106 1724.2 
15178  3-2—4-4 | 162(a) 1589.6 128 1669.6 1669.7) 104 17304 
188 1520.1 15200 | (6) 15004 15906 42 —5s 1669.9} \5a—4: | 108(a) 1733.6 17828 
15203 | 161 1592.3 1582.0 |127 16715 «1671.9 —3-2| (6) 17346 17848 
187(a) 1521.2 1821.0 43—44 1592.3} \3: —40 | 126 1675.6 1675.4 f44—3-2|102 17378 
15212 5a—53/160 1504.7 1508 —4s 1675.6; \5a—5-s| 101 «1740.2 1740 1) (3 
1523.1 «1528.9 Be 1596.5 1680.1  1680.0\ 1740.2; \62 
15232 29—33/159 1596.7 9 40 —54| (b) 16808 1680.0f \5.—6«|100 17432 
186 1525.6 «18256 | 158 1601.6 1601.6 —3-2| 1681.6 9 1744.4 > 
185 1527.2 |157 1608.7 16036 32-2 |123 16885 16834 41-4. | 98 17460 17460 —h 
184 15287 15288 8 |156 16078 16078 3% —41 1683.9) (5 —5a| 97 1747.5 
183 15315 15323 44-3: |155 16102 16103 —61|122(a) 1684.6" 1684.6} 443-33 1748.9) 
15313 60 —6: |152 16174 16174 h —1a 1684.3} 95 1740.6 1740.7} 15: 
182(a) 1583.2 15334 32-3 |150 16233 162933 4: —h | 1686.1 1686.1) 1740.7) —4: 
(6) 1533.9 15840 3 | 149° 16942 16242 16863} —6-:| 93 1751.6 17519 6: 
(c) 1535.5 18851 64—6-2 16283) {44 | 121 1689.9 1600.7 —7 | 91 1787.1 
1538 42/148 1628.3. 16283} {2-s—1o | 120 1690.9 1690.9) [42 —4e | 90 1758.7 
181 1539.2 }1a—24 1628.6} {66 | 17021 
1539.1} (6-2—5o 16359 


* Lines marked with an asterisk are components of a group of very close-lying lines. In such cases the frequency positions of the individual components may be slightly 


in error. 
(W) The character (W) after a line indicates that it is a line of small intensity. In certain such cases the actual frequency position may be slightly in error. 


TABLE II. Rotation-vibration lines in the 3.34 water vapor band. 


Ling TRAN- Ling TRAN- Line TRAN- Line TRAN- 
No. » Oss. » Cate §SITION No. » Oss. »Catc.  sITION No. » Oss. » No. » Oss. »Catc. 
103-0 30115 30112 3 —4 3082.2 30823 82-h 3158.2 3157.9 ? 63 3266.7 44 
—a 
102-a 3016.8 ? 30048 300.8 31788 31784) | 61 32775 32775 —h 
102-5 3017.7 3017.6 6s—7s 6s —7s 3178.9) \22 —3-2 32776 
102-¢ 30188 3018.7 6-7-7 | 90-5 3096.2 3096.2 80 3186.5 3186.1 3279.7 
10l-a 30237 30233 3: —4. | 3006.9 3007.0) 3197.2) | 3281.5 3281.1 
101-d 3028.6 5038.3) —5s 31023 31023 | 79-5 3199.0 31986 | 57 32026  32022| (4: —4e 
100 3082.2 30828 | 3105.6 3105.6 | 56 32042 32042 
98-a 3035.8 Oo 73-a 3220.7 32204) | 33054 3305.2 3: 
98-< 3038.9 ? 3119.0 3119.1 4: 72 3228.5 63-5 3313.2 ‘ 
98-d 3040.3 ? 85-b 3119.8 {4s 7l-a 3230.8" 32312 42—30 52 33143 $3143) {3 —24 
98-e 3043.9 5a 3120.1f \20 —22 7i-b 3231.1 3231.1 32-22 - 3314.8) \64 —7o 
“30674 82-4 31413 ? 65 3258.2 32585 44-45 | 47-0 3346.7 ? 
3077.6 82-e 31435 31433 6: —6s 3261.7 | 47-b 33484 ? 
93-5 30788 3079.2 32-4. | 31456 31466 2% —3. | 64 3261.6 3261.7 43 —h 
92-a 3080.9 } 82-9 3149.7 3149.8 53 —5s 3261.8 
3080.8) 
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by recording the absorptions due to the water 
vapor present in the optical path of the spec- 
trometer itself. In the case of the less intense 
bands, it has apparently been impossible to 
operate the spectrometer with slits sufficiently 


narrow to give the desired definition and 
resolution. 

It has been pointed out elsewhere® that it has 
been found practicable to surround the optical 
parts of the spectrometer in this laboratory with 
an airtight box so that the water vapor inside 
can be reduced to a mere fraction of the amount 
normally present. Under these conditions the 
lines can be almost completely resolved, a fact 
which seems to justify a complete remeasurement 
of the water vapor bands which lie in this region. 
These measurements, reported on earlier,'® are 
now complete and in this paper, Part I, we wish 
to discuss our results on the fundamental band 
v2 and its overtone 2v2, both of which are due to 
oscillations of the electric moment parallel to 
the axis of the intermediate moment of inertia. 
Such oscillations are sometimes spoken of as 
perpendicular oscillations because they are 
perpendicular to both the axes of the largest 
and the smallest moments of inertia which in 
two different limiting cases may be thought of 
as approximate symmetry axes. We shall defer 
to Part II the discussion of our measurements on 
the bands originating with alternations of the 
electric moment parallel to the axis of the 
smallest moment of inertia. 


II. EXPERIMENTAL 


The spectrometer used to make these measure- 
ments is the same used in similar experiments 
carried out in this laboratory. As in the works 
by Sleator and his collaborators the water vapor 
in the optical path of the spectrometer itself 
was used as the absorbing layer. To obtain any 
resolution whatever of the lines in the funda- 
mental band v2 which lies near 6.25y, it was 
necessary to dry the air in the spectrometer 
thoroughly. This was accomplished by placing 
large glass trays containing P.O; in the bottom 
of the box which surrounds the optical parts of 
the spectrometer. The box is made airtight by 
snapping into place the top of the box which is 
fashioned much like a refrigerator door. The 
instrument was allowed to stand from fifteen to 


® D. Cameron, W. C. Sears and H. H. Nielsen, J. Chem. 
Phys. 7, 994 (1939). 
10H. H. Nielsen, Phys. Rev. 55, 346A (1940). 
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twenty hours for equilibrium to be reached 
inside before measurements actually were begun. 
One such charge of PO; was found sufficient to 
maintain the equilibrium for a week or more. 
In spite of these precautions there was still too 
much water vapor present, the peaks of the most 
intense lines causing absorption of nearly 100 
percent. 

For the measurements on the fundamental, 
an echellette grating ruled with 2000 lines per 
inch was used. The spectrometer could be 
operated successfully with slits equivalent to as 
little as 0.3 cm~' and data were recorded at 
settings on the spectrometer circle separated by 
intervals equivalent to about 0.2 cm~'. The 
data were recorded as galvanometer deflections. 
In order to obtain as accurate a representation 
as possible of the relative intensities of the lines 
throughout the band the foreprism of the 
instrument was set carefully beforehand so that 
the galvanometer deflections at settings between 
two absorption peaks were nearly the same at 
the beginning and at the end. The composite 
picture of the data obtained from many individual 
runs over the band is shown in Fig. 1(a) and 
the frequency positions of the lines, identified 
for convenience by the arbitrary numbers 
adopted by Plyler and Sleator, are given in the 
second column of Table I. 

It was unnecessary to dry the air in the 
spectrometer to gather the data on the overtone 
band 2v2. The measurements were made with 
an echellette grating ruled with 4800 lines per 
inch. The spectrometer slits were narrowed so 
that they were equivalent to about 0.4 cm. 


(K | 


Galvanometer deflections were recorded at 
intervals on the circle equivalent to about 0.3 
cm~'. The composite absorption pattern of the 
band 2»: is shown in Fig. 2(a) and in the second 
column of Table II the frequencies of the lines 
of the band will be found, identified again by 
the arbitrary numbers assigned to them by 
Plyler and Sleator. 

In both of these curves, far better resolution 
of the lines has been achieved and much more 
detail has been observed. It is felt that until 
vacuum instruments are constructed for this 
region, little further improvement in the resolu- 
tion of the water vapor bands can be hoped for. 


III. IDENTIFICATION OF THE OBSERVED LINES 


In this section we shall present our identifica- 
tion of the lines measured in the two bands v2 
and 2v2 and from them deduce the positions of 
the vibration-rotation energy levels for these 
states. In this connection it is useful to review 
briefly some facts concerning the theory of the 
vibration-rotation spectra of such molecules and 
for our purpose we shall refer to the formulation 
of it presented by Shaffer and Nielsen." In 
their notation the rotation energies for given 
values of the vibration quantum numbers Vj, 
V. and V; and the quantum number of total 


angular momentum J, are given by the relation: . 


E(V, J, + B(V)) +01] 
(1) 


where e¢, is one of the 2/+1 roots of the secular 
determinant containing the elements 


(K|K+2)={(f—K(K+1) ][f—(K+1)(K+2)]} (2) 
[f—(K+1)(K+2) ][f—(K+2)(K+3) J[f—(K+3)(K+4) 


in which 
+J(J+1)p3 and 


pa], 


"In this paper the following typographical errors have 
been detected: (1) the sign preceding the third term in Eq. 
(35) should be negative; (2) in the last term of the equation 
defining Yo, the sign preceding the coefficient of cos y/w:? 
should be negative; (3) in the second line of the definition 
of Zo, h/ Co? should be replaced by h/Co; (4) in the equation 


A(V), B(V) and C(V) being respectively h/8*c 
times the reciprocals of the moments of inertia 


defining Y, the w: outside the braces should be w; and the 
sign before the second term should be +; (5) in the 
equation defining Y2 the sign before the second term inside 
the braces should be — ; (6) the subscripts associated with 
the quantities R; in (41) and (42) should run from i=1 to 
i=6 and not from 1=0 to i=5; (7) in the definition of k,’ 
in relations (10) the quantity (K,/2) should be replaced, by 
(K,/4); (8) in the definition of k,’ in (10) cos* a should be 
replaced by sin? a. 
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(I,), and (J,)"” effective in the given 
vibration state. The quantities p,; and Rs, Rs 
and R, which occur in (1) and (2) are coefficients 
of the centrifugal distortion of the molecule. 
They are in general complicated expressions 
depending upon the normal frequencies and the 
size and shape of the molecule and their defini- 
tions may be obtained by reference to the work 
by Shaffer and Nielsen. The symbol f is intro- 
duced for J(J+1). 

The effective reciprocals of inertia depend in 
a linear manner upon the vibration quantum 
number 


h 


h 3 
h 3 


where the X,Y; and Z; are quantities depending 
in an involved way upon the potential energy 
constants. The moments of inertia themselves 
obey a law similar to that for A(V), B(V) and 
C(V) and to this approximation they are also 
linear in V;. Concerning these Darling and 
Dennison have pointed out that there exists a 
"relation" 


I, —1,-—I, =A, (4) 


where A is independent of the anharmonic 
potential energy constants and may therefore 
be calculated from a knowledge only of the 
normal frequencies and the size and shape of 
the’ molecule. 

The theoretical ‘evaluation of the roots ¢, is 
done by expanding the secular determinant 


J=0 
J=1 
J=2 ¢=[4R.+16Rs—Ref(f—2)], 


whose elements are given in (2). This expansion 
is considerably more complicated than that of 
the ordinary asymmetrical top secular determi- 
nant“ due to the extra (K|K+4) elements 
arising from the inclusion of centrifugal stretch- 
ing terms in the energy expression. The actual 
expansion of the determinant is facilitated by 
the method suggested in the paper by Shaffer 
and Nielsen. Since the elements are diagonal in 
J and M, the determinant factors into steps in 
each of which J and M are constant; the non- 
dependence of the elements (2) on M shows 
that the steps are identical for all the possible 
2J+1 values of M associated with a given J, 
hence only one step needs to be expanded for 
each J value. Each of the J steps can be factored 
into two substeps, one associated with even and 
the other with odd values of K; thus for each J 
there is a substep of order J and one of order 
J+1. It is interesting to note that when expressed 
in terms of the R;, the (J+1)th order step 
associated with a given J is formally identical 
with the J’th order substep associated with 
J'=J+1; this fact reduces considerably the 
work required to evaluate the roots ¢, since 
only one of the substeps associated with a given 
J needs to be expanded because the other one 
has been expanded already for the preceding 
value of J. The substeps can be further factored 
by taking plus and minus linear combinations of 
the rows and columns associated with a given K. 
For our purpose it has been useful to expand 
the determinant for all values of J up through 
J=6, and the algebraic relations for the ¢, are 
given below. For each value of J only the 
equations arising from the substep of order J+1 
are given, those due to the substep of order J 
being identical with the one occurring for the 
preceding J value: 


To conform with a notation used by many molecular spectroscopists A(V), B(V) and C(V) has been chosen to 
t h/8x*c times the reciprocals of the effective moments of inertia. In the notation of Shaffer and Nielsen the 


re 
effective moments of inertia are A,, B, and C,. To avoid confusion we have replaced this notation respectively by J,, 


I, and I,. 


13 This relation may be verified by using for J,, J, and J, the quantities A(V), B(V) and C(V) expanded to the 


power —1, retaining only that term linear in (V;+}4). 
4S. C. Wang, Phys. Rev. 34, 243 (1929). 
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4 
roots of —[20R2+272R3+f(f—2)Re Je+[64(R2+ 16Rs)(R2+4Rs+f(f —2)Rs/4) 
— Je 
lf(f—2)(f—6)(f—12) =0 
J=5 &—[35R2+707Rs¥ (Ri t+2Rs)f 
(Ri t+2Rs)) — (f—12)(f—20)(Ri+34Rs)* 
— (f—2)(f—6)(Ri+10Rs+ fRe)*— (f—2)(f—6)(f—12)(f—20) Re? Je 
(Rs+2Rs)f) —25(f—2)(f—6)(R2+25Rs) 
X (Ri +10Rs+fRe)? — (Rit+2Rs)f) 
Re?) ]=0, 
J=6 &—[56R2+1568Rs— 
f(f—2)Rs) — (f—6)(f—12)(Ra+20Rs)? 
— (f—20)(f—30)(Rs+52Rs)?— (f—6)(f—12)(f—20)(f—30)Re? Je 
f(f—2)Re) —36(f—6)(f—12) 
X f(f—2)Re) 
+[576(R2+16Rs)(R2+36R3) +4(13R2+388R3) 
X — 
— (f—6)(f—12)(f—20)(f—30) Re? (f—6)(f—12) Re?) Je? 
— —36(f —6)(f—12) 
X (R2+36Rs)(Ri+20Rs)*— (f — 20) 
16(f —6)(f—12) 
X (f—20)(f—30)(R2+16Rs) Re? 
} 
X (R2+36Rs) —4f(f—2)(f—6) 
X Re? + 144f(f—2)(f—6)(f—12) 
X (Ri Re? 
X f(f—2)(f—6)(f—12) +2Ret =0. 


J 


(5) 
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When the centrifugal distortion coefficients pi, 
R;, Rs and R, have been computed and inserted 
into Eqs. (1) and (5) inspection will show that 
it is possible to evaluate the moments of inertia 
effective in a given vibration state if it is possible 
to identify enough lines in the band so as to 
determine the rotational energies for J/=1 and 
J=2. In principle then these constants can be 
inserted into the above equations and used to 
compute the energies for higher values of J. 
This should then facilitate the further identifica- 
tion of other lines in the spectrum. When this 
is tried on the energies given by Randall, 
Dennison, Ginsburg and Weber (RDGW) for 
the molecule in the normal state it is found, 
however, that the energies for higher J values 
computed in this manner are not in extremely 
good agreement with the values found experi- 
mentally. The explanation for this is given by 
RDGW who point out that the centrifugal 
distortion is not adequately taken into account 


by the terms quartic in the rotational quantum 


numbers, but that these represent them only in 
first approximation. 

We have, therefore, proceeded in a somewhat 
‘different manner. We have assumed that the 
rotation energies in the state V;= V2=V;=0 
are known more accurately than those for any 
other vibration state. We have equated the 
experimental values of the energies to the 
corresponding theoretical energies up through 
J=5, using only those, however, where e can be 
expressed explicitly. In this manner we have 
solved for the values of the effective moments of 
inertia and the centrifugal stretching coefficients 
which reproduce the experimental data as a 
whole most faithfully. These quantities will, in 
general, not satisfy the relations which the 
theory demands and no real physical significance 
can be attached to them. Their only justification 


TABLE III. Observed and calculated values of pi, Rs, Rs 


and Rs. 


OBSERVED VALUE CALCULATED VALUE 


pi 1.497 x 1.0495 x 10-2 
p2 4.73 X10? 1.503 x10"? 
pa 4.5 4.871 
R; —3.028x 10? —5.699 x10-3 
Rs 95 10-3 1.67 
Re —6.6 X10 —6.75 


is that they reproduce the experimental data 
quite faithfully over a large range of J values 
and are therefore useful in extrapolating to 
considerable accuracy the energies for still higher 
values of J. It is found, however, that the 
moments of inertia effective for the normal 
vibration state determined in this manner have 
only a small percent deviation from those 
determined by Darling and Dennison by the 
more rigorous method. The percent deviation 
between the centrifugal stretching coefficients, 
which are all small, determined by this method 
and those calculated are often large. They are, 
nevertheless, always of the same order of 
magnitude and have the right algebraic signs. 
For the sake of comparison we insert Table III 
which gives the calculated and observed values 
of the coefficients of centrifugal distortion. Since 
these coefficients are independent of the vibration 
quantum numbers, we shall assume that the 
values determined from the data for the normal 
state are also valid for the states vz and 2v». 


The fundamental band v, 


Each line in this band is due to the molecules 
making a transition from a rotation level in the 
normal vibration state to a rotation level in the 
state v2 and its frequency will be the difference 
between a term value in v2 and one in the 
normal state. The frequency v2 is due to an 
oscillating electric moment parallel to the axis 
of the intermediate moment of inertia. The 
selection rules for this case have been conveni- 
ently formulated by Dennison'® and need not 
be restated here. By trial and error enough lines 
are identified so that the term values for the 
quantum numbers J=0, J=1 and J=2 can be 
established. As a check that the identifications 
are right, combination relations are employed. 
For example, we have identified the following 
lines 19—2_2= 1565.1, 19—22= 1499.1 and 
=1635.2. Using for the term values in the 
normal state those given by RDGW we obtain 
for the rotation level 19 in cm~! the values 
1635.1, 1635.2 and 1635.2, all of which are in 
good agreement. By identifying the experi- 
mentally determined term values for J/=0, J=1 
and J=2 with the corresponding of the Eq. (1), 


% TD). M. Dennison, Rev. Mod. Phys. 3, 280-345 (1931). 
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TABLE IV. Rotation-vibration energies for the states vz and 


2v2 in the water vapor molecule. 
ENERGIES IN CmM™! ENERGIES IN ENERGIES IN 
FOR STATE cm™~! STATE cm™~! STATE 
Jr vez 2v2 Jr v2 2ve Jr v2 2v2 


Oo 1595.4 3152.3) 44 2132.2 3770.6|5-s 2000.8 3574.7 


using for the coefficients of centrifugal stretching 
those determined above we may obtain new 
values for R,, R2 and R, which will be effective 
for this state. With the values of R,(v2), Re(v2) 
and R,(v2) so determined it is possible to compute 
the values of the rotational energies for higher 
values of J and from these predict where 
additional lines are to be expected in the spec- 
trum. Once these new lines have been identified 
definitely we may make use of them to improve 
upon the values of R;j(v2) and then proceed to 
identify additional lines. It was found that this 
method was extremely useful in identifying the 
lines in the band ye. To insure that the energy 
levels were correctly assigned, combination 
relations similar to those referred to above have 
been used consistently and in all cases the 
agreement has been gratifyingly good. Lines in 
the spectrum originating from transitions be- 
tween levels up to and including J=6 have been 
identified and it is found that nearly all the 
important lines on either side of the center of 
the band to a distance of about 100 cm~' have 
been accounted for. The identifications of the 
lines in the band v2 are given in the fourth 
column of Table I. The values of the rotation- 
vibration energies for the state vz as deduced 
from these identifications are given in the first 
column of Table IV for the quantum numbers J 
up through J=6. From these the positions of 
the lines have been recomputed and for the sake 
of comparison are given in the third column of 
Table I. 


As a further verification that the identifica- 
tions of the lines in the band are correct, their 
intensities have been computed and compared 
with those observed. The calculation of the 
intensities is simple in principle, but for a 
molecule like water vapor, it becomes a compli- 
cated procedure in practice. This is because the 
molecule is an asymmetric top for which the 
eigenfunctions x(J, 7, which are necessary 
to compute the matrix components of the 
electric moment, are linear combinations of the 
wave functions ¥(J,K, M) for a symmetric 
top. These may be written: 


x(J, 7, K, M), (6) 


where the coefficients Cy,,,4 are the normalized 
subminors of the secular determinant | (K | /7| K’) 
—Eéxx:|. The evaluation of the coefficients 
C,;,,, is very laborious and to be at all rigorous 
it would be necessary to evaluate them not only 
for the state v2, but for the normal state as well, 
since the effective moments of inertia, and 
therefore also the effective asymmetry of the 
molecule, are different in the two states. This 
was actually done for the lowest transitions, but 
for the remaining transitions the approximation 
method discussed by RDGW has been employed. 
How satisfactory this method is may be subject 
to some discussion, but it is found that in general 
the intensities computed in this manner agree 
well with those observed. Actual quantitative 
agreement could hardly be hoped for under any 
circumstances since even after drying the air 
the most intense lines represent nearly 100 
percent absorption and are still quite broad. In 
order to make an actual comparison with the 
observed spectrum, we have plotted a spectrum 
as predicted theoretically using as term values 
those given by RDGW for the normal state and 
those in Table IV for the state v2. Each line is 
drawn as a triangle the base of which is roughly 
twice that of the spectrometer slit and the 
altitude of which is proportional to the calculated 
intensity. In certain cases the lines lie very close 
together in groups. In such cases a more satis- 
factory comparison is made by drawing in the 
composite envelope of the components. This 
band is shown in Fig. 1(b), below the curve 
observed experimentally. We feel that the 


~~ 
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TaBLE V. Theoretical and observed values of I,®, I,, 
I, and 


1,0 Iz) yo 1, A A 
(Vi, V2, Vs) (Ows.) (Cauc.) (Oss.) (Carc.) (Ops.) (Care.) (Ops.) (Catc.) 


1.907 3.054 3.055 0.247 0.241 
1.892 3.095 3.092 0425 OAll 


agreement between the two is extremely satis- 
factory. A few lines remain which have not been 
identified, but in no case has a line of considerable 
intensity been predicted where none is observed. 
It seems reasonable to ascribe the lines which 
remain unidentified as due to transitions between 
levels of still higher J values. 


The overtone band 2v. 


The overtone band is due also to an oscillating 
electric moment parallel to the axis of the 
intermediate moment of inertia and the selection 
rules are therefore entirely the same as before. 
Our procedure in this case has therefore been 
entirely the same as before. By trial and error, 
enough lines were identified to obtain temporary 
values for the R,(2v2). With these, lines in the 
spectrum were predicted and in most cases these 
were found to be sufficiently close to experi- 
mentally measured lines to aid materially in the 
analysis. Again the identifications were verified 
by using combination relations. The identifica- 
tions assigned to the lines in 272 are given in the 
fourth column of Table II and from these the 
rotation levels for this state have been deduced 
up through J=6. These are given in the third 
column of Table IV. Again the positions of the 
lines in the band have been recomputed from 
the term values and are set down for comparison 
in the third column of Table II. A theoretical 
spectrum has also here been plotted for com- 
parison with the experimental curve and will be 
found in Fig. 2(b). In general, the agreement is 
quite satisfactory. More lines remain unidentified 
in this case than in v2, but most of these can, 
no doubt, be attributed to transitions between 
energies of higher J values. Because of the 
nearness of 2v2 to the fundamental v3, it does 
not seem entirely impossible that some of the 
lines remaining unidentified are lines belonging 
to that band. 


IV. EVALUATION OF THE CONSTANTS OF 
THE MOLECULE 


It is now possible on the basis of what has 
gone before to calculate certain new constants 
of the molecule and verify the values of others 
already computed by Darling and Dennison. 
We shall first consider the constants occurring 
in the vibration energy. This may be written 


E(VYiV2 Va) /he= 
+ (Vit3)(Vit3) Vix, (7) 


where »; are the normal frequencies and the U;, 
are relations depending in an involved manner 
on the potential energy constants. Their defini- 
tions for the rotating-vibrating molecule have 
been given by Shaffer and Nielsen, and as far 
as the potential energy is concerned are entirely 
equivalent to the definitions in the earlier work 
by Bonner (in his notation X;,). The two bands 
ve and 22 arise when the molecule executes a 
vibrational transition from the normal state 
where V;= V2= Vs=0 to the states where V2=1, 
V3;=0 and where V2=2, V3;=0, re- 
spectively, their frequencies being given by the 
Bohr frequency condition. For the first of these 
we obtain v’=v2+2U22 and for the second we 
have v’’ =2v2+6U22. The vibration energies are 
in each case coincident with the rotation- 
vibration state 09 so that we have v2+2U22 
=1595.4 cm and 22.+6U22=3152.3 
These lead at once to the value U2.=—19.25 
cm-'!, in good agreement with that given by 
Darling and Dennison. Taking as the values for 
the other Ux those given by them we obtain for 
the normal frequency the value v2= 1653.8 
From their study of the vibration-rotation 
data of the water vapor molecule Darling and 
Dennison have arrived at the following numerical 
relations as most representative of the correct 
values for the effective moments of inertia: 


I, X10 = 1.0229 + (0.0213)(Vi+4) 
— (0.1010) (V2+3) +(0.0486)(Vs+3), 
1, X10" = 1.9207 + (0.0398) (Vi+4) 
— (0.0249) (V2+3)+(0.0077)(Vs+3), (8) 
X10 = 2.9436-+ (0.0611)(Vi+4) 
+ (0.0385) (V2+3) +(0.0441)(Vs+3). 
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It is of interest to compute from the above data 
what are the values of the moments of inertia 
I, I, and I, effective in the states vz and 
2v2 for comparison with the values predicted by 
the relations (8). To do this we have confined 
ourselves to the use of the rotation-vibration 
energies up through J=3 since in these states 
the centrifugal distortion is small so that the 
energy relations (1) may be regarded as fairly 
rigorous. The values }(1/J,+1/I,@), 
+1/I,)+1/2, and were 
averaged over these states, a little more weight 
being given to the values for which J=2 and 
J=3 since the constants will there be less 
affected by the inaccuracies of measurement. 
The coefficients of centrifugal distortion which 
were used were those calculated on a purely 


theoretical basis. The results of this calculation 
are to be found in Table V. 

It will be seen that the agreement between 
the calculated and observed values of the 
moments of inertia effective for v2 and 2v2 as 
well as for A is consistently better than that 
arrived at by Darling and Dennison on the basis 
of the earlier data available to them on these 
bands. The improved values for the fundamental 
band and the new ones deduced from our 
analysis of the overtone band 2v2 represent a 
satisfactory bit of verification of the theory of 
the water vapor molecule. 

I desire to express my gratitude to Dr. W. H. 
Shaffer of this laboratory who has verified 
several of the steps leading to the results herein 
contained. A grant-in-aid from the Rumford 
Fund is also acknowledged with gratefulness. 
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The periodic deviations from the Schottky line observed 
by Phipps and his collaborators and by Nottingham in 
measuring the thermionic emission of electrons from 
tungsten and tantalum can be interpreted as being due to 
the partial reflections of the electron waves on the potential 
hill at the surface of the metal. These partial reflections 
give rise to interference and thus to a periodic term in 
the transmission coefficient for the escaping electrons. 
The transmission coefficient is obtained by using the 
functions of the parabolic cylinder to establish the con- 
nection between the asymptotic expansions of the wave 
functions to the left and right of the top of the potential 


I. INTRODUCTION 


HEN electrons are emitted from a metal by 

the application of heat, or of an electric 

field, or of both simultaneously, one has the fol- 
lowing expression for the emitted electron current 


f Faw, (1) 


hill. The calculated positions of the maxima and minima 
of the deviations agree very well with the observed 
positions. In agreement with experiment, it is found that 
the positions of the maxima and minima are sensibly 
independent of temperature; and the amplitude decreases 
as the temperature is increased. The calculated amplitude 
of the deviations increases with the field, as does the 
observed amplitude. A dependence upon the work function 
of the emitting metal is obtained; however, since the 
heights of the surface potential barriers of tungsten and 
tantalum are very nearly equal, no experimental data on 
this dependence are available. 


where N(W) is the number of electrons with 
energy W normal to the emitting surface. N(W) 
is given by the Fermi distribution in the energy, 
which holds for the electrons in the metal; and 
D(W, F) is the transmission coefficient for the 
electrons incident on the potential barrier at the 
metal surface with energy W normal to the 
surface. The potential field which the emitted 
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Fic. 1. The po- 
tential barrier trav- 
ersed by the elec- 


hoo trons. The metal 

or =U; 

for x>m, V=W. 


—1/2x—x/2x0. 


electrons must traverse may be, with good ap- 
proximation, taken to be 


V =const. =0 inside the metal, (2a) 


outside the metal 
(but not right up (2b) 
to the surface), 


V=W.—e/4x—eFx 


where W, is the difference in the potential of the 
electron inside the metal and at infinity when 
F=0; F is the electric field applied to the metal 
surface ; x is the distance of the electron from the 
metal surface, which is placed at x=0; and e is 
the electronic charge. The term e?/4x is the 
“image potential’? which results from the at- 
tractive force exerted on the escaping electron by 
the induced charge +e in the metal.'~* The term 
eFx is the potential due to the external electric 
field used to attract the electrons over to the 
collector. Equation (2b) will not hold right up to 
x=0 (i.e., the metal surface), where it has an 
infinity. The true potential outside the metal 
very probably goes to zero and joins smoothly 
to the line V=0, at x=0. We may closely ap- 
proximate the true potential by modifying (2) as 
follows: 

V=const. =0 (3a) 


V=W,—e/4x—eFx x>x1, (3b) 


where x,~e?/4W, is the point at which the 
external potential goes to zero, i.e., where it 
crosses the x axis. The discontinuity which has 


1 We presume here, of course, that so long as the electron 
is some distance away from the metal surface, the surface 
may be considered smooth, so that this simple image law 
holds. Similarly, we assume one work function for the 
surface. Actually, each electron energy has its own po- 
tential barrier (reference 2) and each crystal face at the 
metal surface has its own work function (reference 3), the 
largest contributions to the current coming from the facets 
or spots of lowest work function. Agreement of experiment 
with the Schottky theory, which assumes the image law, 
justifies the first assumption. In all work of this type with 
emission phenomena, it seems reasonable to assume a 
weighted average work function. 

2 J. Bardeen, Phys. Rev. 49, 653 (1936); 58, 727 (1940); 
or cf. F. Seitz, Modern Theory of Solids, (McGraw-Hill, 
New York, 1940), p. 163 and pp. of 

3M. H. Nichols, Phys. Rev. 57, 297 (1940). 


been introduced in the potential at x; does not 


exist in nature, of course, but since in the — 


neighborhood of x; the true potential changes 
appreciably in a distance comparable to the 
de Broglie wave-length of the electrons, it may be 
accurately represented with this discontinuous 
one. The introduction of a more complicated 
potential function which would go smoothly to 
zero at x=0 would not appreciably alter our 
results, but would greatly complicate our calcu- 
lations. Also, for our present purposes we have 
neglected the periodic variation of the potential 
inside the metal ; our inner potential is an average 
over the periodicities. 

If, for convenience, the following system of 
atomic units is adopted : unit of length =a = h?/me?* 
=0.528A = radius of first Bohr orbit in hydrogen ; 
unit of energy =e?/2a = 13.54 ev=ionization po- 
tential of hydrogen, Eq. (3) becomes 


V=0 (4a) 
V=W,-1/2x—x/2x0? x>x, (4b) 


where Xo, the position of the maximum of the 
potential, is determined from the equation 
e?/4x9=eFxo. The potential of Eq. (4) is shown in 
Fig. 1. 

So long as intense (greater than 10° volts 
cm-') electric fields are not applied to the metal, 
the electrons whose energies are less than that 
corresponding to the maximum height of the 
barrier (W,’ in Fig. 1) are not pulled from the 
metal. In this case the electron current obtained 
is due to electrons whose energies are greater 
than W,’; this emission is called thermionic or 
thermionic field emission. Thus, for thermionic 
emission Eq. (1) becomes 


i(F, T)= N(W)D(W, F)dW, (5) 


Wa’ 
where W,’, the maximum barrier height at the 
field F, is given by 

= W.—1/xo. 
The term 1/xo represents the amount the barrier 
has been lowered by the field F. This gives us the 


well-known Richardson equation for thermionic 
emission : 


(6) 
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where x’=x—1/x0; x is the thermionic work 
function for F=0; i.e., x= W.— W;, where W; is 
the width of the Fermi band at T=0; Ap» is the 
emission constant; 7 is the absolute tempera- 
ture; and D(F) is the transmission coefficient 
summed over all energies in the Maxwellian 
distribution, to which the Fermi distribution 
reduces for electrons with energies greater than 
W.’. If io is the current one obtains for a given 
temperature with zero external field, then one has 


D(F) e! 


log ir —log io =log ——-+——F!. (7) 


kT 


Assuming for the moment that the transmission 
coefficient does not change appreciably with the 
field, one has the result that plotting log ip 
against F should yield a straight line; this line is 
known as the Schottky line. Equation (7) may be 
expected to hold rather exactly for fields which 
are not too small or too large, i.e., in the range 
3X10°<F<10° volts For fields below 
F=3000 volts cm~', Eq. (7) does not hold be- 
cause of ‘patchy emission.’’ For fields greater 
than F= 10° volts Eq. (1) does not reduce 
to Eq. (5); in other words, with intense fields, 
electrons whose energies are less than W,’ are 
pulled from the metal and tunnel through the 
potential barrier; so again Eq. (7) does not hold. 
In what follows we therefore expect our results to 
agree with experiment for the range of fields for 
which Eq. (7) is valid. 

Using temperatures up to 2300°K and fields up 
to 6.4X10° volts cm~', Phipps and his collabo- 
rators*® very accurately measured the emitted 
electron current from tungsten and tantalum and 
observed periodic deviations from the Schottky 
line. They found further that the amplitude of 
the deviations increases with the field. These 
results were recently verified by the precise work 
of Nottingham.* Both Phipps and Nottingham 
agree that within the limits of experimental error 
the locations of the maxima and minima are 
independent of the temperature. However, there 
is as yet some disagreement about the tempera- 
ture dependence of the amplitude ; Phipps and his 


a 938). L. E. Seifert and T. E. Phipps, Phys. Rev. 56, 652 
5 > Turnbull and T. E. Phipps, Phys. Rev. 56, 663 


(1939). 
®°W. B. Nottingham, Phys. Rev. 57, 935 (1940). 


co-workers observed that the amplitude decreases 
with increasing temperature, while Nottingham 
observed no temperature dependence. 

Mott-Smith’ suggested that the transmission 
coefficient might be periodic because of partial 
reflection of the electrons on the potential hill 
just outside the metal surface (Fig. 1), but the 
special fields and approximations which he used 
did not give the correct results, because he used a 
parabolic approximation to the potential (4) 
which is not valid over even a small fraction of 
the large range of fields used in the experiments. 

In the present investigation the W. K. B. 
method is applied for the solution of the wave 
equation to the right of x;. In order to obtain 
that solution of this equation which represents an 
outgoing (i.e., transmitted) wave on the right 
side of the barrier, ‘“‘connection formulae’’ must 
be applied. It is customary to set up these ‘‘con- 
nection formulae” by replacing the actual po- 
tential in the neighborhood of turning points by 
straight lines. This procedure is valid, however, 
only if the coefficient of y (i.c., W—V) in the 
wave equation has a zero of first order at the 
turning point. In the present case (W— V) has a 
zero of second order; to determine the ‘“‘con- 
nection formulae’”’ the actual potential is replaced 
by a parabola (instead of a straight line) in the 
neighborhood of the turning point. After “‘con- 
nection formulae” are obtained, it is possible 
to calculate both the transmission coefficient 
D(F, W) and its summation over the electrons 
D(F), in an entirely straightforward manner. 

It is to be pointed out that no arbitrary, 
adjustable parameters are introduced into the 
solution. For this reason, the rather good agree- 
ment with experiment which is obtained serves as 
another check of the fundamental assumptions 
involved in the modern theory of emission of 
electrons from metals. 

According to Eq. (7), deviations from the 
Schottky line will be obtained if D(F) changes as 
F is changed. Leaving the Schottky slope out of 
consideration, the deviations are given by 


A log i=log i—log io=log D(F)—log D(0). (8) 


7H. M. Mott-Smith, Phys. Rev. 56, 668 (1939) and 

al communication. The authors were greatly aided in 

this research by communications and conversations with 
Professor Mott-Smith. 
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The problem of finding the deviations is thus 
reduced to finding D(F). This we now proceed to 
do. 

II. THE METHOD 


We have the following Schroedinger equations 
for the wave function: 


d*p,/dx*+ Wy; =0 (9a) 
/dx? + (W— 
x>x,. (9b) 


The solution of (9a) may be written down im- 
mediately 


=a, exp ((W'x)+a2exp (—iW!x) (10) 


the two exponentials representing the incident 
and the reflected wave, respectively. Assuming 
for the moment that we can obtain a solution of 
(9b) valid in the neighborhood of x:, we can 
determine the transmission coefficient by setting 


and solving for 
F)=1—|a2/a;|?. (12) 


Since the emitted electrons have a Maxwellian 
distribution, one then obtains the energy inde- 
pendent D(F) by summing (12) over all electron 
energies, using the Maxwell distribution. 

Our problem is now to find that solution of 
(9b) valid in the neighborhood of x:;, which for 
x>>xo has the form of a transmitted wave: 


Cort exp] i f 


We may write (9b) as 


/dx*-+oy=0; wi 
W— = e+ (x —x9)?/2xx0?, 
where «=W—W,+1/xo is the height of the 
electron energy over the maximum of the po- 


tential barrier. 
The asymptotic expansion of the desired solu- 


| l 
0.02 0.03 0.04 
€ IN ELECTRON VOLTS 


Fic. 2. The transmission coefficient (22) plotted as a 
function of the height of the electron energy above the 
maximum of the barrier for three different field strengths. 
The values of ¢ are given in electron volts; the values of F 
are in volts per cm. 


tion near x, is of the form 


exp| f (14) 


The solution (14) represents a W. K. B. type of 
approximation for the wave function; it is 
certainly valid in regions where the potential is 
not a quickly varying function of x. However, we 
are interested in the validity of (13) near x, where 
the potential is a rather rapidly varying function, 
so we must use more precise validity arguments. 

Now, because of the smallness of x1, one may 
write the inequality x:/2xo?<1/2x,, so that near 
x, the condition for the validity of the solution 
(14) is the same as the condition for the validity 
of the W. K. B. solution to the problem involving 
the pure image potential, for which V= W,—1/2x. 
At x this condition may be written in the form 


Wp 5S 
We 16W 


It is seen that the condition is obeyed very well 
if W.a«xW. However, it seems reasonable to 
believe that the principal contributions to the 
periodic deviations come from electrons whose 
energies are only slightly greater than W.. If 
W=W,, the left-hand side of the inequality is 
approximately one-half for tungsten, so that for 
these electrons the condition is not very well 
fulfilled. Nevertheless, we know that the exact 
solution of the Schroedinger equation with image 
potential for W= W, yields a transmission coeffi- 
cient which differs from that obtained with the 
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W. K. B. solution by about 0.5 percent. Since we 
are interested in obtaining the transmission 
coefficient for a case which yields the same 
validity condition as the image case, it seems 
reasonable to expect that a solution of the type 
given by Eq. (14) will give a close approximation 
to the correct transmission coefficient. Con- 
sideration of higher order terms in the W. K. B. 
expansion would probably lead to a_ small 
correction to the nonperiodic term and a negli- 
gible contribution to the periodic term. We shall 
assume that the representation (14) will yield a 
transmission coefficient which is accurate enough 
to exhibit the deviations we are seeking. 


III. CALCULATION OF THE TRANSMISSION 
Cc OEFFICIENT D(F, W) 


Assuming (14), we shall determine the con- 
stants 6; and by and obtain the transmission 
coefficient D( F, W) by the use of (11) and (12). 

The constants }; and by are determined as 


II. Let 


Proceeding as in I we have for 2’>0 


We associate 


579 


follows: We expand the potential (4b) about the 
point x» and neglect terms of the third and higher 
orders. Thus, 


V= W,.—1/x0— (x 
If we put 2=(2/x9*)'(x—xo) and 
a= W.+1/x0) = (08/2) 
the Schroedinger equation becomes 
dy (15) 


This equation can be transformed into the differ- 
ential equation of the parabolic cylinder in the 
following two ways :* 
I. Let 
u=se in/4 
Then 
dy /dw+(—ia—u?/4)y=0. 


We are interested in the asymptotic expansion 
for xXxo, i.e., for z«0; so we set 2’ =ze'*. Then 
for® 2’>0 

(22) 


r'(}+ia) 


(16) 


u = ze~ **/4, 


(27) 


—ta) 


i(2’*/4+a log 2’ 


(17) 


ei (2 log 2’) with exp lif 
To 
which means 


The outgoing wave for x>x» should be of the form 


y~cexp lif 


’ The procedure used here of establishing the connection between the W. K. B. expansions for large and small x was 
first indicated by H. A. Kramers and G. P. Ittmann, Zeits. f. Physik 58, 225 ff. (1929). The procedure was extended to 
include zeros of orders higher than the second by S. Goldstein, Proc. London Math. Soc. 33, 246 (1931-32), and by R. E. 
Langer, Bull. Am. Math. Soc. XL, 545 (1934). These authors, however, give the connection formulae in a form which is not 
applicable to our problem. The formulae developed here may be applied to the problem of the H** ion; they also give a 
complete solution to the problem advanced by E. C. Kemble, Fundamental Principles of Quantum Mechanics 


(McGraw-Hill, New York, 1937), pp. 109-112. niga 
*E. T. Whittaker and G. N. Watson, Modern Analysis (Cambridge University Press, third edition, 1920), pp. 348-349. 
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which is associated with e~*(+*/4t+ les #), Tt is to be noted that for the case «=0, 
| (xo—x)/V2xox! x<x0 
(x—x0)/V2xox! x>x0. 
The solution under I possesses the asymptotic expansion for 2>0: 


co-* exp f os 


the outgoing wave. Thus Eq. (16) gives the appropriate solution. Comparing (16) with (14), we have 
by als, 


which corresponds to 


Since I'(r)'(1—r) =2/sin ar, choosing r= }+ia we have 


so that we may write 
b,=exp [—ix/8—i arg I'(3+%a) 


Our solution in the neighborhood of x; is then 


y~exp [—ix/8—i arg I'(4+ia) exp [if ous| 


exp [ -if ous], (18) 


an expansion which is valid at x, for all energies from the top of the barrier to infinity. However, the 
second term rapidly goes to zero for increasing electron energies; this means that as we consider 
electrons higher and higher above the barrier, the reflected current becomes smaller and smaller. 
Thus, with the use of the function of the parabolic cylinder we establish the connection between the 
asymptotic forms for large and small x. The validity of the W. K. B. type of approximation for this 
type of problem has been confirmed numerically. For the case when F=0, the Schroedinger equation 
may be solved exactly with the use of the confluent hypergeometric function. The transmission 
coefficients so obtained"® are in close agreement with those obtained by the W. K. B. approximation 
to the exact solution. 
We are now in a position to obtain the transmission coefficient D(F, e). 


D(F, =1-|— —-—_e-*(1 14 ) cos (0-+tan- ~) 
4g} 84! 
1 ¢”? 1 ¢” —1 
+— 4+ 
32 64 
where 


o'(x)dx—arg T(4+ia); ¢’= o= (x1). 


10 Exact calculations for this case have been made by L, A. McColl, Phys. Rev. 56, 699 (1939); cf. also L. W. Nordheim, 
Proc. Roy. Soc. 121, 626 (1928). 
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We may expand the transmission coefficient as follows: 


1 1 y 
D(F, ~ + 
1 1 (1 (i +e-272)8 
Putting in the value for y and neglecting terms which are — of no consequence, 
1 
— cos (0+ 
49! (1 


329% (1 


The coefficient of cos 2v is already very small when a=0 and goes rapidly to zero as a is increased. 
Hence, only a very small error is introduced by neglecting this term. (Actually carrying through the 
summation including this term verifies this.) It is to be noted that terms involving the factor 
e~**=exp [ — r(xo?/2)4e ] may be neglected save for very small values of ¢; hence, in these terms one 
may take 


¢’ 1 W.? (W,)! 

4g) 
Because the exponential term goes to zero so rapidly, its contribution to the transmission coefficient is 
negligible if a>1, the principal contribution coming for values of a1. Hence, one may take 


arg I'(}+ia) ~ —(y+2 log 2)a 
where y is Euler’s constant: y =0.5772. 
Writing 
71 
z0 70 


and performing the integration of the left-hand member numerically, one sees that even for the 
largest ¢ values which need be considered (because of the factor e~** we need consider only very small 
€ values), 
SOK o¥(x, 0)dx ; 
70 


furthermore, f(€) goes rapidly to zero as e—0. Consequently, for the small « values, f(e) may be neg- 


lected. Thus 
v= log 2)a. (21) 


The error introduced here in the period of the deviations is negligible, but the error in the phase may 
be a few percent. With these simplifications (20) becomes ; 


1 W.} 1 
D(F, —————- cos 


4] 


(22) 


with v given by (21). The transmission coefficient (22) is plotted in Fig. 2 as a function of the height 
of the electron energy above the barrier, for various values of the field intensity. The physical 
interpretation of the fluctuations in D(F, e) is given in the discussion. We mention only that the 
positions of the maxima and minima of the fluctuations are determined by xo, which is a measure of 
the width of the potential barrier, and by W, which represents the height of the barrier. 
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IV. SUMMATION OVER THE ELECTRON ENERGIES 


We shall now sum D(F, e) over all electron energies and thus obtain the energy independent 
transmission coefficient D(F) which is a function only of the temperature, the height of the barrier, 


- and the field. Since we are interested only in electrons for which W>W,, the Maxwellian energy 


distribution may be used, one has 
1 a 
f D(F, de. 
kT J, 


In summing over all electron energies, the last term of (22) may be dropped ; also 


4v2 2 W.} 
x (—1)"e-2"** cos } —xo! -——+tan-! ———(y+2 log 2)a;. (23) 
13 4 


The integral of this expression can be readily carried through and leads to a complicated expression 
which, to a very good approximation, may be written as 


=1 log 2 W,} xe\! 1 \2 


xo 1 
2 kT/ | 
The deviations from the Schottky line are then given by © 
A log i=log D(F) —log D(0). (25) 


It is probable that the integral of (23) does not yield a good value for the transmission when 
F—0. In the case F=0 one may make a separate calculation of the transmission, a calculation which 
has been given by several authors. However, if one uses arbitrary intercepts in plotting the deviations, 
the deviations so obtained are independent of D(0). We are not interested in the case where F is 
different from zero but small, since there the periodic deviations are masked by “‘patch-effects.”’ 


4v2 2 
“Cos } — +tan-! 
3 W,,} 4 


V. Discussion OF RESULTS 


According to Eq. (24) the transmission coeff- 
cient consists of two parts. The first term varies 
monotonically with the field, whereas the second 
term varies in a periodic manner as one changes 
the field. This periodic term may be interpreted 
physically as due to interference arising from the 
superposition of waves reflected from different 


heights on the potential hill, the different heights 
corresponding to different x-values."' For the 


"The rectangular barrier serves as an example of 
reflection of the electron waves from two plane potential 
surfaces. In this case the interference of these two reflected 
waves yields a transmission coefficient which is strongly 
periodic in the energy of the incident particles. If one uses 
the rectangular barrier as an approximation to our barrier, 
assuming that changing the field alters the width and 
height of the barrier, one obtains a transmission coefficient 
strongly periodic in both the energy and the field. Reflec- 
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pure field potential (i.e., neglecting the image 
force), x; and x» coincide and no interference 
occurs. A potential involving no discontinuities 
which exhibits this interference phenomenon is, 
for instance, that given by Eckart.” 

Using (24) one may easily calculate the posi- 
tions of the maxima and minima of the deviations. 
These are shown in Table I along with the 
observed values ;* it is seen that the agreement 
between theory and experiment is rather good. 
The positions of the maxima and minima are 
sensibly temperature independent. The period 
depends only upon xo; i.e., it is the same for all 
metals. Physically this means that the period is 
determined by the image force and external force 
only and does not depend sensibly upon the 
exact shape of the potential near x;. The ampli- 
tude, on the other hand, will depend somewhat 
more sensibly (as formula (24) shows) upon the 
shape (particularly the slope) of the potential 
near x;. To illustrate this we may superpose a 
term b/x? on the potential ; thus, V(x) = W,—1/2x 
—b/x*—x/2x,*?. This potential probably changes 
the image potential more violently than is really 


TABLE I. Calculated and observed positions of the maxima 
and minima of the deviations in terms of xo. The maxima are 
given in italics. The observed values are taken from the work 
of Turnbull and Phipps* and of Nottingham.’ The values of 
Xo are here given in angstrom units. 


OBSERVED (TUNGSTEN) CALCULATED 
NOTTINGHAM TURNBULL AND PHIPPS Ww Ta 
172.5+3.0 172.2 173.5 
146.022.2 141.9 143.1 
xo= 118.6 118.6+1.4 114.5 115.6 
94.9 94.940.9 90.1 91.1 
74.2 71.4+0.6 68.6 69.4 
54.6 54.5404 50.0 50.7 
37.6 37.640.3 34.3 34.9 


tion from the rectangular barrier is analogous to the 
reflection of light at two surfaces where one has abru 
changes of the index of refraction. Reflection from the 
barrier of Fig. 1 is analogous to reflection of light in a 
medium with continuously varying index of refraction. The 
transmission coefficient obtained for the barrier of Fig. 1 
differs in another essential respect from that obtained with 
the rectangular barrier; the periodic term does not disap- 
ad when the energy is just slightly under the barrier 

eight, as it does in the case of the rectangular barrier. 
This is similar to the analogy between the transmission 
coefficients for the image potential and for the single step 
barrier. For the image potential the transmission is not 
zero when the energy is just below the height of the 
barrier; for the single step Caster W=W, yields a trans- 
mission zero. 

® Carl Eckart, Phys. Rev. 35, 1303 (1930). 
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Fic. 3. Plots of the observed and calculated deviations 
for the two extremes of the temperature used by Turnbull 
and Phipps. The dotted curves represent the calculated 
values. in the experiments, the intercepts are chosen 
arbitrarily. 


the case. If b is chosen in such a way that the 
slope of the Schottky line is changed by only one 
percent, the period and phase of the deviations 
from the Schottky line are changed only slightly, 
whereas the amplitude may be changed by as 
much as twenty percent. 

In addition to the experiments on the Schottky 
slope, another evidence for the validity of the 
image force at distances of 10-7 cm is here 
furnished by the good fit of the period of the 
deviations with experimental values. The check 
of amplitude and phase shows further that the 
chosen shape [Eq. (3) ] is also a reasonably fair 
representation of the actual potential near x; 
(inasmuch as the dependency of the amplitude 
upon the slope of the potential is not critical). In 
principle, accurate conclusions could be drawn 
about the nature of the true potential from the 
comparison between theoretical dependency of 
the amplitude upon the shape of the potential 
near x, and the observed amplitude. Practically, 
however, this is not feasible; the concept of a 
fixed potential acting on the electrons and the 
one-dimensional treatment of the thermionic 
problem are, in all probability, good approxi- 
mations, but their precise limitations are not 
known at present. 

In Fig. 3, the dashed curves are obtained from 
(24) ; the solid curves are taken from the experi- 
ments.® As the table and the graph show, the 
agreement with experiment is good, inasmuch as 
no arbitrary parameters were introduced in our 
treatment. The W; values used in the evaluation 
of (24) were obtained from the assumption of one 
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free electron per atom. This assumption gives an 
upper limit for W;; a smaller value of W;, as 
indicated by M. F. Manning and M. I. Chodorov, 
Phys. Rev. 56, 787 (1939), would improve the 
agreement with observation, as was pointed out 
to us by Dr. C. Herring. A smaller value of W; 
also means an increase of x;, which would be 
desirable for physical reasons. 

From Eq. (24) it is seen that for low fields, i.e., 
for large xo, the temperature dependence of the 
amplitude of the deviations is given approxi- 
mately by the factor 1/kT ; however, this temper- 
ature dependence changes somewhat with the 
field. For all fields, the amplitude is increased by 
decreasing the temperature. This result is in 
agreement with the work of Phipps and collabo- 
rators, but is contrary to Nottingham’s con- 
clusion that the deviations are independent of 
temperature. However, this question is certainly 
unsettled from the experimental standpoint. 
From the theoretical viewpoint it is difficult to 
see how these deviations can be independent of 
temperature. It is hoped that further experi- 
mental work will clear up this issue. 

The dependence of the deviations upon the 
barrier height, W., of the metal is worth pointing 
out. The principal change should occur in the 
shift of the maxima and minima due to the change 
in the phase term —2/W,'!+tan-1W,!. The 
barrier height which we have called W, is the 
sum of the Fermi energy W; and the work func- 
tion x. For tungsten one has W;=5.8 ev and 
x=4.53 ev, so that W,=10.33 ev or about 0.76 
of our unit. For tantalum W;=5.2 ev and 
x =4.07 ev, so that W,=9.27 ev or about 0.68 of 
our unit. This difference is not enough to produce 
a great change in the positions of the extrema as 
one can see from Table I in which a comparison 
is made between the calculated positions of the 


extrema for tungsten and tantalum. Since no 
metal with a W, value differing much from that 
of tungsten can be used for this work on 
thermionic emission, no data are available on the 
W, dependence. However, it is hoped that the 
deviations will soon be obtained with photo- 
electron currents for metals with various values 
of W, and for thermionic emission of thoriated 
metals. 

As we pointed out in the introduction, our 


theory does not take into account any intense , 


field emission. However, if the emitting metal 
surface is not perfectly smooth, certain small 
areas of it are subjected to fields larger than the 
apparent or average applied field which one 
records. Thus at fields of 4X10° volts cm™ 
certain small areas of the emitting surface are 
subjected to fields of 10° volts cm~ or greater, 
and electrons are pulled from these areas and 
tunnel through the potential barrier of Fig. 1. 
The current due to these electrons is small for 
fields of 4X 10° volts cm™", but it is large enough 
to affect our small deviations at these fields. This 
fact explains the marked deviation (as shown by 
the curves of Fig. 2) of theory from experiment 
for fields higher than 3.5 volts It is 
very probable that one can use the parabolic 
connection to develop an expression for the 
transmission coefficient of the potential (4) valid 
for all electron energies. One would then obtain a 
theory taking all currents into account. However, 
the influence of any emission save thermionic 
field emission is so small with the fields used by 
the experimenters that the simplified theory 
presented here covers the experiments very well. 

Part II of this paper, discussing the photo- 
effect in electric fields and the transition from 
thermionic to intense field emission, will appear 
in the near future. 
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Discontinuities of Magnetostriction and Magnetization in Nickel 
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When the magnetic vector of a domain reverses direction there will be no magnetostrictive 
changes in the domain. If the vector jumps through some angle other than 180°, a sudden 
redistribution of strains will be caused by the change of magnetostriction and a sound wave 
will travel out from the domain. Simultaneous records have been secured of the sounds and 
of the magnetization jumps when the magnetization of a nickel wire changes. These records 
indicate that reversals occur mainly when the magnetization of the specimen is weak. Jumps 
through other angles than 180° occur in regions of strong as well as weak magnetization. 


HILE the intensity of magnetization of a 

ferromagnetic is changing it is supposed 
that at one stage the magnetization vectors of 
small domains are suddenly changing their direc- 
tions so as to line up more closely with the 
direction of the applied magnetic field. These 
vectors may swing through 180°, in which case 
we have a “‘reversal,”’ or they may swing through 
a smaller angle. It is probable that the angle of 
jump is governed by conditions of strain as well 
as by natural directions of easy magnetization 
in the crystal. Bozorth and Dillinger! have 
measured the transverse (to the applied mag- 
netic field) components of the small sudden 
changes of magnetic moment, as well as the 
longitudinal changes. They state that their ex- 
periments agree with the following conception: 
“When the net magnetization is small, most of 
the changes are reversals in the domains in 
which the initial magnetization is almost anti- 
parallel to the applied field. When the net 
magnetization is large . . . other domains are 
affected in which the initial magnetization is 
inclined at smaller angles to the field and changes 
of 90° occur.” 

The discovery of a transverse component of the 
Barkhausen effect does not actually prove that 
changes of 90° exist. If the magnetization vector 
of a domain initially makes an angle with the 
field and then suddenly reverses direction, the 
transverse component of the vector will also 
reverse and produce a transverse Barkhausen 
effect. There are two ways of distinguishing 
reversals from other types of changes, such as a 


1 Richard M. Bozorth and Joy F. Dillinger, Phys. Rev. 
41, 345 (1932). 


change of 90°. From a consideration of the condi- 
tions of symmetry it is obvious that the electrical 
resistance and the magneto-resistance are not 
affected by a reversal of the magnetization 
vector. On the other hand, any change in the 
vector which is different from 180° will affect 
both of these quantities. Accordingly, if the con- 
clusion of Bozorth and Dillinger, as quoted 
above, is correct, one would expect magneto- 
strictive jumps to be associated with Barkhausen 
jumps only when the net magnetization is large, 
that is, when changes of 90° occur in the mag- 
netization vectors of the domains. The experi- 
ments described below were designed to test this 
hypothesis. 

The apparatus used for detecting simultane- 
ously both jumps of magnetization and of mag- 
netostriction is shown in Fig. 1. The customary 
method was used for amplifying and recording 
the magnetic discontinuities. The oscillograph 
O; was made by coupling mechanically a light 
mirror to the armature of a Baldwin telephone 
receiver. To detect the jumps of magnetostriction 
another Baldwin receiver was used. In this type 
of receiver the mica diaphragm is connected by 
a short, thin rod to an armature which, moving 
in the field of a permanent magnet, alters the 
flux through a pick-up coil around the armature. 
The mica diaphragm was removed and the short, 
thin rod replaced by a wire about 8 cm long 
fastened to one end of the specimen as shown in 
the figure. In order to record the jumps of 
magnetostriction it was necessary to use high 
amplification and to work in a place free from 
noise and vibration. The apparatus was mounted 
on sponge rubber and transverse vibration of the 
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long connecting wire was limited by oil damping. 
The Einthoven string galvanometer O2 was 
equipped with a gold Wollaston wire tuned for 
high sensitivity. In order to secure a smooth, 
continuous variation of the magnetic field two 
rectifying tubes (RCA 5Z3) were connected in 
parallel and the current through these tubes also 
traversed the solenoid. Breaking or making the 
filament current then caused the solenoid current 
to change smoothly and slowly. The deflections 
of both O; and QO, were recorded photographically 
and simultaneously on a moving film. The 
amplifiers were adjusted so that the maximum 
deflections were about the same when a given 
specimen was taken over a certain range of the 
hysteresis loop. This range (on the left-hand, 
downward branch of the loop) corresponded to 
the removal of remanence magnetization and the 
carrying of the specimen into the saturation 
region well past the knee of the curve. 

If the diameter of a domain in the specimen 
is of the same order of magnitude as the diameter 
of the specimen itself, one would expect the 
longitudinal length change of the domain to be 
transmitted suddenly, without much alteration, 
to the armature of the receiver. However, for 
domains which are relatively small, a single 
jump of 90° would cause a local readjustment of 
strains without much length change of the 
specimen as a whole. A sound wave would travel 
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Fic. 1. Apparatus for detecting jumps of magnetization 
and magnetostriction. 
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down the wire and be detected by the receiver.? 
A reversal of 180° does not necessitate any strain 
readjustment and so should not produce much 
sound. However, a very feeble sound could be 
produced if a reversal is due to a magnetic wave 
traversing the domain. This wave may be con- 
sidered as a barrier layer about 30 atoms thick,* 
separating reversed from unreversed atoms. In 
this layer the magnetization vectors of the 
separate atoms will progressively alter directions 
until the total direction change of 180° has been 
produced. While the layer is traversing the 
domain there will be local strains in the layer 
arising from the atoms whose vectors are rotated 
through less than 180°, but the sound produced 
must be very much feebler than that produced 
when the whole domain suffers a 90° change. 
It seems safe to conclude, therefore, that no 
sound will be detected in association with a 
Barkhausen jump unless that jump is less 
than 180°. 

Figure 2 shows two portions from one of the 
several oscillograph records made for an un- 
annealed nickel wire 3 cm long, 0.02 cm in 
diameter, under a tension of 0.9 kg/cm*. The 
top section of Fig. 2 corresponds to a small part 
of the left-hand side of the hysteresis loop at the 
steepest point just where it, crosses the axis. The 
bottom section corresponds to a small part of 
this same descending branch of the loop, but out 
just beyond the knee. The Barkhausen jumps 
were here too weak to show up on the oscillogram 
O,, although weak magnetostrictive jumps are 
recorded by Oz. 

When an iron wire was investigated with this 
apparatus, the jumps of magnetostriction were 
too small to be registered by the string galva- 
nometer. With high amplification the jumps 
were of the same order of magnitude as the 
background disturbances; however, they could 
be faintly heard when a telephone receiver was 
substituted for Ox. 

Examination of a number of films for nickel, 
similar to Fig. 2, led to the following conclusions. 

2 These sounds were first observed by Arkadiew (Comptes 
rendus 184, 1233 (1927)), who merely fixed the end of a 
ferromagnetic wire to a thin diaphragm across a glass tube 
and heard sounds when a permanent magnet was moved 
near the wire. The writer could not secure satisfactory 
results in this way ry for nickel in a state of strain. 


3F. Bloch, Zeits. f. Physik 74, 295 (1932); J. L. Snoek, 
Physica 6, 161 (1939). 


] 
( 
OSCILLOGRAPH 
(02) 
| Cd 
necever | |y 


MAGNETOSTRICTION IN NICKEL 


O: 


8 8 8 8 8 8 8 8 


Fic. 2. Records of Barkhausen jumps (white line, 0, on dark background) 
and magnetostrictive jumps (dark band, O2, on white background) at two different 


places on the hysteresis loop. 


(1) There are many instances of Barkhausen 
jumps which have no associated magnetostrictive 
jumps of appreciable size, and similarly, many 
magnetostrictive jumps have no associated Bark- 
hausen jumps. (2) After the knee of the hysteresis 
loop is passed the Barkhausen jumps become too 
small to be recorded by the oscillograph, although 
many feeble magnetostrictive jumps are recorded. 
(3) The largest jumps of both kind occur on the 
steep part of the loop. Both become smaller as 
saturation is approached. 

These results may be interpreted as follows. 
On the steep part of the loop the magnetization 
vectors swing suddenly through large angles, 
complete reversals of 180° being frequent, but 
other changes of smaller angles, say 90°, also 
occur. As saturation is approached complete 
reversals no longer occur,—in fact, the angular 
jump becomes so small that the Barkhausen 
effect is very weak, although appreciable mag- 
netostrictive jumps occur. 

Another possible view is that as saturation is 


approached the sizes of the domains affected 
become smaller and smaller. However, since the 
Barkhausen effect becomes, relatively, much 
smaller than the magnetostrictive effect, it is 
necessary to suppose that some factor other than 
domain size is changing. It is natural to assume 
that reversals disappear and that the jump 
angle of the vector is smaller as saturation is 
approached. 

According to Kennard‘ a domain, in order to 
be in equilibrium when no external field acts, 
must be a long, slender filament which is mag- 
netized longitudinally. For such a domain we 
expect the magnetization vector to swing through 
180° when it jumps in a weak external field. The 
existence of jumps of smaller angles in weak 
fields is indicated, however, by the magneto- 
strictive impulses of the upper section of Fig. 2. 
At present the explanation of this discrepancy is 
not apparent. 


*E. H. Kennard, Phys. Rev. 55, 312 (1939). 
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The velocity of longitudinal waves in cylindrical bars may be expressed as the velocity at 
infinite wave-length times a function of two variables: Poisson's ratio, and the ratio of the 
diameter of the bar to the wave-length. This function is computed over the domain of the 
arguments which is of physical interest. Asymptotic values for the velocities at very short 
wave-lengths are deduced, and the variation of the displacement as a function of the radius is 
discussed. It is found that a similar analysis can be applied to torsional and flexural waves. 


N the light of recent experiments dealing with 
the propagation of ultrasonic elastic dis- 
turbances in cylindrical bars,' it seems desirable 
to carry out a more complete numerical develop- 
ment of the theory than is at present available. 
Existing treatments of the subject are open to 
the common criticism that, as they approach 
rigor, they become too unwieldy to be of prac- 
tical value, and the numerical computations 
based upon them become prohibitively long. It 
is the purpose of the present article to provide a 
concise and accurate method of analyzing the 
results of experiments dealing with longitudinal 
waves in slender cylindrical bars. Certain facts 
pertaining to the propagation of other types of 
disturbances will be touched upon. The solutions 
upon which the present work is based are the 
familiar ones due to L. Pochhammer,? presented 
in detail by Love,’ and discussed by several other 
authors." ® 
We deal first with the velocity of propagation 
of longitudinal waves in an infinitely long rod. 
The frequency equation for this case is suggested 
by Love.* This equation asserts a functional 
relation between wave-length, frequency, radius 
of the bar, density, and the two elastic moduli 
of the isotropic material under consideration. We 


*Paper No. 71, published under the auspices of the 
Committee on Research in Experimental Geology and 
Geophysics and the Division of Geological Sciences at 
Harvard University. 

t Now at the David W. Taylor Model Basin, Navy 
Department, Washington, D. C. 

1S. K. Shear and A. B. Focke, Phys. Rev. 57, 532 (1940). 

2?L. Pochhammer, J. f. Math. (Crelle) 81, 324 (1876). 

3A. E. H. Love, Mathematical Theory of Elasticity 
(Cambridge University Press, 1927), fourth edition, p. 287. 

4G. 5S. Field, Can. J. Research 11, 254 (1934), etc. 

5 R. Ruedy, Can. J. Research 5, 149 (1931). 

® Reference 3, p. 289. 


thus appear to be dealing with an equation 
involving six variables, which we would like to 
be able to solve explicitly for any particular 
variable as a function of the other five. The situ- 
ation is further complicated by the transcendental 
nature of the functions involved. It is important 
to notice, however, that, by a suitable choice of 
parameters, the equation may be reduced to a 
function of three variables, each of which is to 
be determined as a function of the other two. 
Such a functional relationship is readily ex- 
pressed in tabular form by computing a com- 
paratively small number of roots. The simpli- 
fication thus introduced is doubtless to be 
anticipated from the standpoint of dimensional 
analysis. 

We define the following quantities, and tab- 
ulate their equivalents in terms of Love's 
notation: 


Lovr’s 
SYMBOL DESCRIPTION NOTATION 
v Velocity of the wave b/y 
Vo “Bar velocity” 
B (1—20)/(1—¢); ¢=Poisson’s ratio 
x (v/v0)?(1+¢) 
L Wave-length 2x/y 
d Diameter of the bar 2a 
h y(Bx—1)! h’ 
k y(2x—1)! k’ 


e(y) 


With these exceptions we adhere to Love's 
usage. Clearly x is intrinsically positive. We shall 
also assume that 0<8<1, a restriction which 
implies 0<o0<}. If materials should be found 
with o<0, our analysis would have to be 
re-examined, and important modifications would 
presumably be forthcoming. 
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The frequency equation, derived from Love’s Eq. (54), p. 289, takes the form 


O*Jo(ha) p*pr 
da? 
dJo(ha) 


da 


m 


By means of suitable substitutions, and a 
somewhat space-consuming manipulation, (1) 
reduces to 


(x—1)* (ha) —(8x—1)[x—y(ka)]=0. (2) 


Referring to the definitions, we see that the 
latter equation is of the form F(x, 8, ya)=0, 
and that therefore the transformation to a 
function of three variables has been accom- 
plished. 

In discussing the characteristics of Eq. (2), 
the following properties of the function, ¢(y), 
are important: 


¢(0)=2; 
lim g(iy) =y. 


Zeros and poles of ¢ correspond to zeros of Jo 
and Ji, respectively. 

Let us now suppose Eq. (2) solved explicitly 
for x, so that 


x=x(8, ya). (3) 


The surface described by Eq. (3) consists of 
sheets over the (8, ya) plane. These sheets will 
be designated as the zeroth, first, etc., according 
to their heights above the plane. The case of the 
zeroth sheet is trivial, for Eq. (2) is satisfied 
identically by x=0, independent of 8 and ya. 
The first sheet is the one corresponding to the 
ordinary propagation of longitudinal waves, and 
is doubtless the most important, but the higher 
sheets may be of some interest. 

The first sheet is completely described by 
Table I, over a reasonably wide range of the 
independent variables 8 and ya and the results 
have been plotted in Fig. 1. The parameters 
of Eq. (2) are not particularly convenient except 
for purposes of computation; accordingly more 
conventional variables have been used in the 
preparation of the final table. The dependent 
variable is v/v», while the independent variables 


Jo(ha 


dJ (ka) 
da 


are o and d/L. Numerical uncertainty should not 
exceed 2 units in the last place, and interpolation 
in the table should be accurate to four significant 
figures up to d/L=1. 

Of particular interest are the values of x cor- 
responding to ya=0 and ya= ~. When ya=0, 
Eq. (2) yields-at once 


x(B, 0) 


but some manipulation is required to obtain 
x(8, ©). We must suppose that for some large 
value of ya, the arguments of the ¢ functions 
become and remain pure imaginaries as we move 
out along the sheet, and also that the arguments, 
ha and ka, increase without limit. Under these 
conditions we have g(ha)+—iha and g(ka)= 
—ika, and Eq. (2) reduces to 


x — 4x°+ (6—28)x+(8—2)=0. (4) 


If we make suitable substitutions to eliminate x 
and to express this result in terms of Poisson's 


2uy 
=0. (1) 


TABLE I. v/v as a function of d/L and oe. 


a 
o 
w 
nN 
o 
N 
N 
So 
nN 
n 


1.00 | 0.70196 0.70058 0.69768 0.69344 0.68814 0.68203 0.67537 
1.20 | 0.65419 0.65266 0.65030 0.64712 0.64321 0.63869 0.63368 
1.40 | 0.62836 0.62623 0.62361 0.62048 0.61687 0.61284 0.60844 
1.60 | 0.61393 0.61118 0.60815 0.60479 0.60111 0.59713 0.59289 
1.80 | 0.60565 0.60236 0.59892 0.59526 0.59139 0.58731 0.58304 
2.00 | 0.60080 0.59707 0.59326 0.58932 0.58524 0.58101 0.57664 


~ 
+ 
~ 
~ 
n 
~ 
= 


% | 0.60213 0.59491 0.58804 0.58148 0.57516 0.56903 0.56307 


| 
q 
1 
| 
| 

0.10 0.15 0.20 0.25 0.30 0.35 0.40 

1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 

0.99994 0.99986 0.99975 0.99961 0.99944 0.99924 0.99901 
| 
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Fic. 1. v/t as a function of d/L for various values of o. 


ratio and the velocity, we find that this cubic for 
x expresses a result identical with that deduced 
by Rayleigh for the velocity of surface waves 
(Rayleigh waves) upon a semi-infinite isotropic 
solid. The essential properties of Rayleigh’s 
equation are well known; for the possible values 
of o (or 8) there is always exactly one real root, 
which is positive. Computation of the roots of 
(4) for various values of 8 reveals at once that 
the initial suppositions made in deriving it are 
both satisfied. The results of the computation 
appear in Table I, opposite the entry d/L= ~. 

It is interesting to note in Fig. 1 that the curves 
for v/v9 drop below the asymptotic values de- 
duced from Eq. (4) and approach the asymptotes 
from beneath. We were led at one point to doubt 
that the curves actually approached the proper 
values, but computations carried out for d/ZL = 20 
indicated clearly that beyond a certain point, 
the values of v/vp increased as predicted. The 
possibility of an oscillatory behavior for v/v» at 
large values of d/L is precluded by the imaginary 
character of the arguments of the ¢ functions in 
Eq. (2), and by the monotonic nature of the 
¢ function of an imaginary argument. The 
minimum value thus implied for the velocity of 
longitudinal waves is a somewhat unexpected 
result, and one which would perhaps be difficult 
to verify experimentally. We have not computed 
the values of d/L for which the minimum occurs, 


but a rough estimate indicates that for most 
substances, the wave-length would have to be 
about } the diameter. The determination of the 
position of the minimum by the methods of the 
calculus is rendered laborious by the tran- 
scendental nature of Eq. (2). 

Further examination of the curves of Fig. 1 
reveals the possibility of determining Poisson’s 
ratio from measurements of the dispersion of 
longitudinal waves. If we plot v/v» as a function 
of o for various values of d/L, it will be observed 
that o is best determined by measurements of the 
velocity when d/L =}, or thereabouts, a condition 
which may readily be obtained experimentally. 

The statement is occasionally made’ that at 
very high frequencies, the velocity of propagation 
of longitudinal waves in a bar should approach 
the velocity of compressional waves in an infinite 
medium. In the case of the infinite medium, the 
wave involves no motion in the plane of the 
wave front, and the displacement is uniform. In 
the case of a long bar, motion in the plane of the 
wave front is inevitable, and the displacement 
is far from uniform, as will appear below. It is 
hard to see how a valid analogy can be drawn 
between the two cases. 

Whether or not the Rayleigh velocity can be 
observed experimentally at very short wave- 


7G. S. Field, Can. J. Research 5, 624 (1931). 
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Fic. 2. Axial component of the displacement as a function 
of r/a for various values of d/L, ¢=}. 


lengths depends on the stability of the mode of 
vibration with which it is associated. It is pos- 
sible that. beyond some limiting value of the 
frequency, velocities corresponding to one of the 
higher sheets of Eq. (3) will prevail. At the 
present time these sheets appear to be of purely 
academic interest. It may be shown that they 
all approach asymptotically the plane x= } when 
ya becomes infinite, and that they do so from 
above. For extremely small values of ya, the 
values of x become infinite. Thus each of the 
higher sheets contains a curve corresponding to 
the velocity of compressional waves in an 
infinite medium (as a function of 8). For any 
particular value of 8, this curve determines a 
discrete value of ya, a circumstance which has 
been pointed out by Field.’ The solutions for the 
displacements corresponding to velocities given 
by the higher sheets always involve one or more 
nodes as a function of the radius; the problem 
of verifying the existence of such nodes of 
vibration appears somewhat difficult. 

It seems worth while to give brief attention to 
the displacements corresponding to the first 
sheet. The axial: component, u,, is described by 


the equation 


where W is a function of the radius alone; viz., 
Ci a 

W=AtyJ (hr) +— {rJ,(kr)}. (5) 
r or 


The ratio A/C in (5) was determined in the 
course of obtaining Eq. (1). Utilizing this fact, 
and converting into the notation which we have 
introduced for the purposes of computation, we 
find 


1—x Jo(hr) Jo(kr) 
+(2x—1)*———-},_ (6) 
(8x Ji(ha) Ji(ka) 


where B is arbitrary, and is proportional to the 
amplitude of vibration. It will be noted that the 
signs in (6) are independent of the choice of sign 
in the extraction of the square root, for the 
radicals occur only in conjunction with J;,'s 
having the same radicals for their arguments. 
Similar considerations exclude the possibility of 
a complex solution. The latter point is of interest 
in that it implies that the displacement is in 
phase at all points of the bar. 

We have plotted W as a function of r/a for 
various values of the wave-length in Fig. 2, 
using values of the velocity taken from Table I, 
with o=}. Each curve for W(r/a) has been 
adjusted to make the maximum amplitude +1. 
When d/L=0, we have of course a uniform dis- 
placement. But as the wave-length is decreased, 
the displacement at the center of the bar rapidly 
decreases, and vanishes when d/L is about }. 
Then for a small range of wave-length the dis- 
placements at the center of the bar are of op- 
posite sign to those on the surface, and at some 
point along the radius, we encounter a node. 
This situation is altered as d/L becomes slightly 
greater than 0.4, when the displacement at the 
surface of the bar vanishes. As the wave-length 
is still further decreased, the displacement reap- 
pears at the surface, and when d/L=1, it is 
nearly uniform once more. After this, the dis- 
placement at the center of the bar gradually 
decreases, until finally, when the wave-length 
becomes infinitesimal, the motion is confined 
strictly to the outside surface. The Rayleigh 


| 
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wave velocity obtained for the limiting case of 
Eq. (2) corresponds to these surface waves. 

The reduction of the frequency equation to a 
function of three variables can be carried out as 
above for the cases of both torsional and flexural 
waves. In both cases, the parameters and func- 
tions required are the same as for the longitudinal 
waves. 

For a torsional disturbance, we find at once 


o(ka) (7) 


which yields the velocities not only for the 


2x—1 
g(ha) —2 g(ka) —2 
g(ha)—1 —(x—1)[¢(ka)—1] 


A partial check on the accuracy of (9) was ob- 
tained by substituting the first few terms of the 
series expansion for the ¢ functions, to obtain 
the simple formula for the approximate velocity 
of flexural waves. Construction of a table of 
roots of (9) similar to Table I would involve a 
considerable expenditure of time. If limiting 
values of the ¢ functions of an imaginary argu- 
ment are substituted in (9), we again obtain 
Rayleigh’s equation. The form of the displace- 
ments for this limiting case has not been at- 
tempted, but it seems likely that here also the 
displacement becomes confined to the surface. 

In applying the solutions of Eqs. (2), (7), and 
(9) to bars of finite length, only torsional waves 
yield a solution which permits rigorous satis- 
faction of the boundary conditions on plane 
surfaces at the ends of the bar. For the torsional 
case, then, we may expect that the solution 
applies regardless of the length of the bar, even 
when it is so short as to become a disk. For the 
other two cases, it is possible to determine the 
wave-length so as to make the normal traction 
across the end surfaces vanish, though a small 
shearing traction remains. Experimentally, we 
know that this shearing traction is zero when 
resonance conditions obtain. If we consider a 
slender bar, vibrating so that it contains many 
wave-lengths, we expect Pochhammer’s solutions 
to be exact, except in the neighborhood of the 
ends. Here we may suppose that the wave-length 


BANCROFT 


familiar case where the displacement is propor- 
tional to the radius, but for the more complicated 
cases as well. For these cases we find 


(8) 
o+1 


where v is again (E/p)', and C is a constant ob- 
tained by solving Eq. (7). 

For the case of flexural waves, even the sim- 
plified frequency equation is probably too com- 
plicated to be of much use, but it seems worth 
while recording it ; viz., 


2x—1 
—2[ ¢(ka) —2]—(ya)?(2x—1)| =0. (9) 
1 


becomes slightly altered, and if we attempt to 
compute it by dividing the length of the bar by 
the number of waves it contains at resonance, a 
small error will be introduced. This discrepancy 
has led to considerable confusion as to the 
precision which may be expected of Pochham- 
mer’s theory in any specific case. It is certain 
that the failure makes it impossible to develop 
a solution for a thin disk from the solution for a 
cylindrical bar, for in the case of the thin disk, 
the boundary conditions on the flat surfaces are 
of paramount importance, while those on the 
curved surface are trivial—quite the reverse of 
the situation for the bar. On the other hand, 
experimental evidence indicates that when the 
ratio of diameter to length is less than 0.4, the 
error involved in neglecting these stresses is less 
than 0.05 percent. Table II gives velocities of 
longitudinal waves in a set of specimens cut from 
a single length of 3-inch drill rod. The measure- 


TABLE II. Observed velocity of longitudinal waves in steel 


bars of different lengths. 
LENGTH WAVE-LENGTH RESONANT VELOcITY 
(cm) (cm) FREQUENCY (KM/SEC.) 
15.235 10.157 50473 5.126 
5.083 10.166 50426 5.126 
15.235 7.518 67184 5.118 
3.812 7.624 67162 5.120 
15.235 5.078 100390 5.098 
2.540 5.081 100364 5.099 
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ments were made by the method of Bancroft and 
Jacobs.* It will be seen that, within the experi- 
mental error, the observed velocity as computed 
from the assumed wave-length is independent of 
the length of the specimen, and is a function of 
the wave-length alone. 

Except for the effects of the above-mentioned 
discrepancy, which we have shown experi- 
mentally to be slight in the case of longitudinal 
waves, it is felt that the velocities determined 
from Pochhammer’s solutions may be used with 
perfect confidence. 

In conclusion, it appears that the work of 
Shear and Focke! is explained on a quantitative 
basis. Table I reproduces the behavior of the 


_ §D. Bancroft and R. B. Jacobs, Rev. Sci. Inst. 9, 279 
(1938). 


longitudinal vibrations which they observed 
with remarkable fidelity. The flexural vibra- 
tions behave in much the way one would ex- 
pect in the light of Eq. (9), and it is note- 
worthy that the experimental data suggest 
that a common asymptote for the flexural and 
longitudinal velocities at high frequency is not 
improbable. The observed torsional vibrations 
may possibly belong to one of the sheets of Eq. 
(7). It also seems likely that the unexplained 
points at high frequency lie upon one of the 
higher sheets of Eqs. (2), (7), or (9), for in a 
qualitative way they lie in the region associated 
with these more complicated vibrations. 

It is a pleasure to acknowledge our indebted- 
ness to Professor Francis Birch for his encourage- 
ment and help, particularly in checking most of 
the rather tedious algebraic work. 
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The Coulomb repulsion between protons must give rise to a tendency for the proton density 
to vary within a nucleus from a minimum value at the center toa maximum near the boundary. 
A non-uniform proton density may be expected to create forces which distort the neutron 
distribution and tend to make the two distributions vary in the same manner. If surface 
effects are neglected, it is possible to calculate the energy correction associated with the 
non-uniform densities without making special assumptions about the nuclear forces. The 
neglect of surface effects permits the assumption that the variations in density are small 
departures from essentially constant distributions. It is found that the departure from uniform 
density is appreciable in heavy nuclei, but the energy correction is negligible. 


1. INTRODUCTION 


HE Coulomb repulsion between protons 
must give rise to a tendency for the proton 
density to vary within a nucleus from a minimum 
value at the center to a maximum near the 
boundary. A non-uniform proton density may 
be expected to create forces which distort the 
neutron distribution and tend to make the two 
particle densities vary in the same manner. 
If surface effects are neglected, it is possible to 
calculate the energy correction associated with 
the non-uniform densities as well as the densities 


themselves in a comparatively rigorous and 
simple manner. The neglect of surface effects 
permits the assumption that the variations in 
density are small departures from essentially 
constant distributions. A suitable model for the 
systematic neglect of surface effects is provided 
by supposing the nuclear system enclosed in a 
box of radius R. At the boundary the radial 
derivative of the wave function with respect to 
any radial coordinate must vanish: 


0 
5p =0, =R. (1) 


| 
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This boundary condition makes the particle 
densities have zero slope at r=R in contrast 
with the correct boundary condition which makes 
the particle densities vanish at r= «. When the 
Coulomb interaction is omitted, the boundary 
condition implies essentially constant particle 
densities throughout the nucleus. We postulate 
the existence of a volume energy density, 
&(p», px), Which depends only on the particle 
densities. The total energy is 


(2) 


with 
E,= f f f wate. (3) 


f f f (4) 


(5) 


The first term in Eq. (2) is the volume energy 
which includes the major part of the kinetic 
energy and all the potential energy derived from 
the specifically nuclear forces. No special assump- 
tion will be made in this paper about the nuclear 
forces or the energy density 6 An exchange 


‘term, calculated by the statistical method, is 


included in the Coulomb energy.! The third 
term, E», is Weizsicker’s formula which gives a 
lower limit for the additional kinetic energy 
associated with non-uniform particle densities.! 

To determine the nuclear radius, use is made 
of the auxiliary quantity 


R,=A'R((N—Z)/A) (6) 


obtained by solving the equation 
—E,(RNZ)=0. (7) 
oR 
From the minimum condition, 
—E(RNZ)=0, 8 
(RNZ) (8) 


1C. F. v. Weizsacker, Zeits. f. Physik 96, 431 (1935). 
Equation (50) of this reference lacks a numerical factor 
(3/4)! in the right-hand member of the equation. 


we obtain 
R=R,(1+E,./E,"), (9) 
(10) 
E,. is the Coulomb energy for the radius R, and 
(11) 
In this discussion we take? 
R,=0.5A*(e?/mc*), E,’’=100Amc*. (12) 


In Eqs. (6)-(12), the particle densities are 
assumed to be constant. After the energy cor- 
rection for non-uniform densities has been calcu- 
lated, it will be clear that the associated correc- 
tion to the radius is negligible. 

The perturbation method developed in Section 
2 requires knowledge of the second partial 
derivatives of E,(RNZ) with respect to N and Z. 
Section 3 contains the derivation of relations 
expressing the partial derivatives in terms of 
total derivatives which can be evaluated from 
empirical or semi-empirical data. Numerical re- 
sults are obtained and discussed in Section 4. 


2. THE PERTUBATION METHOD 


The equations 
4rR*p,o=3N, (13) 
(14) 
express the particle densities in terms of functions 


v, and v,. Consistency requires that v, and 2, 
satisfy the orthogonality conditions 


1 1 
ff f vpx*dx =0, (15) 
0 0 


while the additional condition 
dv, (x) /dx = dv,(x)/dx=0, x=0,1, (16) 


follows from Eq. (1) and the continuity proper- 
ties of the wave function. 
It is now possible to express the energy as a 


2 E. Feenberg, Phys. Rev. 59, 149 (1941). 
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power series in v, and v, with the result 


E, Pro) 


3 
+1 ( f ) 
2 0 aN? 
1 
+( f 
0 


00 
+2/ f ) 
0 ON dZ (17) 


E,=3e?Z(Z—1)/5R 


1 
f 
0 


1 z 
+15 f f 
0 0 
1 
(18) 
0 
: 1 
Ey, = 3h?, f (dv,/0x)*x*dx 
0 


+N f | (19) 


All terms in v, and v, of order higher than the 
second have been dropped in keeping with the 
assumption that the particle densities do not 
depart greatly from their average values. 

Now we introduce a set of polynomials which 
satisfy the conditions (15) and (16) and are 
normalized and orthogonal on the interval 
0x1 with the weight factor x*. The functions 
v, and v, are written as linear combinations of 
these polynomials with coefficients determined 
by the condition that the energy be made as 
small as possible. Actual calculations have been 
carried through with only two members of the 
set: 


= 6.1675 {1—3.75x°+2.5x'}, 
v2= 11.028 { 1 — 13.7879x? (20) 
+24.1819x*—11.2424x*}. 


It should be clear from the following table that 
any reasonable variable particle density can be 
represented rather well by a constant plus a 


suitable linear combination of v; and v2. 


x 0 0.2 0.4 0.6 0.8 1.0 


6.167 5.366 3.454 1.172 -—0.740 —1.542 
11.028 6.881 0.603 -—2.121 —0.349 1.672 


A number of integrals involving v; and v2 are 
required. These are 


1 
f = —0.1432, 


1 
f vex'dx = 0.0428, 
0 
1 z 
f Vix f v1(x")x’*dx'dx = 0.02470, 
0 0 
1 z 
f f vo(x")x"*dx’dx = 0.00827, 
0 0 
1 z 
f f v2(x")x"*dx'dx = 0.05400, (21) 
0 0 
1 z 
f Vox f v(x") x’*dx'dx = —0.05201, 
0 0 
1 
f (dv; dx) *x*dx = 20.378, 
0 
1 
f (dv2/0x)*x*dx = 61.302, 
0 


1 
f (dv; / 0x) (Ov2/0x)x*dx = — 2.014. 
0 


We write 
= 4,01: +-),v2, 
(22) 
Ue 
and obtain 


E, = Pro) 


3 2 
— 
2 oN? 


+ (a,2+b,2)Z2— 


+2(a,a,+b,b; RNZ), (23 
(a,a,+b,b;) ), (23) 


wow 
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TaBLe I. Partial derivatives of E,.(R.NZ)/me*. 
24 50 1.070 1.384 —0.773 
44 100 0.481 0.760 —0.379 
62 150 0.295 0.543 —0.248 
80 200 0.212 0.421 —0.183 


+0.3580a,—0.1070b, +0.3705a,? 
+0.1241b,?+0.0299a,b, 

(24) 

Ey =156.2A(R/R,)? 

x [(2N/A)(a,2+3.016,? —0.200,b,) 
(25) 


3. PARTIAL DERIVATIVES OF £, WITH 
RESPECT TO N AND Z 


The only assumption required in this section is 
that E, be given as a function of R, N and Z. 
There is no need for the special assumption (3) 
and the surface.energy need not be neglected 
since it does not appear explicitly in the calcula- 


tion. To avoid any implication of dependence on _ 


the assumption (3) new symbols Eo, Ro are 
introduced here in place of E,, R,. Equation (8) 
is required in the form 


—E,(RNZ) =E,/R. (26) 
aR 
We use a number of relations of the type 
0 
— —E,(RNZ) 
aN OR 
d 
— 
dN dN dR/AR 
(27) 
dN 
dfd dRa 
= ——— — JE,(RNZ). 
dN dN 


These relations when combined with Eq. (26) 


vield 


_ “Cos ~) 
oN? dN dR? 
E. dR dR d E. 


dN\RdN/) dNdNR’ 
d 


az? aR? 
(28) 
E.dR\ dRd E. 
dZ\RdZ/ dZdZR’ 
ae dad aR dR 0°E, 
aNaz. dNdZ dN dZ aR? 


dRd dRd\E, d 


— —+— — 
dNdZ dZdN/R RdNdZ 
To obtain the corresponding equations for 
Eo(RoNZ) we replace R by Ro in Eq. (28) and 
omit the terms containing the Coulomb energy. 
There results 
d*Ey Ey” dRo 3 


ON? R,? 


(29) 


d*Eo 
0Z? 
dd Eo” dRy dRo 


dNdZ RdN dZ 


If is replaced by Eo’’/Re? 
and only terms containing E»”’ to the first and 
zeroth powers are retained, the right-hand 
members of the Eqs. (28) become identical with 
the corresponding right-hand members in the 
Eqs. (29). This means that the second partial 
derivatives of E,(RNZ) in Eq. (23) may be 
replaced by those of E,(R,NZ) without appre- 
ciable error. 


TABLE II. Energy correction and density coefficients. 


ZA a, by a, b, px (R)/p, (0) p,(R)/p, (0) 
24 50 —0.0215 0.0031 —0.0092 0.0006 —0.7 mc? 1.15 1.07 
44 100 —0.0368 0.0063 —0.0167 0.0016 —3.4 mc? 1.27 1.12 
62 150 —0.0481 0.0092 —0.0222 0.0026 —7.8 mc? 1.35 1.16 
80 200 —0.0580 0.0119 —0.0269 0.0037 —14.3 mc* 1.43 1.20 
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To evaluate the total derivatives we use the 
energy formula 


E,(R.NZ)/mc?= —aA 
+[96| (30) 


which contains Wigner’s linear term in | N—Z}. 
The quadratic term comes partly from the 
kinetic energy and partly from the potential 
energy. Wigner’s more general formula* reduces 
to Eq. (30) for even values of N and Z. The 
major contributions to the second derivatives 
come from the quadratic term 25(N—Z)?/A. If 
the numerical coefficient 25 is too small, which 
seems likely, our results for the energy correction 
and the deviations from uniformity will come 
out too large. Another approximation tending in 
the same direction is that involved in the use of 
Weizsicker’s E, for the extra kinetic energy 
term. Thus the energy correction and the devia- 
tions from uniform density given by this calcu- 
lation must be interpreted as upper limits on 
the true values. Table I contains numerical 
values of the partial derivatives computed from 
Eqs. (29) and (30). 


4. NUMERICAL RESULTS 


The energy correction produced by the quad- 
ratic and linear terms in a,, a,, },, 6, can be 
written as 


= — + + pa,?+qb,* 
— — (31) 
with coefficients in mc? units defined as follows: 
Ca=0.430Z(Z —1) / A}, 
=0.128Z(Z—1)/A}, 


#E, 
156A1(22/A) +0.4412(Z —1)/A! 


—0.611Z'(Z—1) A}, 


3 E. Wigner, Phys. Rev. 51, 947 (1937). 


“+470A'(2Z/A)+0.149Z(Z—1)/A! 
—0.611Z"Z—1)/A!, (32) 


2 


3. @PE, 
py 
2 dN? 


3. @E, 
=-N*——+470A*(2N/A), 
2 dN? 


t,=15.4A4(2N/A), 
tr =15.4A1(2Z/A) —0.018Z(Z—1)/A!. 


Equation (31) can be reduced to diagonal form 
by a repeated application of the process of com- 
pleting squares. The numerical results of mini- 
mizing 6E are listed in Table II. 

Certain features of the results in Table II can 
be understood without calculation by examina- 
tion of Eq. (31), remembering that t, and ¢, are 
almost negligible. It is evident that a, must be 
negative and b, positive. Also a,a, and 6,b, must 
both be positive which makes a, negative and 6, 
positive. The negative value of #@E,/dNdZ 
requires the neutron and proton particle densities 
to vary in a parallel manner resulting in a 
tendency toward the formation of a _ hollow 
center within the nucleus. The dependence of 5E 
on N, Z is very similar to that of the Coulomb 
“expansion” energy E,.2/2E,", but the latter 
correction is three to four times the size of the 
former. Finally it must be stressed that approxi- 
mations in this calculation make the computed 
magnitude of 5E somewhat greater than the true 
value. Also the coefficients a,, 6,, a,, 5, are all 
somewhat too large in absolute value. Satis- 
factory convergence is indicated by the small 
values of the },/a,, ratios. 

E. Wigner has estimated 6E by another method 
with results similar to those stated in Table II. 
The writer wishes to thank Professor Wigner for 
information about his calculation in advance of 
publication. 
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The electrostatic energy of electrons in a lattice has been calculated by the Rayleigh- 
Schrédinger perturbation method. This method gives a value which becomes logarithmically 
infinite for a metal. For an insulator, however, the gap in energy above the first occupied 
zone leads to a finite result. Numerical values have been found for the three cases in which 
the gap is one-half the width of the first zone, equal in width, and twice as wide. 


1. INTRODUCTION 


HE customary approximation for the solu- 

tion of problems of atomic or crystal 

physics is the Hartree-Fock approximation. The 

question of deviations from this approximation 

and of an improvement of the underlying picture 
deserves therefore considerable attention. 

We shall be concerned with the case of doubly 
occupied orbits, i.e., with such states in which 
each orbit is either empty, or occupied by two 
electrons with opposite spins. It is well known 
that in this case, for the calculation of the 
electrostatic and kinetic energy, one can use a 
wave function of the form 


| 


The number of electrons is 2, the symbols 
X1, ***, Xny V1, Yn Stand for three-space coordi- 
nates each. The x denote coordinates of electrons 
with one-spin orientation, the y coordinates of 
electrons with opposite spin orientation. The 
¥1, are orthogonal and normalized wave 
functions, representing the m occupied orbits. 
The wave function (1) represents, of course, a 
singlet state of the system. 

Two methods have been tried so far to obtain 
an improvement of the wave function (1). The 
first one is essentially based on the variational 
idea: The sum of kinetic and potential energies 
of more general wave functions than given by 
(1) is calculated and the functions and parame- 


*On sabbatical leave from Allegheny College, Mead- 
ville, Pennsylvania. 


ters which occur in this wave function are 
determined by the requirement that the total 
energy be a minimum. The difficulty encountered 
is that it is, in the case of many electrons, not 
easy to find a wave function for which the total 
energy can be calculated. For a small number of 
electrons, on the other hand, the variational 
method is relatively easily applied and converges 
in many cases rapidly.' Even in the case of many 
electrons, in particular crystals, some success 
has been achieved by this method.? 

The second method which will be tested in the 
present paper is the Rayleigh-Schrédinger per- 
turbation method. The application of this 
method encounters relatively small difficulties 
in the case of only doubly occupied orbits, 
because the original state is not degenerate and 
the successive approximations are easily obtained 
in the form of infinite series.’ It is generally 
assumed that the successive approximations 
converge in the case of a small number of 
particles, that is when m is small. In the case of 
n large, with which we will be concerned, the 
question of convergence is more complicated 
and we first turn to a discussion of this question. 

If the state considered is nondegenerate, 
the first approximation for the energy in the 
Rayleigh-Schrédinger scheme is given by the 
expression 


EO =E,= f (2) 


1 Cf. the well-known papers of Hylleraas and James and 


Rev. 46, 1002 (1934); Trans. Faraday 

Soc. 34, 678 (1938). 

os. Brillouin, J. de phys. et rad. [7] 3, 565 (1932); 
S. Flugge, Zeits. f. Physik 96, 459 (1935); S. Watanabe, 
Zeits. f. Physik 113, 482 (1939). 

The conclusions reached in the present paper are, however, 
in many respects different from those of these authors. 
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where H is the Hamiltonian. This is the energy 
in the Hartree-Fock approximation. In the 
second approximation, the wave function is 
given by 


ai, (3) 
i 


f Ey. (3a) 


Here VW, are wave functions of the form (1) with 
the difference, however, that the orbits occurring 
in the determinants are not the occupied orbits 
of Wo but at least partly other members of a 
complete orthonormal set which contains y¥:, 

-+, wa. In the “excited states’? V,, the two 
determinants are in general different because in 
these states the orbits are not all doubly oc- 
cupied. V,*HY¥, is the first approximation 
to the energy of the excited state /. The energy 
of the lowest state becomes in the second 
approximation 


(4) 


The summation over / must be extended in (3) 
and (4) over all possible excited states, i.e., 
over all wave functions in which the two de- 
terminants contain any orbits y. 

If the Hamiltonian contains, in addition to 
the kinetic energy, only terms each of which 
involves the coordinates of not more than two 
particles, i.e., two-particle interactions alone, a 
considerable simplification can be effected. In 
this case, a; vanishes if V, contains more than 
two orbits not contained in Wo. This can be 
seen directly upon writing down the integral 
occurring in (3a). As a consequence hereof, the 
wave function V»® contains only those states ¥, 
in which not more than two particles are “‘ex- 
cited” to states different from yi, ---, wa. (We 
shall see later that only such states occur in 
which two excited electrons are present in (3).) 
This means, physically, that a measurement of 
the number of electrons present in the states 
v1, ***, Yn gives at least the value 2n—2, and 
possibly the value 2m. For m large this is, of 
course, an unreasonable result as one will expect 


that the number of excited electrons is roughly 
proportional to ”. This is also obtained when 
the variational method is used. In addition to 
this difficulty, in all cases in which calculations 
have been carried out so far, the sum of the 
squares of the a; increases, with increasing n, 
beyond all bounds so that the wave function (3) 
cannot be normalized. Not only does (3) fail to 
contain all excited states which it should, but it 
contains more of some than it should. It is 
evident, therefore, that in the case of n large, 
(3) is not an approximate wave function in the 
usual sense, and it could be used only if the 
statistical interpretation is used in a modified 
form.‘ The reason for this difficulty is that Vo 
is not a good first approximation. Doubtless, for 
the real stationary state, the probability is 
vanishingly small that all electrons shall be in 
the unexcited states yi, ---, wa. Hence the 
coefficient of Wo is vanishingly small in the 
correct wave function while it is 1 in (3). 

It will be seen that in spite of this, under 
certain conditions, the expression (4) for the 
energy gives a reasonable result even in the 
case of m very large. This hardly can be inter- 
preted as the expectation value of the energy 
for the wave function (3)—because this wave 
function cannot be normalized and would give 
an entirely different result if it could be. It can 
be considered, however, as the second-order 
term in a power series expansion of the total 
energy, in terms of the fictitious parameter with 
which the perturbation energy appears to be 
multiplied in the Rayleigh-Schrédinger scheme. 

It remains, of course, questionable whether 
this power series for the total energy converges 
for as large a value of the parameter as is 
assumed in the actual problem. It is quite 
conceivable that the radius of convergence of 
the Rayleigh-Schrédinger series decreases with 
increasing m so that, for m large, it is smaller 
than the value which this parameter actually has. 
It is also possible that the number of terms 
needed in order to obtain a reasonable approxi- 
mation increases with increasing m, while, in 
practice, one can proceed only to the second term. 

In order to test this method it has been applied 
to two problems: the electron interaction in 


‘This has been essentially proposed by Watanabe, 
reference 2. 


600 Cc. W. UFFORD 


metals and the electron interaction in insulators. 
In the first case, it was assumed that the wave 
functions occurring in (1) are plane waves, i.e., 
that there is no lattice field present. Although 
this corresponds to a fictitious uniform distribu- 
tion of the positive charges, it is in reality for 
many cases not a bad approximation to the true 
Hartree-Fock wave function. In this case, the 
series (4) gives a correction which increases, with 
increasing m, more rapidly than n itself. Hence, 
the energy correction per particle, obtained from 
(4), increases beyond all bounds, with increasing 
n—an evidently incorrect conclusion. We believe 
that this result is independent of the approxi- 
mations employed. In this case, the Rayleigh- 
Schrédinger theory cannot be applied. 

In the case of an insulator, the lowest £; is at 
a finite distance from E»—the distance of the 
gap between the first two Brillouin zones. In 
this case Ey— E; cannot become arbitrarily small. 
For the wave functions we shall use 
in this case the expressions given by Peierls’ for 
the case of weak lattice fields. Although in this 
way only a few wave functions are very different 
from plane waves, those matrix elements which 
were very large for free electrons are materially 
reduced. The whole sum (4) becomes propor- 


f J 


tional to m. This gives a constant energy cor- 
rection per particle which is, furthermore, of the 
right order of magnitude. This proves, of course, 
only that our result is reasonable, not that it is 
correct. In fact, it cannot be correct if the 
condition for Peierls’ approximation is satisfied, 
i.e., if the lattice potential is very nearly con- 
stant. For a constant potential, we have again 
the case of free electrons in which (4) breaks 
down. In fact, we obtain an energy correction 
per particle which increases beyond limits if the 
lattice field becomes vanishingly weak, a result 
which cannot be correct. Our result can be 
correct, however, for moderately strong lattice 
fields. 

In justice to the Rayleigh-Schrédinger method, 
it should be pointed out that (4) always gives 
finite results in case of short range forces as they 
occur in nuclei.6 Whether this finite result is 
correct remains, in my opinion, an open question. 


2. THE ELECTRON INTERACTION IN METALS 


First the matrix elements a; of (3a) must be 
calculated for H=e?/r;; to find the electrostatic 
energy of interaction of the electrons. Substi- 
tuting the wave functions (1) in (3a) gives for 
the integral 


- |dxdy, (5) 


where dx =dx,---dx, and dy=dy,:-+-dy,. For plane waves 
¥x(xi) =exp [tk-x,] and E,=h?k?/2m. (6) 


The components of k have the form 2rn,/L, 2rn,/L, and 2rn./L where nz, n, and n, are integers’ 


and L is the length of an edge of the crystal, L}= V. The functions y, are normalized and orthogona 


so that 


(1/V) f = Bau 7) 


As will be shown presently, the matrix elements for which only one electron is excited vanish. 
Therefore only those with two electrons excited will be considered. There will be two types of matrix 


5R. Peierls, Ann. d. Physik 4, 121 (1930). 
®S. Watanabe, S. Flugge, reference 3. 
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elements: those for which two electrons with parallel spin are excited, and those for which two 
electrons with antiparallel spin are excited. As the greater contribution to (4) arises from the second 
type, it will be sufficient to consider it alone. 

Let k and k’ be two wave vectors in the reciprocal lattice which represent unexcited or occupied 
states, and 1 and I’ wave vectors representing excited states. In the sum over kk’ll’ to obtain the 
energy of all the electrons in (4), k and k’ must be vectors to all points inside the sphere of momentum 
space whose radius represents the highest occupied state, and 1 and I’ must be over points outside 
this sphere. The matrix elements (5) may then be written 


e 


Vix’ = 


since the same terms are obtained from the diagonal of the complex conjugate matrices as are 
obtained from each of the other terms. As there are n! terms in each matrix, the (m!)? cancels out. 
Again 
e 1 


| — Ver 


as the integration over all coordinates except x, and y,- gives unity because of the normalization 
of the y,(x;)s. Here k and k’ represent electrons with opposite spins. The parts of the matrix elements 
with parallel spin, where both electrons have the coordinates x or both y, are neglected. Because 
of exchange, these parts contain the difference of two terms which are nearly equal at those points 
where these parts have their largest value. Thus these parts are small compared to the one calculated 
in (9). 

Now if plane waves are introduced and a transformation 


Xi — Yu =20, Xe tyr =e (10) 
of Jacobian 8, then 
4e* dod 
View f exp +1’) - }— (11) 
te* do 
-—f exp [i(k—1+k' +1’) (12) 
i 


because of the orthonormality of the ¢ waves. 

If only one electron is excited, then l’=k’ and the factor 6(k—1) gives a factor zero, since if one 
electron is excited k#1. Thus all matrix elements with a single electron excited are zero. 

The integral over @ gives 


so that 
167e? 6(k +k’ —1—1’) 1 

=— 6(k+k’—1-1’). (14) 
V ji-—k/? 


Vixw = — 


The energy may now be calculated by substituting (6) and (14) in (4). The second-order term is 
Vix [2 5(k+k’—1-—1’) 


(15) 
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Fic. 1. Wave vectors 
k and k’ in the xy 
plane for the normal 
states inside the first 
zone, and 1 and I’ for 
the excited states out- 
side the zone. It will be 
seen that these vectors 
satisfy the condition 
that 1—k=k’-I’. 
is the vector to the 
boundary of the first 
zone. 


This sum can be evaluated as a special case 
of the sum considered later for insulators. It is 
found to approach infinity logarithmically.’ This 
is due to the difference of the wave vectors of 
the normal and excited states in the denominator 
of (15). This denominator approaches zero as 
the wave vectors approach each other at the 
boundary of the occupied states. In insulators, 
on the other hand, there is a gap in energy at 
this boundary with the first energy band com- 
pletely filled. The gap will introduce a term in 
the denominator which will prevent the energy 
from becoming infinite. This question of insu- 
lators will be considered. 

In (15), the energy of the free electron, 
E,=h*k?/2m, has been substituted for E,. This 
is, strictly speaking, incorrect, because the energy 
of the system contains exchange terms in 
addition to the kinetic energy of the electrons. 
Hence, actually® 


E.= (k2 km?) /2m 
In | (15a) 


should be substituted for FE, in (15). (The 
energy is normalized in (15a) in such a way 
that it vanishes for k=k,,). The second term of 
(15a) is unimportant in insulators, since there 
the energy denominator will contain a finite 
term (the energy gap between the zones), 
compared with which the second term of (15a) 
is negligible. However, in a metal, the second 
term of (15a) is actually more important than 
the first one, because it predominates in the 
region k~k,, which gives the largest contribution 


71 am grateful to Dr. Conyers Herring for a discussion 
of this case. ; 

§P. A, M. Dirac, Proc. Camb. Phil. Soc. 26, 376 (1930). 
Also L. Brillouin, J. de phys. et rad. 3, 565 (1932) and 
J. Bardeen, Phys. Rev. 49, 653 (1936). 


to the summation. If one wants to be quite sure 
that the summation gives a divergent result, 
one must substitute (15a) for (15) and carry out 
the summation. This can be done, but the 
calculation is not given here in detail. If the 
summation is replaced in the usual way by 
integration, the integrand is an expression the 
leading term of which, in the neighborhood of 
x=0, is (xInx)~-'dx. This integral diverges, 
although the divergence is much weaker than 
that of x~'dx which is obtained if the second 
term® in (15a) is neglected. This means that, 
strictly speaking, the second approximation gives 
an energy correction which is, per electron, 
proportional to InInm rather than to Inn. 
Although this is not a very bad divergence, it 
shows, nevertheless, that the Schrédinger per- 
turbation theory becomes increasingly worse for 
metals as the number of particles increases. 


3. THE ELECTRON INTERACTION IN INSULATORS 


In order to introduce the gap in energy found 
in insulators into the matrix elements of electro- 
static interaction used to calculate the second- 
order term in the energy of the electrons, a 
method due to Peierls’ will be followed. Because 
of Bragg reflection of the lattice waves by the 
crystal, the first-order wave functions are de- 
generate. Therefore a linear combination of 
first-order wave functions must be used. These 
may be obtained by modifying the plane waves, 
with the period of the whole crystal, by waves 
with the period of the lattice constant. Here 
only two waves will be used which correspond 
to values of the wave vector in the two lowest 
zones. The normal state is represented by the 


Fic. 2. x’ is to be integrated 


Y * throughout the cross-hatched 
Kea volume. «x and x—d lie inside the 
d . sphere with its center at 0. 


koa «’+x—2 lies outside this sphere. 
x’ lies inside and 4 outside the 
sphere with its center at 0’. 


*T am very much indebted to Professor John Bardeen 
for pointing out to me that it is necessary to include the 
exchange term when the calculation is carried out for 
the metal. 
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wave vector k in the first zone combined with that for which the x component is k,—2k,. The 
y and z components are the same. Waves with k’ and y component k,’—2k,» are combined for 


the other occupied state. Similarly the x components of the two wave vectors of one excited’ 


electron differ by 2k, and the y components of the other excited electron differ by 2k,. The 
vectors for a simple cubic lattice are shown in Fig. 1 together with points which correspond to 
the wave vectors with which they are combined. It will be seen that these vectors satisfy the 
requirement (12) that k+k’=1+I!' or 1,—k, =k,'—l,’, l.—k.=k,'—l,’. The z 
components are all equal. 

Thus for the normal electrons 


(1/V*) {cos $8 exp [7k-11]+sin $8 exp ], 


vir = (1/1) {cos 36" exp [tk’- ]+sin 36’ exp [i(k’- r2—2kmy2) J} 
and for the excited electrons 
¥i=(1/V?) {cos exp }y exp J}, 
= (1/14) {cos exp [il’-r2]+sin exp 
The energy corresponding to yx is 
= (h?/2m)[k? cos? 38+ sin? $8)]—2 sin $8 cos (17) 
= (h?/2m)[k?+ 2kn(Rm —kz)(1 —cos ]—7 sin B, (18) 


where 


is the first Fourier coefficient of the periodic lattice potential, P’. Similar expressions hold for E,-, 
E, and Ey. 8 is now chosen to make E; a minimum 


2 


n cos B=—Rk»(k» —Rz) sin B, (19) 
m 
mn 
an p= (20) 
Rz) 
h*m(Rm — Rz) 
2 mn 
from which 
B= (Rm —kz)/my approximately. (22) 
Similarly 
2 mn 2 mn 2 mn 
B h*Rm(le—kz) h*k»(l,’ — ky’) 
m\ez— = “Rm 
2 2 2mn 2 2 2mn 


The wave functions may be substituted in (9) and the matrix elements Vix: calculated. Only 
those four, of the possible products, for which k,, cancels out in the exponent give elements different 
from zero. Identical integrals as given in (12) are obtained for each of these four products. With 
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the values of these integrals from (14), the matrix elements of (9) become 


Vix = — (cos $8 cos $y sin $8 sin }y+cos 38 sin 38’ cos $y sin }y’+sin 38 cos 38’ sin }y cos 
4 2 
sig 38 sin 38’ cos cos }y’) (25) 
5(k-+k’ —1—1’) (26) 
2 2 V 
*h*km? (le—kz)(ly— Ry 
( M — 1-1, (27) 


where, for the last line, (24) has been used. As the matrix elements are complicated, (4) is difficult 
to integrate. It may be seen that the maximum value of (/,—k-)(l,—k,y)/|l—k|? is }, so that the 
maximum value of the matrix element is 1e*h*k»,?/2Vm*n?. This is the value when /=8, i.e. on the 
boundary of the first Brillouin zone. For large values of /, the matrix element is given by (14). 
Hence an approximation which is valid at /=& and large / is 


4re?- 6(k+k’—1-1’) 


(28) 
Vo 


= 


(4) will now give the second-order term in the energy. The energy difference between the normal 
and excited electrons may be found from (18). 


Eo— (h?/2m)(k? + 2Rm(km— kz)(1 —cos 8) +k’? + (Rm — ky) (1 —cos —1° 
— 2Rm(m—12)(1—cos y) — 2km(km —ly)(1—cos y’) ] 
— B+sin B’—sin y—sin y’) (29) 
= — (h?/m)(1—k) -(1—k’) —4n, (30) 


since on the boundary of the first Brillouin zone where 1,=1,=k.z=ky=km, (23) shows that B=,’ 
=2/2, y=y'=-—7/2. This is a valid approximation since it is only for small values of Eo—£, 
that the terms in the energy sum of (4) are large, and Ey)— E, is smallest on the boundary. Therefore 


V%2( (1—k) (1—k’) +4mn/h2] 


This value of the energy reduces to (15) for metals if »=0, i.e., when there is no periodic lattice 
field. In order to carry out the summation of (31), the simple cube of Fig. 1 will be replaced by a 
sphere of radius k,,. The sum over the wave vectors is then replaced by integration over k and k’ 
inside, and over | and I’ outside this sphere. The change from sum to integral gives a factor V*/(27)°. 
Let 


(31) 


E=% 1 — 


= —167°e'm V/h?(27)*, a=8mn/hk,,”, (32) 
then 
dkdk’'dl 
f (33) 
- (1—k’) +akm?/2] 


The integration of (33) may be simplified by setting kn» equal to 1, and integrating with the unit 
sphere as the boundary of the occupied states. This may be done by setting 


x=k/km, and 2%=1/km, (34) 
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which gives 


dxdx'd2 
f (35) 
+a/2] 


There are three cases to be considered: 


(a) (b) 1=|A—«|=2, and (c) 2={a— x]. (36) 
The integration will be carried out for case (a) first with respect to x’. Let 
|A—x|=a, and =c. (37) 
Then 
and (38) 


* is a vector to points inside a sphere of radius 1 in Fig. 2. x—2 will also be a vector inside the same 
sphere since |2—«| =1. 2, however, must be a vector lying outside a sphere of radius 1. This second 
sphere may be drawn with x—d as a center. From the end of x—4, x’ will be a vector lying inside 
the second sphere. Finally since 4’ = (x’-++x—2) must be a vector lying outside the sphere containing 
«x, the volume of integration for x’ is outside the first sphere and inside the second, and is shown 
cross hatched in Fig. 2. Using for x’ elliptic coordinates x’ and ¢ with the volume element 2mx’dx'cdc/a 
gives for the part of (35) depending on *’ 

dr 1 «’+a cdc 

f f 

“14 


= (x /4a*) | [4a?— (b?+a)?] In (2a+b?+a) In 
+ In +a+ 2a — 2a?) + =f(b?+a,a). (39) 
For 1=|A—«) =2, (39) becomes 

a+k’ cdc 1 a+k’ cdc 

=(f +f cds’ a) (40) 


a—1 


as before. For — 


= (x /4a*) |[4a?— (6°?+a+2a) —In (b?+a—2a) ]+4a(b?+a) } =g(b?+<a, a). (41) 
Thus 


dkd2 
E=Akn? f f(b?+-a, a)+2(b?-++a, a) ], (42) 
(a?+a?/8)? 


where f is used for 0=a=2 and g for 2=a. 
The integration with respect to 2 must be carried out outside the sphere O’ of Fig. 2. Using elliptic 
coordinates a and X again gives 


l4+« 
x“ ada a+« 
b?,a)\dd}. (43 
+f g(b*, a) (43) 


f 
) 
| 
l 
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Now let 
f(b2, F(x, 2), f g(b2, a)\dd=G(«, d). (44) 
Then 
ada 
E=2rAk,* AS. 1)] 


The integration with respect to x is next carried out within the unit sphere. 


1 ada 1 
f a+x)— F(x, 1)] 


2 eda 1 
F(x, x) — 


Fla. xk) — F(x,a—k 
+f [F(«, a+%)— a—*)] 


1 
xdx[ G(x, x) —G(x,a—x)]}. (46 


Finally let . 
a) = [CF — F(x, 1) Jxdx, 


F2(x, a) = a+x)— F(x, a—x) ]xdx, (47) 


Gi(x, a) = ftv a+x)—G(x, a—x) 
which gives for the energy 


1 ada 
E=8nr°Ak,,* (1, —F, 


[Fi(1, a)+F.(a—1, a) — Fi(a—1, i(1, a) (48) 


3 ada 
+f ——. 
1 (a?+a?/8)? 2 (a?+a*/8)? 


since F,(0, a)=0 and G,(0, a) =0. Now it turns out that 
F\(1,a)—Fi(1—a, (49) 
so that the energy becomes 


The integrals of (47) have been evaluated by elementary integration and depend on a and a. The 


_ integrals were divided into two parts, those containing logarithms and those containing polynomials 
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only. The latter were integrated except for certain terms which it was convenient to combine with 
the logarithmic integrals to prevent their becoming too large to integrate graphically. This gives 
E in the form 


da 
60 4a‘(a*+a?/8)* 
+ (2a?+4a + a)*[16a? — 12a(a+2a*) + (a+ 2a*)*] In (2a*+4a+a) 
12a(a—2a*) +(a—2a?)?] In (—2a?+4a+a) 
+20a*(a?+2a+a)*(—a?+6a—a) In (a?+2a+a) 
+20a?(—a?+2a+a)*(—a?—6a+a) In (—a?+2a+a)+80a‘a' In a} 

+ { —256 —16a+(3v2 — 3)a?+(a?+32)-"( Sat —8a?+232a?— 1280a+1536) 
+169Vv2 tan-! (4v2/a)+(232+3a2/2) In (a?+32) /a?} 


E=A 


| In (4a+a) 


2a*+a)*{ (2a*+a)?— 80a? | In [1 — 16a?/(2a?+ a)? 


1+4a/(2a?+a) 
In +128a! 

1—4a/(2a?+a) 
—tan-! (4v2 ‘a) ]+32/a 


—32(a?+32)—"[ —3a—32+(160/a)]} }. (51) 


If in (51) a=0, then the value of the energy as given by (15) for metals may be evaluated. The 
integration may be carried out and gives 


E=Akp}(x*/3) { —32[1—(3) In 2] lim In a—(32/5)(In 2)*— (448/15) In 2+4n*—204/5}, (52) 


which is logarithmically infinite as has been mentioned in Section 2. 

The integrals still remaining in (51) were evaluated graphically for »=h?k,?/2m, h®kn?/4m, and 
h?km?/8m. These values of » correspond to gaps between the first and second Brillouin zones equal 
to twice, once, and one-half of the energy width of the first zone, respectively. From (32) the corre- 
sponding values of a are 4, 2 and 1. The integrand for a=0 at a= 40 differed by only 12 percent from 
that for a=2, and was approaching it. The integral from 40 to « made up 40 percent of the result. 
Therefore the integrals from a=40 to © were computed for all values of a by setting a=0. 

The energy is thus made up of a series of terms :; 


E = f f f terms) +| f f terms)| (53) 
0 2 40 0 2 


= — 366.36 — 17.58 +539.00 +399. 12) (54) 


for a=2. From (32), A = —162°e'mV/h?(2x)*= V/167' if the energy is to be in Rydberg units. If r, 
is defined by V = (4/3)xr,'2n where 2n is the number of electrons, then &,,’=92/4r,? and 


E= — (4/3) (ar3/162") (92 /4r,*) 290.748 = —0.0142 Ry (55) 
per electron. 
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Thus the following values of the energy have been obtained : 


E=—0.0304 Ry 
a=2, E=—0.0142 Ry (56) 
a=4, E=-—0.0073 Ry. 


The correlation energy as given in (35) depends only on a=8mn/h*k,? rather than on both » 
and k,, since V which occurs in A is inversely proportional to &,,’. This shows that (35) cannot be 
quite correct, because the correlation energy will be independent of » when 7 is very small, but will 
still depend on k». This is, of course, the limiting case of metals where, as shown in Section 2, the 
method breaks down. Again if k,, is very small, i.e., 7, very large, the correlation energy should 
approach zero, because the whole electrostatic energy (the partial absence of which is represented 
by the correlation energy) goes to zero. This is not manifest in (35) or (54). 

Wigner? gives —0.58/(r.+5.1) Ry as the correlation energy. If r,=4 Bohr units is taken as a 
reasonable value, this gives — 0.064 Ry so that the values in (56) are of the same order of magnitude. 
The values in (56) are nevertheless somewhat smaller than the —0.045 Ry estimated by Wigner 
to be correct. This is reasonable because one would expect the correlation energy to increase further 
with decreasing » even if it did not increase as drastically (to infinity) as (52) for this limiting case 
indicates. 

I am grateful to Professor Wigner for suggesting this problem and for his assistance in solving it. 
I also wish to thank the Physics Department of Princeton University for their hospitality and 


Allegheny College for a leave of absence. 
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A Determination of the Elastic Modulus s,; of Beta-Quartz 


A. W. Lawson 
Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 
(Received January 23, 1941) 


The adiabatic Young’s modulus of beta-quartz along a direction at 45 degrees to the optic axis 
has been measured as a function of temperature. These data, in conjunction with that of 
Atanasoff and Kammer, and of Perrier and de Mandrot, yield a value for the adiabatic elastic 
modulus s;; of beta-quartz equal to —0.226X 10" cm*/dyne at 600°C and indicate that the 
value reported by Osterberg and Cookson for c4, is inconsistent with the stability of beta-quartz. 


HE crystal structure of beta-quartz is The fact that the piezoelectric modulus d,, is not 
hexagonal holoaxial."* In virtue of the evanescent has been established by Osterberg and 
particular symmetry characteristic of such a Cookson‘ who succeeded in exciting shear vibra- 
structure, beta-quartz may be piezoelectric. The tions in beta-quartz by applying an alternating 
possible nonvanishing components of the piezo-electric field parallel to an axis of diagonal 
electric tensor expressed in Voigt’s* notation are: symmetry (X axis). From a knowledge of the 


000 di 0 O fundamental frequency of free vibration of a 
0000 -dy 0 thin crystal plate and the dimensions and density 
000 0 0 oO of the plate, Osterberg and Cookson, employing 


TW. Bragg and R. E. Gibbs, Proc. Roy. Soc. London # theory analogous to that developed by Mason,* 


A109, 414 (1925). Se 
2R. W. ¢. Wockoff, Am. J. Sci. 11, 112 (1926). 4H. Osterberg and J. W. Cookson, J. Frank. Inst. 220, 


3W. Voigt, Lehrbuch Kristallphysik (Teubner, 361 (1935). 
Leipzig, 1928), pp. 830-831. 5 W. P. Mason, Bell Sys. Tech. J. 13, 446 (1934). 
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calculated a value for the elastic constant C4, of 
beta-quartz. 

Recently, Atanasoff and Kammer® have also 
reported a value for cq. This value is based on 
the observations of higher harmonics of the 
fundamental frequency of shear vibration in X 
and Y cut plates. The calculation assumes that 
the eigenfrequencies of a finite plate will differ 
inappreciably from those for an infinite plate 
provided that the wave-length is small compared 
to the lateral dimensions of the plate. 

Unfortunately, the value of cy reported by 
Atanasoff and Kammer differs from that re- 
ported by Ostérberg and Cookson by nearly a 
factor of two. In discussing this discrepancy 
Atanasoff and Kammer comment on the in- 
adequacy of the mathematical methods used by 
Osterberg and Cookson and also on the possi- 
bility of an error in the identification of the 
modes of vibration. 

The object of this paper is to report some 
measurements of Young’s modulus in quartz 
along a direction making an angle of 45° with 
the axis of hexagonal symmetry (Z axis) and in a 
plane containing the axes of hexagonal and 
diagonal symmetry (XZ plane). These data in 
conjunction with the data of Perrier and de 
Mandrot’ prove that the value of cs, reported by 
Osterberg and Cookson is inconsistent with the 
stability of beta-quartz. Furthermore, the value 
reported by Atanasoff and Kammer is found to 
be consistent with the stability of beta-quartz, 
and is used to calculate a value for the elastic 
modulus $13- 

The data herein recorded are actually a by- 
product of an attempt made about a year ago 
to excite longitudinal vibrations in beta-quartz. 
Consider a set of axes X’ Y’ Z’ oriented in an 
arbitrary fashion with respect to the principal 
crystallographic axes X Y Z. Let ai, a2, a3, be 
the direction cosines of X’, 81, 62, 83, the direction 
cosines of Y’, and 71, Y2, ys, the direction cosines 
of Z’, referred to the old axes X Y Z. If now the 
possibility of exciting longitudinal vibrations 
along the Z’ direction by means of an alternating 
electric field along the Y’ direction be investi- 


( ml Atanasoff and E. Kammer, Phys. Rev. 59, 97 
1941). 

7A. Perrier and R. de Mandrot, Arch. des Sciences 4, 
367 (1920). 


gated, it is necessary to compute the value of the 
piezoelectric modulus d’:;. The value of this 
modulus in beta-quartz is given by :* 


For convenience choose 62=1, 8; =83;=0. Then 
the effective modulus will be a maximum if 
1=73=0.707, or, in other words, if the new 
primed axes are obtained from the old by a 
rotation of 45° about the Y axis. The value of 
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Fic. 1. The observed frequency of free longitudinal 
vibration of a (90°, 45°) cut quartz crystal as a function 
of the temperature. 


d's; in this case is —d,,/2. Since dy, has a finite 
value in beta-quartz, it is possible to excite 
longitudinal vibrations by an alternating electric 
field perpendicular to a hexagonal axis of sym- 
metry in a crystal rod whose length is parallel 
to a direction making an angle of 45° with the 
axis of hexagonal symmetry. Moreover, when 
such a rod is in the alpha-phase it is still possible 
to excite longitudinal vibrations since the value 
of d’s3 is given by :* 


23 = Bo( = —dia/2. (2) 


In practice, a crystal rod 3.756 cm X 0.2870 cm 
X 0.6275 cm is used. The smallest dimension is 
along the Y axis. Platinum electrodes are burned 
on the faces of the crystal perpendicular to the 


§ W. Voigt, reference 3, p. 849. 
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Y axis. A mount, constructed of stainless steel 
and Lavite, clamps the crystal lightly at the 
middle and simultaneously supports the crystal 
in the furnace and makes electrical contact with 
the electrodes. 

The power is supplied to the crystal through 
two amplifier buffer stages from a transitron 
oscillator whose frequency may be varied con- 
tinuously by a universal vernier condenser. The 
resonant frequency of the crystal is detected by 
the anomalous variation of the impedance of the 
crystal in the neighborhood of resonance as 
determined by a peak voltmeter. The frequencies 
corresponding to the maximum and minimum 
impedance of the crystal are measured by beating 
the unknown frequency against the standard sub- 
frequencies from a calibrated quartz-controlled 
oscillator driving two stages of multivibrators. 
The resonant frequency may then be calculated.® 

The temperature of the crystal was determined 
by means of a Chromel-Alumel thermocouple 
used in conjunction with a student potentiometer. 
It is estimated that the temperature is ascer- 
tained in this manner to within +2°C. 

Figure 1 shows the observed resonant fre- 
quency as a function of temperature. Values 
corresponding to the points plotted on the curve 
are tabulated in Table I. Measurements are not 
taken above 835°C for fear of breaking the crystal 
by transforming the quartz to trydimite. Obser- 
vations are taken with both rising and falling 
temperature and it is found necessary to stabilize 
the temperature for as long as half an hour in the 
neighborhood of the alpha-beta transformation to 
avoid the occurrence of an apparent temperature 
hysteresis. The alpha-beta transformation was 
successfully traversed at least six times. 

If the crystal rod employed were infinitely 
thin in cross section the resonant frequency at 
any temperature and the length / would be 
related to 1/s’33, the Young’s modulus of the 
crystal along the length, by the simple relation: 


1/s'33 = 4f*l*p, (3) 


where p is the density of quartz at the tempera- 
ture in question. Unfortunately, the finite cross 
section complicates the calculation somewhat. 
The effect of the finite cross section is to load the 


*S. Siegel and S. L. Quimby, Phys. Rev. 49, 663 (1936). 


crystal oscillator and depress the frequency. It is, 
however, possible to correct for this depression 
of frequency by calculating the kinetic energy of 
the crystal associated with motion in the lateral 
direction on the assumption that the mode of 
vibration is unaltered. The consequent change in 
frequency may then be calculated by using 
Rayleigh’s principle. The assumption that the 
mode of vibration is unaltered is justified if the 
fractional change of frequency is small. The 
result of such a calculation is: . 


1/s'ss=4 p(1+K), (4) 


where K is the term representing Rayleigh’s 
correction and is given by: 


where x’, v’, and 2’ represent the dimensions of 
the beta-quartz crystal and s’;3; and s’s3 are 
elastic moduli which determine the lateral con- 
tractions in the X’ and Y’ directions respectively 
for an extension in the Z’ direction. The values of 
the quantities s’;; and s’23 in terms of the 
principal elastic moduli are most easily found by 
cyclic permutation of the indices in the equations 
given by Voigt.'® Thus: 


(6) 


For a (90°, 45°) cut quartz crystal, Eqs. (6) and 
(7) reduce to: 


S’13=43(Sis +512), (8) 
(Sir +533 +2513 —S44) =S'33— 3544, (9) 


respectively. The above equations hold only 
when the crystal is beta-quartz. The expressions 
for alpha-quartz are quite analogous but some- 
what more complicated due to the lower sym- 
metry involved. In order to calculate the value 
of K it is necessary to know 5'33, 511, S33, $13, S12, 
and S44. 51, and S33 have been determined by 
Perrier and de Mandrot. sa is known from 
Osterberg and Cookson and Atanasoff and 
Kammer. Since K<1, it suffices to calculate an 
approximate value of s’3; from Eq. (3). At this 


10 W. Voigt, reference 3, p. 593. 
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juncture we anticipate a later calculation to 
remark that an approximate knowledge of s’3; in 
conjunction with the known values of si, S44, 
and $33, permits the calculation of an approxi- 
mate value for si; and the limits between which 
the value of si2 must lie. 

At 600°C the value of K, determined by 
averaging the values corresponding to maximum 
and minimum possible values of 512, is found to 
be 0.003. This value is certainly small compared 
to unity and thus justifies the calculation. 

The value of p at 600°C is 2.518 g/cm*." The 
length of the crystal at that temperature, calcu- 
lated from the thermal expansion data of Jay” 
is 3.808 cm. The observed frequency of free 
longitudinal vibration is 79.98 kc. These data, in 
conjunction with the value of K quoted above, 
yield a value for s’3; at 600°C. We find, in fact, 
from Eq. (4): 


$’33= 1.067 X cm?2/dyne. (10) 


To obtain further information concerning the 
elastic moduli of beta-quartz it is necessary to 
relate this quantity to the principal elastic 
moduli. The simplest relation is according to 
Voigt 
S’33=S11(1 — ys")? 
(11) 
which in the case of a (90°, 45°) cut reduces to: 
$'33= (Sur +S33 +544 +2513). (12) 
Since the compressibility x is given by: 
X = +2513) +533, (13) 
we may write Eq. (12) in the form: 


$'33=4(x+533 — 2512+ 2544). (14) 


In order to calculate the value of sis, it is 
necessary to know the values of 541, S33, and S44. 
By linear interpolation of the data of Perrier 
and de Mandrot the values of the isothermal 


1 A. L. Day, R. B. Sosman and J. C. Hostetter, Am. J. 
Sci. 37, 16 (1934). 

"A. H. Jay, Proc. Roy. Soc. London A142, 237 (1933). 

18 W. Voigt, reference 3, p. 740. 


moduli s;;'*° and s33'*° at 600°C are found to be: 


=0.9124 cm?/dyne, (15) 
= 1.029 K cm?/dyne. (16) 


In order to calculate the corresponding adiabatic 
values, use must be made of the relation: 


= T Vana / Cp. (17) 


By linear interpolation of Jay’s data on the 
thermal expansion of quartz, the values of a, 
and a; are found to be: 


+21.4X10-*/°C, (18) 
a3= —6.9X10-*/°C, (19) 


According to Kelley,“ the value of C, may be 
calculated from the equation: 


C,=10.95+0.0055T cal./mole—°C. (20) 
From Eq. (17) we find: 

=5,,=0.898 cm?/dyne, (21) 

=S33 = 1.027 cm?/dyne. (22) 
Using the value s4,=5.165 X10-" cm*/dyne re- 


ported by Osterberg and Cookson and Eqs. (14), 
(21), (22) and (10) we find: 


2512 = (x +2.821X10-'*) cm?/dyne (23) 
and, since x >0, this requires that: 
Sig > +1.411 cm?/dyne. (24) 


However, since s¢g=2(si1—S12) and is a positive 
quantity, we must find: 


S12 < +0.898 X cm?*/dyne. (25) 


Since the quantity s;, obviously cannot satisfy 
conditions (24) and (25) simultaneously, we are 
forced to conclude that 51, S33, Sis OF S44 is in 
error. In view of the discrepancy between the 
values of s44 reported by Osterberg and Cookson 
and Atanasoff and Kammer, s44 seems the most 
likely suspect. Indeed, using s44=2.795X10-" 
cm?*/dyne as reported by the latter, we are led 
to the conclusion that: 


2512 > — 1.919 10-"* cm?/dyne, (26) 


4K. K. Kelley, Bur. Mines Bull. 371 (1934). 
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TABLE I. The observed frequency f of the (90°, 45°) cut 
quartz crystal as a function of the temperature T. Frequencies 
are given in kilocycles and temperatures in degrees centigrade. 


f r f T f 

25 82.43 | 482 77.32} 569 71.23 | 660 80.42 
175 81.25 | 500 76.65] 574 66.00 | 710 80.56 
185 81.13 | 513 76.10} 579 74.46 | 765 80.65 
275 80.36 | 536 75.13} 580 76.25 | 795 80.79 
357 79.37 |545 74.33] 588 79.64 | 825 80.91 
405 78.68 | 552 73.70; 595 79.92 | 835 80.95 


466 77.47 | 557 73.00| 620 80.18 


a condition which is quite compatible with the 
conditions for stability of the lattice. Further- 
more, by means of Eq. (13) we arrive at the 
value: 


$13 = —0.226X10-" cm?/dyne. (27) 


Perrier and de Mandrot also measured Young’s 
modulus in a (0°, 50°) cut quartz crystal. From a 
somewhat questionable linear interpolation of 
their data, we find at 600°C: 


E(0°, 50°) =0.951 X10! dyne/em?, (28) 


which, in conjunction with Eqs. (11), (17), (21) 
and (22), and the value of sy, reported by 
Atanasoff and Kammer, leads to a value: 


$13 = — 0.235 X10-" cm?/dyne. (29) 


The agreement between these two values (27) 
and (29) is certainly within the experimental 
error and tends to confirm the conclusion that 
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the result reported by Osterberg and Cookson 
for C44 is incorrect. 

Due to the numerous sources of data employed 
in the calculation, an actual estimate of the 
accuracy of the result here reported for si3 is 
difficult. It is doubtful that the precision measure 
is less than about +10 percent, inasmuch as the 
quantity in question is a small difference between 
large quantities. 

If, indeed, Osterberg and Cookson correctly 
identified the modes of vibration which were 
studied by them, then the conclusion must be 
made that their method for obtaining the solu- 
tion to the differential equation governing the 
propagation of elastic waves is unsuited to the 
calculation of the eigenfrequencies of a thin 
plate. 

It may be of interest to note that the par- 
ticular type of piezoelectric drive here described 
should be the most efficient manner of exciting 
longitudinal vibrations in all piezoelectric crys- 
tals belonging to the symmetry classes D4, C,, 
Va, and 73. 

In conclusion, the author desires to express his 
gratitude to the Faculty Research Committee 
of the University of Pennsylvania for a grant 
used in connection with this experiment, to the 
members of the Works Progress Administration 
project number 25796 for much routine assist- 
ance, to Dr. G. P. Harnwell, Director of the 
Randal Morgan Laboratory of Physics, for his 
kind interest, and to Dr. Frederick Seitz for 
several stimulating discussions. 
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Letters to the Editor 


PHYSICAL REVIEW 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
lo this department. Closing dates for this department are, for 
the first issue of the month, the eighteenth of the preceding 
month, for the second issue, the third of the month. Because 
of the late closing dates for the section no proof can be shown 
to authors. The Board of Editors does not hold itself responsible 
for the opinions expressed by the correspondents. 
Communications should not in general exceed 600 words 
in length. 


Vibrational Structure of the *X,+—'Z,*+ 
Rydberg Series of N; 
Y. TANAKA AND T. TAKAMINE 


Department of Physics, University of California, Berkeley, California 
March 14, 1941 


N 1938, Worley and Jenkins' found a new Rydberg 
series in N2 in absorption, namely the (0,0) bands of a 
series of electronic states converging at *Z,* of N.* 
(Worley-Jenkins series). Since there are a considerable 
number of other absorption bands of nitrogen which are 
not yet analyzed,? especially below \800A, a further study 
was undertaken in Berkeley with a 3-meter grazing- 
incidence vacuum spectrograph. The helium continuum 
was used to produce the background, the pressure of 
nitrogen in the spectrograph being about 1/100 mm with 
the path length 130 cm. 

One important outcome was the finding of (1,0) bands 
of the above series, forming another Rydberg series. 
Naturally, the newly found (1,0) bands appear with much 
less intensity than the (0,0) bands, but their existence 
seems to be beyond doubt and constitutes a clear evidence 
as regards the reality of the Worley-Jenkins series. 

In the region between \938A and 784A, corresponding 
to nm=2 and n=17, altogether twelve members were 
observed, four members n=5, 7, 8 and 13 being overlaid 
by other strong bands. With the data given in the above- 
cited report by Worley for the Worley-Jenkins series, 
values of AG; (in cm™) were calculated for different elec- 
tronic states and found to agree within a few units. 

The mean value of AG; of eleven members observed 
comes out as 2171 cm™ and agrees closely with that for 
N2* *E,* state, namely 2174.75 cm™. The limit of the (1,0) 
bands lies at 127,834 cm™ (A782.26A) or at 15.77 volts, 
while the addition of AG; for N,* *2,* to the reported limit 
of the Worley-Jenkins band gives 127,845 cm™. 

A more detailed report of this experiment, including the 
confirmation of the Nz emission bands* which converge to 
the same limit as that of Hopfield’s Rydberg series (con- 
verging at *Z,,*), as well as the results for Hz, O2 and CO 
absorption bands, will be published in the Scientific Papers 
of the Institute of Physical and Chemical Research (Tokyo). 

It is a pleasure to thank Professor F. A. Jenkins for his 
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continued interest and advice, and to express our gratitude 
to Professor R. T. Birge for allowing us to work in Le Conte 
Hall. 


1 R. E. Worley and F. A. Jenkins, Phys. Rev. 54, 305 (1938). 

2 In the region between 1000A and 800A, many bands were recently 
anal by Dr. R. E. Worley in Berkeley. 

3T. Takamine, T. Suga and Y. Tanaka, Sci. Pap. Inst. Phys. Chem. 
Research (Tokyo) 34, 854 (1938). 


Mean Life of Slow Mesotrons 


FRANCO RASETTI 
Department of Physics, Laval University, Quebec, Canada 
March 15, 1941 


HE author recently described! an experiment which 
enabled him to detect the disintegration electrons 
emitted by mesotrons at the end of their range. 

This is a preliminary report of new results obtained 
with an improved arrangement which was designed to 
measure a decay curve of mesotrons at rest. 

The set-up may be briefly described as follows. A fourfold 
coincidence counter set defines a beam of mesotrons, which 
impinges upon a block of iron 10 cm thick and 2.5 cm wide. 
A battery of anticoincidence counters, placed below the 
iron, selects the events in which a mesotron is absorbed. 
It is found that a fraction of these absorption processes is 
associated with the emission of a particle from the ab- 
sorber. This particle may be a disintegration electron, or 
possibly a scattered mesotron. 

A system of circuits simultaneously records the numbers 
n,, M2 and ns of such particles emitted within, respectively, 
36, 3.1 and 1.2 microseconds after the passage of the 
primary mesotron. It is found that many of the particles 
are delayed, which is expected to be the case for the disin- 
tegration electrons, but not, of course, for scattered 
mesotrons. 

The conditions of the experiment are such that m, m2 
and ms are affected by a “background” due to undelayed 
processes (scattered mesotrons+showers) which, however, 
is identical (not statistically) for all three; hence the dif- 
ferences are significant. Since we can safely assume that all 
mesotrons have decayed within 36 microseconds, we may 
write: 

nz) =exp (—1.9X10~*/r). 

The results for 207 hours of counting are m,=170, 

n2= 142, n3=119. This yields: 


7=3.1+41.5 microseconds. 


So far the accuracy of the present measurement is 
rather poor; its interest lies rather in affording a deter- 
mination of the mean life that is more direct and less 
dependent upon accessory hypotheses than the one 
deduced from the atmospheric absorption effect. 

The measurements are being continued in order to 
increase the statistical accuracy of the result. An inves- 
tigation of the important point—whether all slow meso- 
trons or only the positive ones emit a disintegration elec- 
tron*—is also being planned. 


1F. Rasetti, Phys. Rev., in press. 
2S. Tomonaga and G. Araki, Phys. Rev. 58, 90 (1940). 
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On the Fine Structure Pattern of Cosmic 
Rays at Mexico City. II. 


E. J. Scuremp, Washington University, St. Louis, Missouri 
AND 


ALFREDO BANos, Jr., Institute of Physics, National University of Mexico, 
Mexico City, Mexico 


March 15, 1941 


HE survey of directional cosmic-ray intensity at 
Mexico City which was reported in a previous com- 
munication! has since been extended to include the four 
additional azimuths a=45°, 135°, 225°, 315° (NE, SE, 
SW, NW). The resulting zenith angle intensity distribu- 
tions J(z,a), together with those previously found, are 
shown in Fig. 1. From these eight curves it is now possible 


Fic. 1. Normalized directional cosmic-ray intensity J(s, a) at Mexico 
City, as a function of zenith angle sz, for the eight azimuths S-N (top 
curve); SE, NW (second curve); E-W (third curve); and NE, SW 
(bottom curve). Each vertical scale division represents 25 percent of 
the vertical intensity. The average probable error is about 1 percent 
of the vertical intensity. 


to deduce the principal features of the intensity pattern 
I(s, a), and to infer therefrom some interesting properties 
of the primary radiation. 

First, the W and SW azimuths, together with their next 
neighbors, show the existence of an absolute maximum of 
I(s, a) at about z=12°, a=250° (20° S of W). Similarly, 
the E and SE azimuths and their neighbors reveal a second 
prominence at about z= 15°, a=110° (20° S of E). Unlike 
the prominences observed in Missouri,? these two are 
isolated from one another, and do not join to form a 
continuous ring of prominences concentric about the 
zenith. With respect to position, they are very nearly 
mirror images of one another in the N-S plane. Their 
E-W positional symmetry is of the kind which has previ- 
ously been considered by one of us,’ and by reason of which 
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it was concluded that the primary radiation contains 
positive and negative particles of the same e/m. The 
relative magnitudes of the western and eastern prominences 
are about 130 percent and 95 percent of the zenith in- 
tensity, respectively. If each prominence were due to 
primaries of but one sign of charge, these intensity values 
would give the relative abundance of the positive and 
negative primary particles concerned. The ratio 130/95 
must be considered a lower bound for the positive-negative 
ratio of abundance, which actually may be appreciably 
greater. 

A second observation, which correlates well with the 
preceding one, is that the intensity of the N-S plane falls 
off less rapidly at large zenith angles than the intensity in 
adjacent azimuths. This fact indicates that there is but a 
slight overlapping of the allowed cones for positive and 
negative particles, and that this overlapping occurs in the 
limited region to either side of the N-S plane, thereby 
producing the observed excess intensity in that region. It 
appears, therefore, that the spectra of positive and negative 
primaries both possess a band of low energies which falls 
off rapidly within an energy range in the vicinity of 0.42 
stérmer, or about 11 billion ev for electrons, corresponding 
to an allowed cone which, for positive particles, opens just 
to the east, and for negative particles just to the west, of 
the N-S plane. This observation does not, of course, 
exclude the existence of lines or bands in the primary 
spectra at still higher energies. 

That such additional higher energy lines or bands do 
exist is indicated by the further observation that, in all 
the eight curves of Fig. 1 excepting the southern curve, 
there is evidence for a circular ring of prominences at 
z=42°, concentric with the zenith. This type of locus of 
prominences is similar to those observed in Missouri,? 
which have been explained in terms of atmospheric ab- 
sorption of primary spectral lines or bands too energetic 
to be influenced by the earth’s magnetic field.‘ On the 
basis of equal atmospheric path length, the position z = 42° 
of the ring at Mexico City corresponds quite closely to the 
position z=7° of the innermost ring observed in Missouri 
when account is taken of the difference in altitude of the 
two stations. In the case of the 42° ring at Mexico City 
and the 7°, 20° and 35° rings in Missouri, all of which are 
insensitive to the earth’s magnetic field, there appears, so 
far, no way of distinguishing the sign of charge of the 
associated primary particles. 

The present survey is now being further extended in 
order to improve and interpolate the data shown here. 
We wish again to express our thanks for support from 
the Penrose Fund of the American Philosophical Society, 
the Rumford and Permanent Science Funds of the 
American Academy of Arts and Sciences, and from 
Washington University and the National University of 
Mexico. 


1E, J. Schremp and Alfredo Bafios, Jr., 4g Rev. 58, 662 (1940). 
In this reference ¢ designates azimuth angle, but henceforth a will be 
so used, in order to coincide in nomenclature with the Lemaitre- 
Vallarta theory of the allowed cone. 

2D. M. Cooper, Phys. Rev. 58, 288 (1940). 

3E. J. Schremp and H. S. Ribner, Rev. Mod. Phys. 11, 151 (1939). 

4 Reference 3, p. 150. 
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LETTERS TO 


The Nature of the Primary Cosmic Radiation 
and the Origin of the Mesotron 


MARCEL SCHEIN, WILLIAM P. JESSE AND E. O. WOLLAN 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
March 13, 1941 


URING the past year, our counter measurements of 
the vertical intensity and the production of meso- 
trons at high altitudes have been continued.' In these 
measurements various arrangements of counters in three-, 
four- or fivefold coincidences have been used. The vertical 
intensity has been obtained with lead thicknesses of 4, 6, 
8, 10, 12 and 18 cm interposed between the counter tubes. 
The combined results of these experiments are plotted as 
Curve A of Fig. 1. It is seen from these measurements that 
the intensity of the hard component increases continuously 
to the highest altitudes reached. In our previous paper a 
single point of low statistical weight indicated a maximum 
in the mesotron curve which is not confirmed in our sub- 
sequent experiments. In the original experiments with 8 
and 10 cm of lead we considered the possibility of a con- 
tribution at very high altitudes to the observed intensity 
from primary electrons with sufficient energy to penetrate 
the lead (E>10"* ev). This consideration led us to carry 
out the experiments with greater and with smaller lead 
thicknesses. The close agreement between the points 
obtained with the various thicknesses to pressures of 2 cm 
Hg (less than 1 radiation unit from the top of the atmos- 
phere) is evident from the figure. This shows that measure- 
ments of the hard component made at very high altitudes 
with lead thicknesses even as small as 4 cm are not appre- 
ciably affected by electrons. 

Further evidence that the traversing particles are not 
electrons was obtained by an arrangement of side counters 
registering showers generated in the lead by the traversing 
particles. A typical arrangement is shown in the figure in 
which counters 1, 2, 3, 4 and 2, 3, 4, 5 register the vertical 
intensities for 4 and 6 cm of lead, respectively, and counters 
1, 2, 6, 4 and 2, 6, 4, 5 register particles accompanied by 
showers. If the traversing particles are electrons, there is 
a high probability of generating in the first 2 cm of lead 
a shower of many particles. In no case were more than a 
few percent of the traversing particles accompanied by 
shower counts in the side counters. 

Because of the constancy of the penetrating power of the 
particles which we measure, and because they are not 
shower producing, we conclude that there are no electrons 
of energies between 10° and 10" ev present at the highest 
altitudes reached. Since the energy required for electrons 
to penetrate the earth’s magnetic field of 51° geomagnetic 
latitude is about 310° ev, and since our measurements 
were carried out to within the first radiation unit from the 
top of the atmosphere, it seems difficult to assume the 
presence of electrons (E<10" ev) in the primary cosmic 
radiation and, hence, they must be replaced by some 
penetrating type of charged particles. The mesotrons 
themselves cannot be the primaries because of their 
spontaneous disintegration. Hence, it is probable that the 
incoming cosmic radiation consists of protons. The fol- 
lowing facts support this assumption. 


THE EDITOR 615 


Lead 

4 4¢m 

Pfotzer 
v wo 

o l2" 

x 18 


Fic.1. Curve A: Intensity of the hard component for various lead 
as a function of pressure in cm Hg. Curve B: Total vertical 
intensity of cosmic rays obtained by Pfotzer as a function of pressure. 


1. Another experiment which we have performed at 
high altitudes shows that mesotrons which can penetrate 
18 cm of lead are produced in multiples mainly by ionizing 
nonshower producing particles.* 

2. The number of incident particles as determined by 
Bowen, Millikan and Neher® from their ionization chamber 
measurements at high altitudes is approximately the same 
as the number of penetrating particles which we observe 
close to the top of the atmosphere. 

3. The measurements of the east-west asymmetry of 
cosmic rays have led Johnson‘ to suggest that the primaries 
of the hard component are probably protons. (In order to 
make it certain that all the incoming cosmic rays are 
positively charged, an east-west experiment for pen¢- 
trating particles should be carried out at high altitudes.) 

We hope to continue these observations with lower 
thicknesses of lead to compare with the measurements 
without absorption screen made by Pfotzer. Since we have 
assumed that the primary cosmic radiation consists of pro- 
tons, the electrons known to exist in large numbers in air 
at high altitudes (curve B) must be of secondary 
origin. Furthermore, as seen from our experiments, 
the average energy of these electrons is low (E<10* ev). 
These facts suggest that the electrons in the atmosphere 
arise mainly from the decay of the mesotron and knock-on 
processes. 

The writers wish to express to Professor A. H. Compton 
their appreciation for his support of these experiments and 
his continued interest in them. 


1M. Schein, W. P. Jesse and E. O. Wollan, Phys. Rev. 57, 847 


(1940). 
? This proaees is in otfiticn to the already reported of the 


process 
of mesotrons non-ionizing radiation. 
assay S. Bowen, R. A. Millikan and H. V. Neher, Phys. Rev. 53, 217 


4T. H. Johnson, Rev. Mod. Phys. 11, 208 (1939). 
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Angular Momenta of Gravitational Fields 


MARIO SCHOENBERG* 
University of S. Paulo, Brasil, South America 
March 5, 1941 


XPRESSIONS for the components of the angular 

momentum density of the gravitational field can be 

obtained by taking into account the invariance of the field 
equations under the infinitesimal linear transformation: 


= +x, (1) 
the w," satisfying the conditions 


which characterize an infinitesimal Lorentz transformation.' 
The calculations can be carried out by expressing the 
invariance of the Lagrangian: 


1 


the I being the Christoffel symbols of the second kind. It 
is convenient to consider & as a function of G*” and G,*": 


Starting from the relation: 
L(G"(x"), = Gat”(x)) (3) 
a straightforward calculation leads to the identity 


0 


av ar ) 
S Bap — py av). (5 


Ta’ and ta’ are, respectively, the components of the energy 
momentum density of matter and the gravitational field. 

As a result of (4) S3-@! and S*!? can be interpreted 
as the components of the spin density of the gravitational 
field and and 


= [eata®x? ept,®x*] Shap 


as the components of the total angular momentum density. 
This is not a tensor density but it behaves as a tensor 
density under linear transformations of coordinates of the 
Lorentz type. 

The space integrals M*? of the expressions 


e,(T,°+ t,*)x*]+E% 


do not depend on x° as can be seen from (4). With the use 
of Cartesian coordinates in the flat region surrounding the 
world tube of the matter and going over continually into 
arbitrary coordinates inside the tube, it can be seen that 
the values of the M*? do not depend on the choice of the 
coordinates in the interior of the tube so long as the Car- 
tesian coordinates outside remain the same; under a trans- 
formation of these Cartesian coordinates the M®? transform 
as the components of a tensor of the second order. 

In the case of special relativity the conservation of 
angular momentum results from the invariance under the 
group of rotations; in general relativity the same connection 
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holds if we consider the invariance under the rotation group 
of the flat space surrounding the world tube of matter. 

I thank Professors Gamow and Pauli for discussions. 

* Now at the George Washington University with a fellowship of the 
Guggenheim Foundation. 


1 This is a particular case of the general method discussed by Belin- 
fante, Physica 6, 887 (1939). 


On Time Lags in Coincident Discharges 
of Geiger-Miiller Counters 


M. E. Rose anv W. E. RAMSEY 


The Bartol Research Foundation of the Franklin Institute, 
Swarthmore, Pennsylvania 


March 12, 1941 


N specific ionization measurements with Geiger-Miiller 
counters it has been observed that the efficiency of 
oxygen counters is abnormally low, indicating about 30 
percent of the expected ionization by cosmic-ray particles." 
As an explanation of this anomalous result it has been 
suggested that, as a result of electron capture, time lags** 
are introduced between the passage of the ray through the 
counter and the initiation of the discharge. 

To test this hypothesis and to study these time lags we 
have used a vertical arrangement of three counters in coin- 
cidence with the two outer counters sufficiently small so 
as to define the path of the cosmic-ray particles which 
were used as the primary radiation. By moving these 
exploring counters relative to the central counter (see 
Fig. 1) the (local) efficiency of the latter could be measured, 
as a function of resolving time of the amplifier, for various 
distances between the point of origin of the electrons and 
the counter wire. The resolving time, during which the 
central counter must be discharged to produce a count, 
was measured with a pulse generator which supplied 
known pulses separated by known time intervals.2? The 
amplifier circuit, which was designed for great sensitivity, 
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Fic. 1. Counter efficiency as a function of resolving time (logarithmic 
scale) for oxygen and A —Oz (0.94 A, 0.06 O2). The curves i, #¢ and iii 
were obtained with the peers of the exploring counters as shown 
in the upper portion of the figure. Standard statistical errors are indi- 
cated. For the A —O: counter the crosses refer to position i and the 
circles to position it. The exploring counters were 1 cm in diameter, 
20 cm long; the central counter was 2.7 cm in diameter, 29 cm long. 
In each case the gas pressure was 9 cm Hg. 
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was subjected to the most stringent tests to insure that 
all spurious effects were eliminated. The runs were repeated 
in random order and the results were found to be fully 
reproducible. 

The measurements were made for two extreme cases; 
namely, for oxygen which has a comparatively large capture 
cross section and for an A—O, mixture (6 percent O2) which 
is known to be very efficient and for which very little 
capture effect may be expected. The efficiency of the O» 
counter (Fig. 1) is seen to be rather low and to exhibit a 
very pronounced time effect which is readily interpretable 
in terms of capture. At the pressure used (9 cm Hg) the 
total ionization of cosmic rays in the counter should be 
rather small' (2-5 ion pairs) and there is an appreciable 
chance that all electrons formed will be captured before 
they reach the wire. The time lags, between the discharges 
in the central and exploring counters, will be statistically 
distributed about a most probable value (~10-° sec.). 
The resonance-like form of the curves in Fig. 1 may then 
be understood. At small resolving times very few ions 
reach the vicinity of the wire in time for a coincident dis- 
charge whereas at long times all of them reach the wire in 
time so that the efficiency is then mainly determined by 
the probability for cosmic-ray ionization. The distribution 
in time lags is obviously given by differentiation of the 
curves. The effect of changing the path length for cosmic- 
ray ionization and for electron capture is shown by the 
three curves marked 7, ii and iti and the results are clearly 
in accord with expectations. 

Further confirmation of the capture hypothesis is ob- 
tained from more quantitative considerations. The prob- 
ability g that an electron will be captured before reaching 
the wire is given in terms of the efficiencies at low and high 
resolving times (¢, and ¢,, respectively) by 


qg=1—log (1/1—«,)/log (1/1—«,) (1) 
and the probability of capture per collision is 


1 

p=1-—(1-g) ~ ims’ (2) 
where N=total number of collisions. From curve ii one 
finds g=0.85 and p~10~. From the mean time lag 
(4X10-5 sec. from ii) and the mean distance traversed 
before capture (calculated from 2) the ionic mobility is 
found to be 2 cm?/sec./volt at atmospheric pressure. These 
results are in as good agreement with independent capture 
and mobility measurements as could be expected.‘ 

In contrast to the results for O. the A—O, mixture 
shows an efficiency (97 percent) which is independent of 
resolving time and relative position of the counters. 

In the interest of higher efficiency it may, therefore, be 
suggested that the use of gases which capture electrons 
easily should be avoided. Otherwise long resolving times 
are required and care should be exercised in changing 
circuit constants. 

'W. E. Danforth and W. E. Ramsey, Phys. Rev. 49, 854 (1936) and 
unpublished results. 

2C. G. Montgomery, W. E. Ramsey, D. B. Cowie and D. D. Mont- 
nT. Du Phys. Rev. 56, 635 (1939). 

we Nature 144, 152 (1939). 


i B. Loeb, Fundamental Processes of Electrical Discharge in Gases 
(John Wiles, Sons, 1939). 
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Relative Importance of Different Elements 
for Neutrino Production 


G. Gamow 
The George Washington University, Washington, D. C. 
March 7, 1941 


N the article appearing in this issue, the author and 
M. Schoenberg propose the theory of neutrino energy 
losses in contracting stars, and develop the formula for 
the rate of these losses as a function of the temperature and 
density. 

In view of application to the theory of stellar collapse, 
it is interesting to investigate the relative importance of 
different elements for such processes, taking into account 
their actual abundance in stellar matter, and the known 
values of the decay energies and decay periods of the 
corresponding unstable nuclei. From the table of artificial 
8-active nuclei published by G. T. Seaborg,! and from the 
table of relative abundance of elements in the cosmos by 
V. M. Goldschmidt? the diagram shown in Fig. 1 was 
constructed. 
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Fic. 1. A diagram showing relative importance of different elements 
for the energy losses through neutrino emission. The fractions given 
on the abscissa axis represent the contraction of stellar radius [for 
peter (©)] for which the corresponding central temperatures are 
at 


For the sake of simplicity, each isotope has been charac- 
terized schematically by the temperature at which 
the neutrino-producing urca-process reaches saturation: 
T=Ewnax/k, and the rate of energy loss at that point: 
[W]=Emax'%-c (¢ being atomic concentration of the 
element in question). It must be remembered that actual 
energy loss begins at a somewhat lower temperature than 
that corresponding to saturation, and that the energy 
carried away by neutrinos continues to increase slowly for 
still higher temperatures. (See the exact curve for Fe**, 
as given in the article.) The curve for He*® was calculated 
for an assumed concentration of 10~, which follows from 
the equilibrium conditions between the He*-producing 
processes: (H'+H'—H?+ $+; H?+H*~He'+n) and He’- 
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destroying process: (He*+He*~Be’). It is probable that 
the concentration of He? is much lower. 

It is interesting to notice from the diagram how the 
rate of energy losses is affected by different factors. Thus, 
for example, urca-process in Fe gives, at 10® °C, very 
large energy losses, in spite of the comparatively small 
abundance of this particular iron isotope. On the other 
hand, the abundant element N*™ leads to a much smaller 
energy loss because the 8-decay of C™ belongs to the 
transitions prohibited by the selection rule; if this transition 
were permitted, the energy loss through neutrino emission 
would be as high as 10® or even 107 erg/g sec. It is thus 
possible that, between as yet uninvestigated §-active 
elements, one might find one or two which will give very 
high energy losses even at 10° °C. 

But, using only the tabulated nuclei, we find that energy 
losses as high as 10" or 10" erg/g sec. will take place at the 
temperature of 10!° °C, Such a temperature will be attained 
in the interior of a contracting star of 5 sun masses when 
its radius is reduced to one percent of its original value. 
It is quite possible that this stage of contraction cor- 
responds to the beginning of the catastrophic collapse 
characterizing the processes of stellar explosion. 


1G. T. Sea , Chem. Rev. 27, 199 (1940). 
2V. M. Goldschmidt, Skrift. Norsk. Vid. Akad. N. 4 (1937). 


The Absolute Cross Section for the Photo-Disin- 
tegration of Deuterium by 6.2-Mev Quanta 


James A. VAN ALLEN* AND NICHOLAS M. SMITH, JR.* 


Department of < Institution of 
ngton. 
a 1, 1941 


HE photo-disintegration of deuterium is of special 
interest because of the relative simplicity of its 
theoretical treatment. We have been engaged in deter- 
mining the absolute cross section of the photo-disintegra- 
tion with the 6.2-Mev radiation produced when fluorine is 
bombarded with protons. An accurate measurement is 
made possible by the recent calibration of the yield from 
a thick target of crystalline calcium fluoride! as 8.9+0.5 
X10 quanta per microcoulomb of protons at 370 kev and 
3.74+0.2<10* quanta per microcoulomb at 1050 kev. 

A high pressure ionization chamber filled with deuterium 
gas was used to count visually the pulses due to the 
photo-protons. The chamber consisted of a high voltage 
electrode cylindrical in shape with a grid electrode along 
its axis, the whole being placed inside a steel pressure 
vessel with auxiliary pressure gages and valves. The 
volume of the chamber was about 36 cc. 

Blank runs made with ordinary hydrogen in the 
ionization chamber gave null results, showing that neutrons 
produced in the magnetic deflection box of the accelerating 
tube by deuterium contamination in the beam were not 
present in sufficient numbers to make an appreciable 
background. 

With deuterium in the chamber the pulse counting rate 
was observed to parallel closely the quantum yield curve, 
both with respect to varying the incident beam current 
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at constant bombarding voltage, and to varying the bom- 
barding voltage. In order to increase the heights of the 
pulses the stopping power of the gas in the ionization 
chamber was increased by the addition of about ten percent 
of argon. Check runs taken with hydrogen and argon 
showed that no photo-disintegrations occurred with the 
mixture. 

Inasmuch as actual counting conditions required that 
the pulses exceed a certain minimum height to be detected, 
a fraction of the pulses was lost. These short pulses were 
those from disintegrations occurring near a wall of the 
ionization chamber with the proton heading into it. The 
apparent cross section, calculated assuming that no 
counts were lost, was measured at various pressures of the 
deuterium-argon mixture. The true cross section was then 
determined by extrapolating to infinite pressure (where the 
wall effect is obviously zero) the curve of the apparent 
cross section as a function of the reciprocal of the pressure. 

Two such curves with different geometrical conditions 
and with different gas mixtures have been taken during 
January and February, 1941. In one the axis of the chamber 
was perpendicular to the line from the target, and in the 
second the axis of the chamber passed through the target. 
The extrapolation in the first case gave 11.3 cm? 
and in the second case 12.3X10-** cm*. The weighted 
average of these two results taking all known errors into 
account is then 


11.6+1.5 X 10-*8 cm?. 


In these curves the wall effect at the highest pressure 
amounted to about 20 percent. We are preparing to repeat 
these experiments under conditions in which the wall effect 
is considerably less. The ultimate accuracy with our present 
experimental arrangement can be expected to be about ten 
percent. 

Our value is approximately equal to that of the Bethe- 
Peierls theory* for zero force range (12.6 1078 cm*) and 
is therefore distinctly lower than any present theory® 
predicts when based on a reasonable value for the force 
range. 

In preliminary experiments with the 17.5-Mev gamma- 
rays produced when lithium is bombarded with protons, 
no effect which is attributable to deuterium has yet been 
observed. The counting rate with the chamber filled with 
hydrogen has in fact exceeded several fold the expected 
rate from deuterium. This is not surprising since the y-ray 
energy is in excess of the binding energy of the last proton 
or neutron of nearly every nucleus. A carbon-walled 
ionization chamber, however, shows good prospects of 
success in our latest trials. 


ie Institution Fellow. 
an M. Smith, Jr., Foes. Rev. 59, 501-508 (1941). 
- A. Rev. 


Mod. Phys. 8, 124 (1936). 
. W. “aha O. Mohr, Proc. at A148, 
2066225 reit and E Condon, Phys. Rev. 49, 904-911 


Hf G. Breit, J. R. Stehn and E. 
1937); H. Fréhlich, W. Heitler and BO 
85-102 (19 40). Professor E. Teller to us last 
the results of the calculations of Fréhlich et al. were incompatible with 
the general theory of the interaction of the electromagnetic field with 
a nucleus. See A. J. Sent, Phys. Rev. 52, 787-789 (1937) and 
W. E. a: Jr. and L. I. Schiff, Ph wy 53, 651-661 (1938). 
Recently [W. Rarita, J. Yt. a H. A. Nye, Phys. Rev. 59 
Perr (1941)] a detailed recalculation of the meson theory has 
appeared; numerical results were given only for the cross section and 
angular distribution for 17.5-Mev radiation. 
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